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The Effect of Boundary Conditions and Mesh Size 
on the Accuracy of Finite Difference Solutions of 
Two-Point Boundary Problems 


By Davip BorweEIn and ANDREW R. MITCHELL, St. Andrews, Scotland i) 


1. Introduction 


In a recent communication, Fox and MItcHELt [1]?) showed that in certain 
problems, provided the finite difference interval is chosen with sufficient care, 
boundary value methods give good results when step-by-step solutions are 
unstable. The purpose of the present paper is to examine critically the accuracy 
of boundary value techniques based on finite difference methods as applied 
to the numerical solution of ordinary differential equations with two-point 
boundary conditions. 

Fox [2], using the method of matrix inversion, examined this problem in 
some detail. In the present paper, matrix methods are not used. Instead, exact 
solutions of the finite difference replacement of the chosen differential equation 
are obtained and compared with corresponding exact solutions of the differen- 
tial equation for a variety of mesh lengths and two-point boundary conditions. 


In addition, the accuracy of approximate numerical solutions of the differ-. 


ence equation is examined. The latter solutions can be obtained either by an 
iterative process such as relaxation or by direct methods of solution of the 
associated set of simultaneous equations. 
2. Conditions for Which the Difference Equation Has No Solution 
Consider the differential equation 
2 
eee ee 8) (1) 
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PD Davip BorWEIN and ANDREW R. MITCHELL me ; ZAMP 


over the range 0 < x < L, with j and & constants and g(x) an arbitrary func- 
tion of x. If the range is divided by N internal nodes, distance / apart, a simple 
finite difference replacement of (1) is 


Vn 49+ 0 Vnty + Or = Sr [7 =, 1:2, cic Nene ae (2) 


where a=1+j7h/2, b=kh?—2,c=1—jh/2, f, = h* g {(r + 1) h} and y,, 
Very» erg are the values of y at x=rh, (r+ 1)h, (r + 2) h respectively. 


Throughout the paper, the values A = + 2/7, for which (2) becomes a first order 


equation, are excluded. A solution of (2) is 


y = A+B, (3) 
1 1 us 
Yasir A amt Byntt + oo Da — Ul) fee 5 Oe ee 
7r=0 


where A, are distinct roots of the auxiliary equation ax*+ bx%+c=0, and 
A, B are constants to be determined by the boundary conditions. 


It is convenient at this stage to derive certain relationships involving A and 


4, which are used throughout the paper. We select 


_ —b+ (0?—4ac)1# _ —b— (P—4ac)1P 
A 24a HA ASS 2a : 
where A + « = —b/aand/ wu = c/a. Three cases are considered: 


(1) A, wcomplex: Put A = (c/a)? e4, uw = (c/a) e-*®, to obtain the result 


A [i= 24 (=)” sins 0 , (4) | 


where s is a variable. 
(2) A, wreal (Aw > 0): Put Aju = c?9, to get 
2(<)"sinhs6 (A>0,u> 0), 
M— w= (5a) 


c 


(—1)°2 (<)” sinhs 0 (A<0,"<0). 


(3) A, w real (Aw < 0): Put — A/u = e?®, to obtain 


c 


2 (-<)" sinhs 6 (seven) , 


A — p= (5b) _ 


c 


2 (-<)" coshs 9 (s odd). 


See ae a eee 
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Two different types of two-point boundary conditions are considered: 
(I) y=0, Y at x = 0, L. In this case, the solution given by (3) becomes 


N 


ret et awe 
bene re (A" — gw") fy_, 
Vas1= gNeT _ Nat (An+t — yn+t) 6) 
1 1 n 
gia oar aes wi 
where 0 <n < N. No solution of (6) exists in general if 
{Net —pNt1= 0, (7) 


Using (4) and (5) with s = N + 1, it is seen after some manipulation that (7) is 
satisfied only if A and yw are imaginary and 


kh? =2[1—(1- vai ne) cosa 7 | (Reecoshs 2 yah sy Non stnere) 


Mesh lengths for which the difference equation has, in general, no solution are 
said to be critical. 


(II) y= 0 at x=0, and dy/dx+ay= 6 at x=L. The derivative dy/dx +ay=B 
at x = L is replaced by the finite difference expression 


2 yy + (b—2aah) yy,, = g{(N+1)h}—2aBh. 


The solution given by (3) then becomes 


1 1 wy 
Yass = A (Mt — pret) + oa DG He far (9) 
where 
N N-1 
ty —2aBh— 7 (228) SO HD tar (0g) Oo tee 
r=0 r=0 
é (b= 2aab t= Mt") 4 (ao) (AY = 2”) 


and 0 <n <N. Solution (9) does not in general exist if 


(6 —2aah) (AN*! — pN+1) + (atc) (AN —p®) =0. (10) 


_IfA and yu are imaginary, (10) is satisfied if the mesh length h satisfies the rela- 


tionship 
Rhe == 2 [2 — (1 - fp” ee 


“Ti | (Ke, Zw Ny, (0D 
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where {h (4 — kh?) — 72312 
tang = 1 i ; . 
es ae oa Bee) 
2(1 Titht)a+ > (kh ) 


If 2 and w are real, then (10) takes the form 


b— (b?—4ac)t21N+1 2Zacah+b(1—c) — (b®—4ac)\? (12) 
area ate 2Zacah+b(1—c) + (b?—4ac)t?’ 


where b? > 4ac. In order to illustrate the existence of mesh lengths satisfying 
(12), consider the case 7 = 0. If kh? > 4, (12) can be solved to give 


2 sinhy (13) 


oS tanh2 (N+ 1) v’ 


where 


20 


Ril h 
= Fare and y= cosh! 


2 ’ 


and where for a prescribed value of N, there is a value of 6 > 1/(N + 1) cor- 
responding to any value of k? h > 2. If k < 0, (12) can be solved to give 


2 coshy 
ea (N+ 1)v’ (14)9 
where 
2a a ERT 
= he Bd 9 me SUIS eas 


and where for any value of N, there is a negative value of 6 corresponding to 


any value of (—R)# h > 0. 


3. Correspondence Between Exact Solutions of the Differential and 
Difference Equations 


In the previous section, critical values of the mesh length are given for 
which the difference equation has in general no solution. These critical lengths 
depend on the coefficients of the difference equation and the boundary conditions 
of the problem. Now the range of the problem is given by 


L=(N+1)h (15) 


and so corresponding to each critical mesh length there is a critical range of 
problem given by (15), for which the difference equation, applied at N internal 
nodes, has in general no solution. In the present section, the differential equa- 
tion is examined for critical ranges, and the latter compared with the values 
obtained in the previous section for the difference equation. 

The differential equation (1) has solution 


ye ei #2 [A elil4—h)t? x ah B Cm tea. Os 4 = p(x) ‘ (16) 


EEE 
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if 72/4 > k, and 
y = esa! [4 sin (a = - a x+B cos (h - “s ae x| apa NET) 


if 72/4 < k, where in each case (x) is a particular integral of (1) and A, B are 
constants to be determined by the boundary conditions. 

(I) If y= 0, Y at x = 0, L, the differential equation has no solution in general, 
if 72/4 < k, when 


pla 7p — K II 


1 1/2 
et AN 
(1- <r) 


Now returning to the difference equation, if (15) is used to eliminate h, (8) can 
be solved to give 


ee ppg: PEP Nae Ki { K II 8k (2k 12 
gs (1 4(N+ a be, Ne eo mab GY aiend: 7 Pe be 1)| (19) 


(niet alee Sp, (18) 


K I (4 


eee ae RI 
= C08 ar F = 1) — sin’; \ : (20) 


TJ 
If k > 77/2, a unique solution exists for L for each pair of values of K and N. 
This is obtained from (19) with the positive sign. For fixed K, as N tends to 
infinity, each value of L tends to the corresponding critical length of the differ- 
ential equation given by (18). If k < 72/2, it follows that 0< (4 R/j72?-—1) <1 
and so from (20), for prescribed K and N, the number of values of L is 0, 1, 
or 2 depending on K and N. Provided N is sufficiently large, however, there 
are two values of L for each K. The value given by (20) with the positive sign 
again tends to the corresponding critical length of the differential equation, 
and the other value tends to infinity, as N tends to infinity. 

The following table illustrates the manner in which a critical range of the 
difference equation approaches the corresponding critical range of the differen- 
tial equation as N tends to infinity. For convenience, the values; =0,k=9 
are chosen. 


Values of Critical Range 


1 
baie 3 | g 15 31 - KIT 
1 1-0204 1-0408 1-0458 1-0464 1-0472 
2 1-8856 2-0408 2-0810 2-0915 2:0944 


3 2-4636 2:9632 3-0963 B130Z 3-1416 
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(II) If y= 0 at x =0 and dy/dx +ay= fat x=L, the differential equa-— 
tion has no solution in general when 


pulaey Se ee ae (Kies 2; cee (21) 
(Fra) 
where 
Tle ee 
fecal 
tang = icin: 
— al 


if 72/4 < k, and no solution in general when 


jeanne _ 2a—j— G4 aye 
2a— f+ (72-4)? 


(22) 


if 72/4 > k. By comparison with (I), it seems probable that a critical range of the 
difference equation computed from (11) approaches the corresponding critical 
range of the differential equation given by (21) or tends to infinity as N tends 
to infinity. Consider now critical ranges of the difference equation obtained 
from (12). Suppose first 72/4 << k, then there are values of « for which (12) 
yields no critical ranges for any N. For any other value of «, it seems likely 
that there is a corresponding critical range Ly for all sufficiently large N. It 
then follows that 


hes 2 oe. 1) (i - : py”, 


and hence Ly tends to infinity as N tends to infinity. If 72/4 > k, however, the 
position is more complicated. Again there are values of « for which (12) fails 
to yield a critical range. For any other value of « it seems that for each suffi- 
ciently large N, there are two possible critical ranges, one of which tends to the 
corresponding critical range of the differential equation given by (22), and the 
other tends to infinity as N tends to infinity. 


4. Behaviour of Solutions Near Critical Points 


Three examples are now given of corresponding solutions worked out with 
ranges which are nearly critical for either the difference or the differential 
equation. In all examples, the coefficients of the difference equation are exact, 


and the value N = 3 is chosen in order to reduce the computation. The differ- 
ential equation in all examples is 


d’y 
ga tky=kx (k>0), (23) 


——— 
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which, although simple, is adequate to illustrate some of the results of the 
previous section. 


(1) In the first example the boundary conditions are y = 0, Y at x = 0, L. 
The solution of the difference equation is given by (6) with a = 1, where 
: = 1 - 2479 7, [{ 1 1/2 
A,w=1— > hii hh h(1— q hie) 
and 
i= Gr1)ee (s=0,1,..., N). 


Choose k”? h = 0-765, a value close to the critical value 2 sinz/8 for the case 
of three internal nodes, and the solution becomes 


Vy = —444-77 L + 445-02 Y, 
Vo = —629-10 L + 629-60 Y, 
Vs = —444:97 L + 445-72 Y, 


where y,, Va, Ys are the values of y at the three internal nodes. The solution 
of the differential equation is given by 


Wie. ; 
Y= spa Sink a+ x, 


with k'/? L = 3-06, and so at corresponding points, 


Y= — 82477L + 8-4977Y, 
Ye = —11-7598 L + 12-2598 Y , 
ys = — 84409L + 9-1909Y. 


(2) The boundary conditions of the previous example are retained, but this 
time a range is chosen which is close to a critical range of the differential equa- 
tion. Choose k!/? L = 3-14, which is close to the critical value z, and the solu- 
tion of the differential equation yields 


y, = —443-16 L + 443-41 Y, 
Yo = —626:85 L + 627-35 Y, 
Vz = —443:38 L + 444-13 Y. 
The solution of the difference equation with k'/* h = 0-785 gives at correspond- 


ing points 
y, = 8-7352 L— 84852 Y, 


yo = 12-2416 L —11-7416 Y, 
V3 = 85126 L — 7-7626Y. 
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In the above two examples, it should be noted that although the solutions 
disagree violently for general Y, there is exact correspondence if Y = L. The 
reason for this becomes evident when the solutions of the difference and differ- 
ential equations are rearranged for general N to give 


Y — Lae (Ant — pnt) +(n+1)h 


Vn = gN+i Net 
and 
Pein (a & Ly el? hone ee 


Vn+i Po sin k1/2 1s 


respectively. Thus in examples (1) and (2), although the range of the problem 
may be nearly critical for either the difference or the differential equation, a 
fortunate boundary value Y may still enable reasonable agreement to be ob- 
tained between the exact solutions. 

(3) In the final example, the boundary conditions are y = 0 at x = 0 and 
dy/dx +ay=B6 at x=L. The solution of the difference equation is given by 
(9) witha =1,b=k h? — 2, 


Ap=1- Zk LAR (G RAI), 
and 


jae (54.1) bes =.0 eee 


Choose k"? h = 2:5 and 6 = 1:5, values which almost satisfy (13) and so this 
example is worked out near critical conditions of the difference equation. The 
solution of the difference equation for N = 3 is 


Y= = 217593 33655 
Y,= 9250 — 10888, 
Y, = — 37125 + 43688, 
Y,= 148550 —17476 B, 


where Y,= 4 y,/L (ry =1,..., 4). The solution of the differential equation (23) i 


B-—(1+aL) 
a sink? 7 + R12 coskl2 LT 


be sinkl? % + x, 


and so the corresponding values are 


Y,= 263023 —0-1917918, 
Y, = —0-61205 + 0-3073 B, 
Y;= 5-55507 —0-3006 £, 
Y,= 2:51808 + 0-17434 B. 


i Gia) a a 


ete 
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In general, the lack of agreement between the two solutions is even moremarked 
than in the previous examples. However, the solutions correspond exactly if 
B = 8-5, and sq even although the values of kh and 6 are near critical values 
for the difference equation, a fortuitous value of f may still enable reasonable 


agreement to be obtained between the exact solutions of the difference and dif- 
ferential equations. 


5. Errors in Approximate Solutions of Finite Difference Equations 


So far, numerical solutions of the difference equation (2) have been obtained 
directly from the exact solution (3). In practice, however, an exact solution of 
the difference equation is rarely available, and approximate numerical solutions 
are obtained either by an iterative process such as relaxation or by direct me- 
thods of solution of the associated set of simultaneous equations. The accuracy 
of such numerical solutions of (2) is now investigated. 

The error equation corresponding to (2) is 


eer. he, + Ce. = Ry r= 0,1, 2).n5(N =1)I, (24) 


where €,, &43, €ry_ are the errors in y,, Yp11, Yri2 respectively and R, is the 
residual at the (ry + 1)-th node. Residuals exist for any approximate numerical 
solution of (2), whether it is obtained by an iterative process such as relaxation 
or by a direct method. A solution of (24) is 


&=A+B, (25) 


=r Rs OZnS Ny, 


r=0 


1 
eee A tt ee 


A— pb 


where A, B are constants to be determined by the boundary conditions. 
Consider first the case where the value of the function y is given at both 
ends of the range. If the boundary values require no rounding off, then 


Ey = Eng => 


and so (25) yields 


N 
n Dv = #) Finer 
a (A a L) Eni — PC ‘ai pe") Rn ares BS ore Cal pny?) (26) 
r=0 


(O0<”n=sN-1). 


If A, w are imaginary, it follows from (4) that the coefficients of the residuals 
may be large in modulus when either 6 approximately satisfies (7) or h is near 


ZAMP X/15 


C= 
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the value 2/7. In the following two examples, the errors given by (25) are evalu- 
ated for 7 = 0 and N = 3. 

(1) In the first example, & h = 0-765, and the errors at the internal nodes 
are respectively 


€, = —[445-73 Ry + 629-61 Ry + 445-02 Rs], 
& = —[629-61 R, + 890-75 Ry + 629-61 Ro} , 
&3 = —[44-502 Ry + 629-61 Ry + 445-73 Ry]. 


The errors are now examined in relation to the corresponding values of y ob- 
tained from the exact solutions of the difference and differential equations in 
the previous section, and two distinct cases arise. If Y is approximately equal 
to L, the exact solutions are in good agreement, but approximate numerical 
solutions of the difference equation are likely to involve substantial errors. On 
the other hand, if Y is not nearly equal to L, the exact solutions are in poor 
agreement, but the ratio e/y in an approximate numerical solution is likely to 
be considerably reduced. 
(2) In the second example, kh = 0-785, and the errors are respectively 


=\—[/.7-7623 R, + 11-7418 Ry 84852 Role 
€) = —[11-7418 R, + 16-2475 R, + 11-7418 R,] , 
é, = =| 84852 R, f e7418 Ry eaten 


These errors are comparatively small, and so the approximate numerical and 
exact solutions of the difference equation are likely to be in good agreement 
for any boundary value Y. 

If A, w are real, it follows from (5) that the errors given by (26) are unlikely 
to be large, provided / is not near the value 2/7. 

Consider next the case where the values of the function and its derivative 


are given, one at either end of the range. The appropriate error equation, from — 


(25), is 
a (A = Lb) En = oy (Ar = LU") (ae ue G Hes es pnt) ‘ (27) 
r=0 
where 
ala ) R. (b= 2 h) > N-1 
== it ioe 2 aa t+ ax qr Paver 
. e.. r=0 ( ve ( # ) N-1-r 


oe 


and R, (0 = 7 s N) is the residual at the (7 + 1)-th node. The magnitudes of 
the coefficients in (27) are obtained using (4) if A, us are complex and (5) if A, 
are real. The coefficients are in general considerably greater in modulus when 


a 
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A, w are real, and may differ widely in magnitude, particularly when the range 
of the problem is near a critical range of the difference equation, or h is near 
one or other of the values + 2/j. As an illustration of this unequal ‘weighting’ 
of the residuals, consider an example where N = 3, j = 0, k¥2h = 2:5,6 =155, 
and the boundary conditions are y= 0 at x =O andy’ +ay= Bat x=L. The 
values of A, 4 are —1/4 and 4 respectively and the errors, using (27), are given 
by 
& = iy = IOK + 64K, — 256 Ry: 


R= 346k, + Erne are'K L 1088 Ry 
é= 64R, — 272R, +1092R, — 4368R,, 
a= OR,  20S8 Rk, — 4368 R, + 17476 R,, 


where the coefficients are exact. It is clear that the value of R; has a much 
greater influence on the errors than the value of Ry, and consequently in any 
approximate numerical method of solution of the difference equation, the 
emphasis must be placed on reducing A, as much as possible. Unequal ‘weight- 
ing’ of the residuals in the expressions for the errors, as illustrated in this 
example, may complicate matters considerably as far as direct numerical 
methods of solution of the difference equation are concerned, since in such 
methods it is difficult to attempt to eliminate residuals.in a prescribed order. 
In iterative methods such as relaxation, however, a prescribed order of pro- 
cedure for eliminating residuals is an advantage. In this particular example, in 
which 7 = 0, the residuals which must be eliminated occur at the nodes near 
the end where the derivative is specified. This is likely to be so for general 7 
provided the mesh length is not near the value — 2/7. 


6. Concluding Remarks 


In the present paper, it is shown that there are several major difficulties 
which may confront the computor who is attempting to solve an ordinary dif- 
ferential equation with two-point boundary conditions by difference methods. 
Although the differential equation examined here is of a particular type, never- 
theless, the difficulties described are certain to exist for many other general 
types of equation. 

In the expressions obtained for the errors in the numerical solution in terms 
of the residuals, some of the coefficients are large in modulus while others are 
relatively small, when the range of the problem is near a critical range of the 
difference equation. This unequal ‘weighting’ makes it essential that the re- 
siduals with large coefficients should be reduced as much as possible, a task 
which is more easily accomplished by an iterative rather than a direct process. 
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Finally, in the solution of initial value problems by difference methods, it 
is known that the chance of serious errors arising in numerical solutions is in- 
creased considerably if the differential equation is replaced by a higher order 
difference equation (Topp [3], RUTISHAUSER [4], MircHELr and Craces [5)). 
It seems likely that the use of higher order difference replacements will also 
produce additional complications in the numerical solution of two-point boun- 
ary problems. 
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Zusammenfassung 


Es werden inhomogene Randwertaufgaben bei einer gewohnlichen linearen Dif- — 


ferentialgleichung 2.Ordnung mit konstanten Koeffizienten untersucht. Es gibt 
dann bekanntlich Ausnahmefalle, wo die Lange des Grundintervalls kritisch ist, 
das heisst die Randwertaufgabe im allgemeinen keine Lésung besitzt. Auch wenn 
die Differentialgleichung durch eine Differenzengleichung approximiert wird, hat 
letztere kritische Intervalle, und es zeigt sich, dass in der Nahe eines kritischen 
Falles die Lésungen der Differential- und der Differenzengleichung erheblich von- 
einander abweichen kénnen. Es wird ferner bewiesen, dass bei Verfeinerung der 
Maschenlange die kritischen Intervalle der Differenzengleichung gegen diejenigen 
der Differentialgleichung streben oder iiber alle Grenzen wachsen. 
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Zut Losung spezieller Warmeleitungsprobleme 


Von Rupo.tF Hormann, Berlin, Deutschland 4) 


1. Einleitung 


Die Lésung instationérer Erwarmungsprobleme unter Beriicksichtigung 
von Warmequellen mit sehr allgemeinen Anfangs- und Grenzbedingungen be- 
reitet grundsatzlich keine Schwierigkeiten. Die Durchrechnung fiihrt wegen 
der allgemeinen Auffassung des Problems jedoch auf komplizierte und un- 
handliche Ausdriicke, die sich fiir die Auswertung konkreter Beispiele kaum 
eignen. Die Behandlung spezieller Erwarmungsprobleme erscheint daher zweck- 
massig. Die vorliegende Arbeit bringt als Ergebnis die Erwarmung von Zy- 
linder, Kugel und Quader, wobei die Anfangs- und Grenzbedingungen auf Falle 
der technischen Praxis abgestimmt sind. Eingehendere Rechnungen sollen 
dem Praktiker die Anwendung der Ergebnisse erleichtern?). Zur Lésung der 
Aufgaben wird beziiglich der Zeitkoordinate die Laplace-Transformation ver- 
wendet. 


2. Das zylindrische Erwarmungsproblem 


Auf die Mantelflache eines zylindrischen Kérpers vom Radius R, der vorher 
im Temperaturgleichgewicht mit seiner Umgebung gestanden hat, wirke vom 
Zeitpunkt ¢ = 0 an eine Warmequelle, die mit der konstanten Winkelgeschwin- 
digkeit w rotiert. Ihr Vorhandensein sei die Ursache einer vorgegebenen Tem- 
peraturverteilung auf dieser Mantelflache. Der Warmefluss auf der gesamten 
Oberflache sei in jedem Augenblick proportional der dort herrschenden Uber- 
temperatur (Newtonsches Abkiihlungsgesetz). Dieser Annahme entspricht 
etwa ein auf einer Ebene rollender Zylinder (zum Beispiel Autoreifen), dem an 
der Beriihrungsflache Warme (heisse Strasse, Walkarbeit, Reibung) zugefihrt 
wird. Aber auch zylindrische Werkstiicke, die in Lagern laufen, werden durch 
diese Bedingungen erfasst. Die Temperaturverteilung im Zylinder und ihre 


~ zeitliche Anderung ist zu berechnen. 


1) OSRAM Studiengesellschaft. sib 
2) Alle auftretenden héheren Funktionen und deren Relationen finden sich ausfiihrlich behan- 


delt bei A. SomMERFELD, Vorlesungen iiber theoretische Physik, Band V1, 2. Aufl. (Leipzig 1948), 


und JAHNKE-EmpE, Tafeln héherer Funktionen, 4. Aufl. (Leipzig 1948). 


ae ee 
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Die gesuchte Temperatur O(r, ¢, 2, y) muss die Differentialgleichung der 
Warmeleitung in Zylinderkoordinaten 


020 i eo 1 00 00 00 (1) 


ore | vy Oy ' v2 Og? * Oz? oy 


befriedigen. Die Variable y enthalt die Zeit t in der Form y = at, wobei a die 
Temperaturleitzahl bedeutet. Die gesuchte Temperatur © muss ausserdem der 
vorgegebenen Anfangsbedingung 


@G@=0 tir y=0 (2) 


geniigen. Es wurde fiir die Anfangstemperatur O(r, ¢, z, 0) = 0 gesetzt, um die 
Rechnung mit unschwer zu berechnenden Zusatzgliedern nicht unndtig zu be- 
lasten. Die Grenzbedingungen sind: 


oO +hO=0 fur 2 == fy 


und (3) 
16 tnOCAR fanelyere 
gal =F (p, 279) Tire PRG 
mit 
| fir 0< p< Aly), M=LNy, w=a, | 
F(p,z,y)=\AT(y) fir 1 <yp<gly), ~=Qy+D, (4) 
| fir g<p<2a. | 


Fir O(r, y, z, y) und F(g, z, y) wird die Existenz der entsprechenden Laplace- 
Transformierten &(7, g, z, s) und f(g, z, s) und die Entwickelbarkeit in Fourier- 
sche Reihen vorausgesetzt. Mittels Laplace-Transformation gehen (1), (2) und 
(3) in die Bedingungen im Bildbereich 


079 Oe. 1 06 076 


or? a Y Or “7 ye Og? AYO se=0, (5) 
Oo 
a7 OO = 0 (2%) (6) 
und 
Oo 
ar wns aa (y= R) (7) 


liber. Wie sich leicht nachrechnen lasst, ist ein partikulares Integral 3,,,(r, , 2, S) 
von (5) durch den Ansatz 


Fin(S) a [41n(S) sin ye bin(S) COSM, Z| 1,(i /s =e or r) ein? (8) 


ne ary 


‘ 


SY 


s 


i OA Di dee et 
AVE? Cay 


RAN ees HEA ya 
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gegeben, wobei [,,(i Vs + «? 7) die Besselsche Funktion m-ter Ordnung bedeu- 
tet. Die Gréssen @,,(s), ,,(s) und «, sind durch die Grenzbedingungen (6) und 
(7) festgelegt. Das allgemeine Integral ergibt sich dann als Summe iiber die #,,,: 


le. °) +0o 
Br, p, 2,5)= Dd” S'Oulr, @, 2, 8). 


l=1 n=—0o 
Die Bedingung (6) liefert die Gleichungen 
Ayn, [&, COSH, 2 + A SiN a, 2] — Diy (a, SING, Zy — A cOSa, 2%] = 0, 
Ayn [Hy COSA; Zp + A SING, Zo] + Diy [H, SING, Zp — h cosa, %] = 0. 


Daraus folgt: 


&, SIN a, Z% — h COSa, 29 (9) 
Im “4% COSH, Z9 + ASIN A, 2p 


an = 
und 


(h? — a?) sin2 «2 + 2a,h cos2 a, 2% = 0. (10) 


Die Werte «, sind demnach die positiven Wurzeln von (10). Damit geht (8) in 


B= DS) D buals) Inés + a7) Hale) 9 (11) 
=e 
iiber mit 
H,(z) = a, COSA, (% — 2) + Asina, (% — 2) . 


a, COSA, 2% + hSina, 2 


Es gelten die Orthogonalitatsrelationen 


+39 


[ Hie) Hyle) dz = 0 era ee a 


ie 
N,= / H,(2) Hy(2) dz 
1, L (42) 
2 a, Zo (a2 + h®) + (a% — h?) sin 2 x, Z COS2 a % + 2h a, sin? 2 a, 2 
- 2 x, (a, COS) Z + A Sina, 2)? ; 
+2 Pre 
sina, 2 
K, = [ Hy(2) dz = 2, 
: 1 
— 
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die fiir die folgende Entwicklung notwendig sind. Mit 


27+ 2 


[[ He. 2s) Hie)" ae dp 


Pay DS Cin(8) H(z) e°"® ,  Cin(8) = oa 2aN, (13) 


folgt fiir den Koeffizienten b,;(s) aus (7) der Ausdruck 


R Cy(S) ; 
(n+hR) I, (ifs +o? R)—ifst of R14 (is + of R) 


Din(S) SF 


Die vollstandige Lésung des Erwérmungsvorganges lautet daher im Bildbereich 


o=RY S a (i V/s + of v) Hy wets 
fel na=co(*#@+h R) TL oer is + a? R) oy 


Bevor # in den Originalbereich riicktransformiert wird, ist es zweckmassig, fiir 
# eine Entwicklung nach Besselschen Funktionen hinsichtlich 7 und eine Ent- 


wicklung in eine Partialbruchreihe hinsichtlich der Wurzeln des a vor-_ 


zunehmen. Nullsetzen des Nenners liefert mit 


2 


iV Sink +o = Big, und s,,,—=— (ay + Bink) cast 


die Gleichung 


(n +h R) [Base k) er Bink R Trelis R) = 0. (15) : 


Die Werte f;,,;, sind demnach die positiven Wurzeln von (15). Es gelten auch 
hier Orthogonalitaétsrelationen, namlich 


R 
| EalBine?) In(Bins 7) 7dr = 0 fie + 7(,7=1,2,3,...), 
= R 
Minx = | TnlBinn 1) I,(Bink 1) 7 ar 


R 2 
= ZB, Bins R Tn(Bink R) im PP (Bink R\ mA! 2 n L(Bink k) T+1(Bine®)} , 


R 

— R 

Miox = | Bane vr) rdr = ies L(Bior R 
é lok 


(1 


ere are, “oe 


~rs 


soe) eds 


—— 


ee ee 
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Wegen (16) gilt die Entwicklung 


ie J (i Vs +o a? 7) 
(n+hR)I, Gs + oR) — ips + on} ar ae (is + of R) 


foe] 
c3% Bink Ina (Bink Fe I (B 7) 
tesatl (s + oF nn) (n + h R) i hes Ink ) » 


(17) 


ee 


und (14) geht in 


lo, ©] Cc 


~ Cin nk In in in 
B= RS SS any ences) 1 Bian?) ile) OP (18) 


t=1 n=—co k=1 


uber. Die Anwendung des Dampfungs- und Faltungssatzes der Laplace- 
Transformation liefert das gesuchte Ergebnis im Originalbereich: 


Bin qi, (8 n k) n nk” H 43 
et g dad re Re 


y (19) 
x e-ink” | cl) (r) emi dr , 
‘i 
wobei in (19) wegen (4) und (13) fiir 
hk, —inD =e Ny, 
nly) = + ee 1 iy) extn 


zu setzen ist. In der technischen Praxis interessiert meistens nur der stationare 
Endzustand des Erwarmungsvorganges, wobei in vielen Fallen die zeitliche 
Unverdanderlichkeit von T(y) angenommen werden kann. Mit T(y) = Ty folgt 
dann aus (19) unter Beriicksichtigung von (17) 


¢RhTy cr iS 
AS eran ay 


I, (i Vo? — in Qr) ——_ (20) 


ie : —e 
(n+h R) I, (iva a2 —in® R) —ifo®—inQ RIyy, (i i Vo? —inQ R) 


Bie — 1) 


in(y—2y) 
ime H(z) e ; 


Der Wellencharakter von @,,,, ist leicht zu erkennen. 


3. Das Erwairmungsproblem bei der Kugel 


Auf der Oberflache einer Kugel vom Radius R, die vorher im Temperatur- 
gleichgewicht mit der Umgebung gestanden hat, sei vom Zeitpunkt ¢ = 0 an 
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eine Temperaturverteilung 
| 0 fir 0<e<al), mOo)=2y, 
=/AT(O fir (<9 oO, (= eyTe, | (21) 


G(d, @, ¥) 
0 fir g,— 9 < 22 


vorgegeben. G entstehe auf dieselbe Weise wie beim zylindrischen Problem. 
Ferner gelte auch hier das Newtonsche Abkiihlungsgesetz. Diese Annahmen 
erlauben zum Beispiel die angenaherte Berechnung des Erwarmungsvorganges 
bei einem rotierenden Satelliten im Strahlungsfeld der Sonne. Aber auch Er- 
warmungsprobleme bei Kugellagern sowie die Bestimmung der Temperatur 
im Erdinnern kénnen durch diese Annahmen behandelt werden. Die Berech- 
nung der gefragten Temperaturverteilung geschieht auf dieselbe Weise wie im 
2. Abschnitt. 

Die gesuchte Temperatur O(r, 3, y, y) muss die Differentialgleichung der 
Warmeleitung in Kugelkoordinaten 


1 0%(r @) 1 0 [. 9 00 1 20 00 

ice Riabr eV Perle Gc eree Ss cena: igh By ee 
befriedigen und die Anfangs- und Grenzbedingungen 

@=0 fir y=0, 92 420=6 fir r=R (23) 


erfiillen. Auch hier wird die Existenz der Laplace-Transformierten 


L {O(r, 3, p, v)} = I(r, 9, @, 5) 
und 


L{G(d, v, y)} = g(9, 9, 5) 


sowie ihre Entwickelbarkeit in Fouriersche Reihen vorausgesetzt. Die Bedin- 
gungen (22) und (23) lauten im Bildbereich 


a en, 
y ore = 2 sind | 08 


Oo 1 Oe 


al ng resin? ® dg? sd=0 (24) 


(sino 
und 


5 thi=g fir += R> (25) 


Eine lineare Kombination partikularer Integrale von (24) liefert den Lésungs- — 
ansatz: 


i = a> oy SIRS) p,(i /sr) P"(cosd) e'™? , 


n=0 m=—oo 


Vol. X, 1959 Zur Lésung spezieller Warmeleitungsprobleme 239 


Die F,"(cos#) sind die sogenannten zugeordneten Kugelfunktionen m-ter Ord- 
nung, sie erfiillen die Differentialgleichung 


1 d 


| apm 
sind do (sino © = 


+ {n(n ug? ees | Per =0. 


‘sin? ® f-” 
Es gelten die Relationen 


+1 


[Pi(x) P(x) dx = 0 fir Da 0 (lym = 12,3, «.)5 | 
icy +1 i \! | 
at Y Mf m ae 2 +m 26 
t ASS =| fos (x) oa (x) dx = 2 nti : (n ai m)! , | ( ) 
ay 
—m (1) @—™)! pm po j : 
B," = (-1) oS Pe P=, Pr=0 fir man. | 


Die y,(¢ 7) erfiillen die Differentialgleichung 


1 Ory, n(n+1 
ep + {ot |p, = 0. 


Y or 
Ferner gilt 
ae Ynl) = T,41)2(9) (halbzahlige Bessel-Funktion) , 
_ sing __ sing — g cose eg Sena an 
Yar ap gs ies 02 > Yn= (-@) laa) Yo» ep 
an 


Sag | gat 
do mE Wnt 


- Die Integrationskonstante d,,,,(s) ist durch (25) festgelegt. Aus (25) folgt, wenn 


ES re ee AY wh Ee AN 


fiir g die Entwicklung 


) 


g= s S gam(6) P™(cosd) e'"? 
n=0 m=—0o 
j mit ; ‘ 128) 
(8, p, s) e'™® Pm(cos 9) sin b db dp 
Eam(s) = eG 
- angesetzt wird, : 
dnm(s) = R EnmlS) 


(n+h R) Wn ( is R) —iVs R Yam is R) 
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Die Lésung des Erwarmungsvorganges lautet daher im Bildbereich 


CO 


248 ~ &nm(s ) 22 
eae ee (nth R) y, (i ys R) = 4/5 B vaya (i ays R) (29) 
x Yn (i ys 1) Pm ere 


Fiir die Riicktransformation wird hier & nach den W» hinsichtlich v und in eine 
Partialbruchreihe hinsichtlich der Wurzeln des Nenners entwickelt. Nullsetzen 
des Nenners liefert mit 


die Gleichung 

(n +h R) Vn(Oin R) — Fin RWn41(Gin k) =0, (30) 
woraus die o;, als die positiven Wurzeln berechnet werden kénnen. Es gelten 
die folgenden Orthogonalitatsrelationen : 


R 


Hits Vr(Oin 7) VnlOjn A) dy= 0. fra ] (I, 1 = Sees se) P 


0 


R 
Kin = | Pvalint) Palin?) dr = 
é so (31) 
x { Cin R [vi(Orn R) is Use One k)] hy (2 n+ 1) Vn(Oin R) Yn+1(Fin R) } ? 
R 
— ” R2 
N,= [r Wo(O107) ar = rae W1(07, R) . 
% O 
Daraus folgt 
“s R Yn ( iV/sr) 
(n +h R) Yn (i Vs R) —iys R Vnit 0 Vs R) 
(32) | 


RE Fin Yn+1 (im FR) c 
b= CG + ofn) ("+h R) K,, PnlOin 7) ; 


und fiir 9: 


T= R237 $7 Snmls) Opn Yas (Orn B) ve 
a 4 ~. (s + of,) (n+ hR)K,, Pm) Pr oh 


: 
; 


a 


Ss ee PN CS NU EAE, Ow aa 
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Die Riicktransformation von (33) liefert das gesuchte Ergebnis im Original- 
bereich: 


QO = R3 $3: ae s Sin Vn+1 (Cin ¥) | 


l=1 n=0 m=—co (n +h R) Kin 


y : (34) 
x Palin) Pa (C088) em 0219 [Gyy(e) etm" de, | 
0 
wobei in (34) wegen (21) fiir 
G _¢ ,-imQy ‘ hile ee ty F m 4 
nm) = ~nm © . mit a m cp ay may mM a T(9) yo (cos?) sind dd 


0 


zu setzen ist. Fiir den stationaren Endzustand @,,,,; ergibt sich aus (34) mit 
Beriicksichtigung von (32) 


Osa = ¥ =R D> Sra nm so (cos @) e 


ae Sim Qr) eWim? (35) 
liegt past . Vn (i V—imQr)e —imQy 


(n+hR) yp, (iV—imQ R)—iV~—imOR y,,, (iV~—im@R) 


_ ——_— 


Der Wellencharakter von @,,,; ist auch hier leicht zu erkennen. 


4. Ein spezielles Erwa4rmungsproblem der Elektrotechnik 


Es werde eine Metallplatte betrachtet, die vorher im Temperaturgleich- 
gewicht mit der Umgebung gestanden hat und vom Zeitpunkt ¢= 0 an von 
dem konstanten Gleichstrom J durchflossen wird. Im Innern der Platte wird 
dann pro Volumen- und Zeiteinheit die Warmemenge W, und an den Stirn- 
flachen (| x | = x,) durch den als gleichmassig verteilt angenommenen Uber- 
gangswiderstand pro Flachen- und Zeiteinheit die Warmemenge Q, erzeugt. 
Auf den restlichen Teilen der Oberflache gelte das Newtonsche Abkithlungs- 
gesetz. Fiir die Anfangstemperatur sei auch hier O(x, y, z,0) = 0 gesetzt. Die 
Temperaturverteilung in der Metallplatte und ihre zeitliche Anderung ist zu 
berechnen. 

Die gesuchte Temperatur O(x, y, z, 7) muss die Differentialgleichung der 
gees mit inneren Warmequellen 


(36) 
02 02 02 yi 
| A es ae t 3 +55, t=a t, ==, A= Warmeleitfahigkeit) 


242 Rupotr Hofmann ZAMP 
erfiillen mit der Anfangsbedingung 
O=0 fir ne (37) 


Als Grenzbedingungen ergeben sich 


00 ss 
Be i f6 fir x= +%9 (4 - &o) (38) 


(da auf den Stirnflachen die Warmestrémung durch die erzeugte Warmemenge 
QO, gegeben ist) und 


we EhO=0 fir y=,  hO=0 fir z=1%, (39) 


Mit L{O(x, y, z, t)} = Hx, y,z,s) und den bekannten Regeln der Laplace- 
Transformation gehen (36), (37), (38) und (39) in die Bildbereich-Bedingungen 


Ads ieee (40) 
Oo qd fz 
aes Gs fiir) 4 = =b %5 (41) 
Oo 
72 +hO=0 fir y=+y, se +hO=0 fir z=+2 (42) 


iiber. Eine lineare Kombination partikularer Integrale der homogenen und 
inhomogenen Differentialgleichung (40) liefert den Lésungsansatz: 


be As (43) 
ae) % [c, ,(S) sin? s + e,x + d,,(s) cost iene, «| K,,(y) Px(2) 


n=1 k=1 


K,(y) = Bn COSBy (Yo — Y) + hsinB, (Yo — ¥) 
RK Bn COS By Vo + h sin By Vo 


und 


P,(z) _ VrCOS7p (49 — 2) + Asiny, (% — 2) 
Yu COSY, 29 + ASiny;, 2 ; 


Um (42) zu geniigen, miissen £,, und y, als positive Wurzeln aus den Gleichungen ~ 
(h? — B2) sin2 By Vo + 2h B, cos2 By Vo = 0, 


445 
(A? — yh) sin2 y, 2 + 2h y, cos2 yp 2 a 


| 


| 
| 


~Y Yee aS 
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bestimmt werden. Es gelten die Orthogonalitatsrelationen 


+ Vo 


[Kul Kjly)dy=0 fir n+ 4-7 =1, 2,3,...), 

Ho 

+¥o 

x aA 
oe =| K,(y) dy =— ene , 

n 

+ Yo 

Ry = | Kyly) Kaly) dy 
zo 

= 2 Bn Yo [Bh mi h®) + [Bn Jx h*| sin 2 Bn Vo COS 2 Bn Vo + 2 h By sin? 2 me Vo 
A - Bn (Bn cos By Vo = he sin By, Yo)? ‘ 


++ Bq 


| Pxl@) Pyle) dz = 0 fir k+j (k,j7=1,2,3,...), 


2siny, 
V; =| P,(z) dz = 25D Ye *0 
Vie 
ie 
+ 


2Vx Ye COSY; 29 = hsiny;, el? 


(40) liefert, wenn fiir 


und fiir 
Ene = Bn + Vi 
gesetzt wird, 
Wo 
s (s a Eni) ; 


bn «(S) — Ane 
Schliesslich folgt aus (41) 
és 0 und dny = A fo, 


a sVs+ e2, sini /s + 2, % 


tL 


Die vollstandige Lésung lautet daher im Bildbereich: 


cost Vs + one x 


[oe] oc) | Wo ome i a = 
a xe py ARS. Be K,,(y) P.(2) | s (s + €2,) 7 sVs+2 sini //s + 2p % 


n=1 k=1 


(46) 
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Zur Riicktransformation in den Originalbereich ist eine Umformung von (46) 
zweckmassig. Auf Grund der Entwicklung 


cosi /s + 2, ¥ 2 a 1)? cos x (+/*9) | 
i Vs+ ef, Sint V/s + en % 0 j=0 (S + OFyx) (47) 
2 if mt \2 2 2 
loins a fS a Bn al A , 
geht (46) in 
n=1 k=1 
Wo ee (1) cosl x (4/%9) | 
‘e s (s sia En) Xo i=0 Ss (s Ale OlnR) j 
iiber. Wegen 
el Ee ieee ee ee 
ei Ge a es 
und bei Beachtung von (47) bekommt man das gesuchte Ergebnis im Original- 
bereich: 
O= Dats K,(y) P,(2) 
m=1 h=1 fe | (49) 
7 { qo cosh €p 5 a | Wo. —en kT 2 Go il cos/ z (~/%o) vag 
Soiscesinlic, give | 62% (bese ) tig fod (1) a. a ; 


Fiir die Temperatur @,,,; im stationaren Endzustand (t > oo) folgt aus (49) 


Of er ee = Waar oe cosh é,, x K ly) Plz een) 

An speziellen Fallen der technischen Praxis wurde gezeigt, dass auch bei 
schwierigen Erwarmungsproblemen der Durchrechnung mit klassischen Me- 
thoden keine wesentlichen Hindernisse im Wege stehen. Die Lésung instatio- 
nadrer Erwarmungsprobleme unter Beriicksichtigung von Warmequellen und 
sehr allgemeinen Anfangs- und Grenzbedingungen kann immer additiv aus einer 
Anzahl von Teillésungen aufgebaut werden, die alle grundsatzlich auf die 
gleiche Weise gefunden werden, wie die in dieser Arbeit abgeleiteten Ergebnisse. 


Summary 


The heating of cylinder, sphere and cube is computed, the starting and boundary 
conditions being specified for cases of technical importance. The solution of non- 
stationary heating problems (considering heat sources, and starting and boundary 
conditions of very general nature) can always be synthetized by superposition of 
a finite number of partial solutions, which can, in principle, be found in the same 
manner as those worked out in this paper. 


(Eingegangen: 21. Oktober 1958.) 
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Die Ausgleichsrechnung nach der Methode der kleinsten Quadrate gehort 
heute zu den klassischen Stiicken der numerischen Mathematik und ist langst 


jahrigen Gauss kaum mehr richtig zu wiirdigen wissen. 


1) Institut fiir angewandte Mathematik der ETH. 


ZAMP X/16 


_ so selbstverstiandlich geworden, dass wir deren Entdeckung durch den siebzehn- 
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Wenn man allerdings die Praktiken verfolgt, nach denen vorher zum Bei- 
spiel im Vermessungswesen Widerspriiche enthaltende Beobachtungen in miih- 
samer und undurchsichtiger Weise notdiirftig zum Stimmen gebracht wurden, 
und wenn man sieht, welche zihen Widerstande bei der Einfiihrung der Gauss- 
schen Methode zu itiberwinden waren, dann wird doch klar, dass damit ein be- 
trachtlicher Schritt getan wurde. ; 

Die Theorie, 1794 von Gauss ersonnen, beinahe gleichzeitig von LEGENDRE 
ebenfalls gefunden und von diesem als erstem 1806 publiziert, feierte ihren er- 
sten Triumph bei der Berechnung der Bahn des neuentdeckten Planeten Ceres 
durch Gauss im Jahre 1801. In der folgenden, mehr als 150jahrigen Entwick- 
lung wurden Theorie und Art der praktischen Durchfiihrung auf eimen hohen 
Stand gebracht. 

Die Methode der kleinsten Quadrate kann durch wahrscheinlichkeitstheo- 
retische Uberlegungen sauber begriindet werden, und sie ist von den tiberhaupt 
sinnvollen Méglichkeiten der Ausgleichung die mathematisch zuganglichste. 
Die gesuchten Gréssen treten als Unbekannte eines linearen Gleichungssystems, 
der sogenannten Normalgleichungen, auf. Zu deren Auflésung gibt es viele 
Methoden; wohl die bekannteste ist diejenige der sukzessiven Elimination der 
Unbekannten nach dem GauBschen Algorithmus bzw. einer moderneren Va- 
riante desselben Verfahrens. 

Allein schon von Gauss wurde ein ganz anderer Weg vorgeschlagen und 
auch praktisch erprobt: die iterative Auflosung von Gleichungssystemen durch 
Naherungsverfahren, welche in jedem Zwischenschritt eine bessere Approxi- 
mation der Unbekannten fiefern. Da nun auch die exakten Lésungen der Nor- 
malgleichungen doch nur einen Kompromiss zur Ausgleichung der widerspruchs- 
vollen Messungen darstellen kénnen, liegt gewiss die Frage nahe, ob nicht schon 
die angenaherte Lésung des Gleichungssystems annehmbare Resultate liefern” 
kénnte. 

Die Entwicklung der iterativen Lésungsmethoden wurde besonders in den 
letzten fiinfzig Jahren auf breiter Basis betrieben; unter dem Stichwort «Re- 
laxationsrechnung» hat sich bereits eine betrachtliche Literatur angesammelt. 
Mit der Anwendung dieser Methoden auf die Ausgleichsrechnung erhebt sich 
weiter der Wunsch nach Verfahren, welche die Aufstellung der Normal- 
gleichungen umgehen {1]?). Ein Ziel unserer Arbeit ist die Diskussion der An- 
wendung gewisser Relaxationsmethoden auf die Ausgleichsrechnung (2. Teil). 

Mit dem fritheren Abbrechen eines Naherungsverfahrens stellt sich von 
selbst das Problem, gewisse Abschatzungen fiir die erreichte Genauigkeit zu 
erhalten, und zwar méchte man im Falle der Ausgleichsrechnung insbesondere 
Schranken fiir die Quadratsumme der Verbesserungen, welche zu einem Mini- 
mum gemacht werden soll. Bekanntlich unterscheidet man zwischen vermit- 


) Die Ziffern in eckigen Klammern verweisen auf das Literaturverzeichnis, Seite 280. 
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telnder und bedingter Ausgleichung; die beiden Arten verhalten sich bei einer 
geometrischen Interpretation dual zueinander. Nun liefert im Falle det ver- 
mittelnden Ausgleichung jede Naherungslésung fiir die Unbekannten eine tri- 
viale obere Schranke fiir die erwahnte Quadratsumme, wahrend bei der be- 
dingten Ausgleichung im allgemeinen iiberhaupt keine Schranken bei vorzei- 
tigem Abbrechen entstehen und nur bei speziellen Iterationsverfahren wahrend 
des Verlaufs der Rechnung untere Schranken angegeben werden kénnen. Um 
beidseitige Schranken zu gewinnen, miisste man die beiden Ausgleichungsarten 
kombinieren. Es wird sich zeigen, dass dies unter Umstanden nicht schon beim 
Beginn der Rechnung, sondern erst vor dem Abbrechen getan werden muss. 

Das zweite Ziel der vorliegenden Arbeit besteht in der Ubertragung einer 
von SYNGE [2] und anderen besonders im Zusammenhang mit partiellen Dif- 
ferentialgleichungen der Physik entwickelten Idee auf die Ausgleichsrechnung 
zwecks Gewinnung solcher Schranken (3. Teil). Dies geschieht zundchst in all- 
gemeiner Form, das heisst ohne Bezug auf eine bestimmte Iterationsmethode, 
sodann aber speziell im Hinblick auf das in der zitierten Publikation [1] be- 
sprochene Verfahren. 

Es ist mir eine angenehme Pflicht, an dieser Stelle Herrn Professor STIEFEL 
dafiir zu danken, dass er zu dieser Arbeit den Anstoss gab und mir auch in der 
Folge immer wieder mit Anregungen beistand. Ebenfalls zu Dank verpflichtet 
bin ich Herrn Professor Kosoip, der mir im Gesprach den Standpunkt des 
Geodaten etwas naherbrachte. 


1. Formulierung der Ausgleichsrechnung in Matrizenschreibweise 
1.1 Allgemeine Erklérungen 


In der Ausgleichsrechnung geht man von einem Satz von n Werten aus, 
welche den Charakter von Messresultaten haben. Wenn nun diese Anzahl n 
grosser ist, als zur eindeutigen Bestimmung des fraglichen Systems notwendig 
ware, so werden sich infolge der Messungenauigkeit gewisse Widerspriiche ein- 
stellen. Zu deren Beseitigung hat man an den gemessenen Gréssen Verbesse- 
rungen derart anzubringen, dass alle Nebenbedingungen erfiillt werden. Gleich- 
zeitig soll man méglichst «in der Nahe» der Messungen bleiben. Diese letztere 
Forderung kann auf verschiedene Arten interpretiert werden ; in der bekannten 
GauBschen Methode der kleinsten Quadrate sucht man ihr dadurch gerecht zu 
werden, dass man die Quadratsumme der Verbesserungen zu einem Minimum 
macht. 

Die erwahnten Nebenbedingungen treten bei der vermittelnden Ausgleichung 
dadurch in Erscheinung, dass sich die » Messungen als lineare (bzw. lineari- 
sierte) Funktionen der eigentlichen Unbekannten von geringerer Anzahl y aus- 
driicken lassen miissen («vermittelnde Beobachtungen»). Bei der bedingten 
Ausgleichung hingegen sind die ausgeglichenen Grdssen selber die Unbekannten, 
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welche einen Satz von f linearen (bzw. linearisierten) Bedingungsgleichungen 
zu erfiillen haben (8 < »). Die gemischte Variante der vermittelnden Ausglei- 
chung mit Bedingungsgleichungen zwischen den Unbekannten werden wir in 
dieser Arbeit nicht beriicksichtigen. Ferner wird durchwegs mit Messungen 
vom Gewicht 1 gerechnet, um eine Uberlastung des Textes zu vermeiden. Der 
Fall ungleicher Gewichte lasst sich leicht auf diesen Spezialfall zurtickfiithren. 

Wir fiihren unsere Betrachtungen im u-dimensionalen euklidischen Raum 
R,, durch und verstehen unter einem Punkt oder einem Vektor einen Satz von 
n reellen Zahlen; die Verwendung der beiden Begriffe geschieht durchaus im 
iiblichen, geometrisch anschaulichen Sinne. Die «Ebene» definieren wir dann 
etwa als Punktmenge, die mit zwei verschiedenen Punkten auch deren Ver- 
bindungsgerade enthalt. 

Fiir die algebraische Beschreibung werden wir von der Matrizenschreib- 
weise Gebrauch machen, welche sehr bequem und der Komponentendarstellung 
entschieden vorzuziehen ist (vgl. [3]; dieses Buch enthalt ein Kapitel tiber 
Ausgleichsrechnung). 

Vektoren und Punkte sollen im folgenden mit kleinen lateinischen Buch- 
staben bezeichnet werden. Durch das Skalarprodukt 


(%, Y) = %1 V1 Xe Votes +40 Vn 


wird eine Metrik eingefiihrt. Der Absolutbetrag eines Vektors (bzw. Abstand — 


eines Punktes vom Nullpunkt) ware dann 


|x| =), x) - 


Die lineare Transformation schreiben wir in der Form der Multiplikation — 


mit einer Matrix, das heisst y = A x als Abkiirzung fiir 
Ve Oy Me Ogg Ng a een a Lee 


wobei im allgemeinen yw + v, das heisst die Matrix A nicht von quadratischer 
Gestalt ist. Die transponierte (durch Vertauschung der Zeilen mit den Spalten 
entstandene) Matrix bezeichnen wir mit 47, Wenn man einen Vektor als Matrix 
mit nur einer Spalte auffasst, kann man das Skalarprodukt auch in der Form 


(x, y ) =x" uy 
schreiben. Hingegen ergibt die Multiplikation 
xy? = D 


eine Matrix D mit den Elementen d,, = x, y; als sogenanntes dyadisches Pro- 
dukt der beiden Vektoren x und y, k § yi es Pro 


ese es <ccnantnna~s 
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Wir werden im allgemeinen Matrizen mit grossen lateinischen, skalare 
Gréssen mit kleinen griechischen Buchstaben bezeichnen. Indizes sollen im 
folgenden, soweit nichts Besonderes vermerkt ist, verschiedene V ektoren, und 
nicht deren Komponenten unterscheiden. 


1.2 Vermittelnde Ausgleichung 


Die Variablen x (Vektor mit » Komponenten), welche gemessen werden, 
seien lineare Funktionen der Unbekannten y (y Komponenten): 


Gime Gl ef (1) 


Dabei sind die y in einfacheren Fallen die ersten y Komponenten des Vektors x, 
das heisst so viele, als frei verfiigbar sind. Das bedeutet, dass der obere Teil von 
C die y-reihige Einheitsmatrix ist. 

Liegt fiir x ee Messung / vor, so sind die Unbekannten y so zu bestimmen, 
dass der Betrag von v = x — / minimal wird. Man nennt 


v=Cy+f-l 


die Fehlergleichungen (” Gleichungen, y Unbekannte, wobei y < 1; sonst lage 
gar kein Ausgleichungsproblem vor). 

_ Bei der vermittelnden Ausgleichung kénnen wir / immer als Null annehmen 
bzw. durch eine Koordinatenverschiebung zu Null machen und somit die 
Fehlergleichungen in der iiblichen Form schreiben: 


v=Cy—l. (2) 


Die gesuchten Werte erhalt man durch Auflésen der GauBschen Normal- 
gleichungen 
aes Saag 3 ie (3) 


(y Gleichungen, y Unbekannte). Sehr oft interessiert auch das ausgeglichene x, 
das wir mit s bezeichnen. v hat die Bedeutung einer Verbesserung, welche zur 
Messung / zu addieren ist. 

Geometrische Interpretation im R,,: / soll «méglichst gut» als Linearkombi- 
nation der y Spaltenvektoren von C dargestellt werden, das heisst, der Abstand 
\v| des Lésungspunktes s von / soll minimal sein. Das bedeutet, dass s die Pro- 


3 jektion von / auf F’ ist, wobei F’ die von den Spaltenvektoren c; von C auf- 


gespannte Ebene ist, oder v=s—TZJ | F’. Daraus folgen unmittelbar die 
Normalgleichungen 
Gim= 0) oder eben C7 Cy=C" 4. 
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In Figur 1 ist der Fall n =3, y=2 dargestellt, wo F’ eine gewohnliche 
zweidimensionale Ebene wird. 


Figur 1 


Beispiel zur vermittelnden Ausgleichung. 


1.3 Bedingte Ausgleichung 


Fiir den Vektor der Unbekannten x liege die Messung / vor. x ist so zu be- 
stimmen, dass (v,v) = Min mit v = x — 1, wobei x die Bedingungsgleichungen 


Bx+b=0 (4) | 


(6 Gleichungen, » Variable, 8 < n) zu erfiillen hat. Es liegt somit eine Extremal- 
aufgabe mit Nebenbedingungen vor, welche man in bekannter Weise mit Hilfe 
der Lagrangeschen Multiplikatoren lést. Wenn wir diese mit — 2#,; («= 1, 
2,...,B) bezeichnen und die ¢; zum sogenannten Korrelatenvektor ¢ zusammen- 
fassen, ergeben sich ausser (4) als weitere Bedingungen fiir die Unbekannten ¢ | 
und x die Korrelatengleichungen 


%= BELT. (5) 
Durch Einsetzen von (5) in (4) erhalt man die Normalgleichungen 
BBY t= —(BY8) (6) 


(8 Gleichungen, 6 Unbekannte). Deren Lésung #, in (5) eingesetzt, gibt endlich © 
die eigentlichen Unbekannten ~x. 

Geometrische Interpretation im R,,: Die Bedingungsgleichungen (4) defi- 
nieren eine (7 — B)-dimensionale Ebene F’, und zwar jetzt als Schnittgebilde von — 
Hyperebenen, welche durch ihre Normalvektoren gegeben sind. / soll «méglichst — 
gut» durch einen Punkt aus F’ approximiert werden, das heisst, gesucht ist der — 
Punkt in /” mit dem kleinsten Abstand |v| von J. Das ist natiirlich wieder der 
Fusspunkt des Lotes von J auf F’; er wird aber diesmal bestimmt als Schnitt- — 


| 
| 
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punkt von F’ mit einer Ebene F”, deren Parameterdarstellung 
KS Bt + 
= 1+ {Linearkombination der Zeilenvektoren von B] 


lautet. F” ist total senkrecht zu F’, was man einsieht, wenn man (4) in der 
Form B(x — f)=0 schreibt, wobei f ein beliebiger Punkt aus F’ ist, also 
Bf+6=0. (Wir nennen zwei Ebenen im R,, «total senkrecht», wenn jede 
Gerade in der einen senkrecht auf jeder Geraden der anderen Ebene steht und 
die Summe ihrer Dimensionen gleich » ist.) 

Figur 2 zeigt das Beispiel n = 3, 8 = 2; F’ ist die durch die Ebenen E, und 
E, bestimmte Gerade. F” wird durch die in / abgetragenen Zeilenvektoren }, 
und 6, von B aufgespannt. 


Figur 2 
Beispiel zur bedingten Ausgleichung. 


Diese Formulierung der bedingten Ausgleichung ist der in der Geodasie eher 
gebrauchlichen mathematisch vollkommen aquivalent. Man betrachtet dort die 


 Verbesserungen v als Unbekannte und schreibt die Bedingungsgleichungen in 
_ der Form Bv + w = 0, wobei w = B/1 + b die sogenannten Widerspriiche sind. 
_ Korrelaten- und Normalgleichungen lauten dann: 


v= BTt und BBTt=-w 


_ Bei unserer Form werden offenbar die Beziehungen des nachsten Abschnittes 
: etwas durchsichtiger. 


Man beachte, dass die Unbekannten bei 1.2 Koeffizienten einer Linear- 


_kombination von n-dimensionalen Vektoren sind, bei 1.3 hingegen Koordinaten 


oY ins RR lie ea ee 


yw 
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1.4 Beziehung zwischen den beiden Arten von Ausgleichung 


Da die beiden betrachteten Falle dem gleichen geometrischen Problem 
aquivalent sind, den Punkt s in einer Ebene F’ zu finden, der den kleinsten 
Abstand von einem bekannten Punkte / ausserhalb F’ hat, liegt der Versuch 
nahe, das eine Problem auch algebraisch auf das andere zuriickzufihren. 

Liegt zum Beispiel eine Aufgabe der bedingten Ausgleichung vor, so handelt 
es sich darum, y = » — B Vektoren zu suchen, welche zu den Zeilenvektoren von 
B orthogonal sind, und mit ihnen als Spaltenvektoren die Matrix C zu bilden. 
Ist f ein Punkt aus F’, so lautet die Transformation 


w= Cytf. (7) 


Durch Einsetzen von (7) in (4) sieht man, dass die Bedingungsgleichungen in der 
Tat fiir beliebige y identisch erfiillt sind (B C=0, Bf+6b=0). 

Praktisch hat dieser Ubergang allerdings nur dann einen Sinn, wenn man 
(1 — B) der Unbekannten als neue Unbekannte y nimmt und dann fiir die Kon- 
struktion von C die Bedingungsgleichungen auf einfache Art nach den tibrigen 
f Unbekannten auflésen kann; oder es sei, dass man durch zusatzliche Uber- 
legungen das Problem von vornherein auf beide Arten formulieren kann. 

Bei einer iterativen Behandlung der Normalgleichungen bedeutet nun: 

Vernuttelnde Ausgleichung: sukzessive Verbesserung der Unbekannten y. 
Der Punkt x bewegt sich in F’, bis er auf den Schnittpunkt mit F” fallt. 

Bedingte Ausgleichung: sukzessive Verbesserung der Unbekannten ¢. Der 
Punkt x bewegt sich in F'", bis er auf den Schnittpunkt mit FP’ fallt. 


1.5 Die mittleren Fehler der ausgeglichenen Unbekannten 


In diesem Abschnitt sollen einige Formeln iiber mittlere Fehler zusammen- 
gestellt werden. 


Der mittlere Fehler der Messung / ist bekanntlich bei vermittelnder Aus- 
gleichung 


-/S, 


m = 742 : 


(Handelt es sich um ein und dasselbe Problem, so stimmen die beiden Aus- 
driicke wegen n — y = f tatsdchlich iiberein.) 


Fir den mittleren Fehler M einer Funktion q (J) gilt das Fehlerfortpflan- 
zungsgesetz (1; = Komponente von 1): 
is 


Mat SF 


bei bedingter Ausgleichung 


ee ee ee ~ 


Pees 


ee 


Vol. X, 1959 Iterative Lésung und Vehlerabschatzung in der Ausgleichsrechnung 253 


Bei den ausgeglichenen Unbekannten handelt es sich nun um eine lineare Vektor- 
funktion der Messungen /: 


a(l) = Al 4+- const 


mit konstanter Matrix A. Eine einfache Rechnung ergibt fiir den mittleren 
Fehler M; der 7-ten Komponente von a: 


M; = 4;m’, 
wo q;; das 1-te Diagonalelement der Gewichtsmatrix 
QG=AA™ 


ist. Die ausserdiagonalen Elemente von Q, die gemischten Gewichtskoeffi- 
zienten, braucht man ebenfalls, sobald fiir Funktionen der Unbekannten mitt- 
lere Fehler angegeben werden sollen. 

Wir wollen nun fiir die beiden Falle der Ausgleichsrechnung die Gewichts- 
matrizen aufstellen. Das jeweilige A findet man durch Auflésen der Normal- 
gleichungen nach den Unbekannten. Um anzudeuten, ob es sich bei den aus- 
geglichenen Variabeln um die Unbekannten der Normalgleichungen oder um 
Koordinaten im x-Raum handelt, werden wir die Bezeichnungen Q und P ver- 
wenden. 

Vermittelnde Ausgleichung: 


Chtiywcry, 
= (CO CY Pers, 
A 
Oieided anion. (8) 


Fiir die mittleren Fehler der ausgeglichenen Messungen s brauchen wir 


aCe CC OMe 
A 


Pes asece Orc. (9) 


Bedingte Ausgleichung: 
BBA t= —=(Bl-+d), 
t= —(B B’)~' (BI +5), 
= B? ¢+d =[E-— Bt(B B*)-' B}i + const , 


—~— 


P=AA?=E~— BT(B B)-! B (E = Einheitsmatrix). —_—(10) 
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2. Numerische Behandlung der Ausgleichungsaufgabe 
2.1 Ubersicht tiber verschiedene Relaxationsverfahren 


Nachdem im 1.Teil die theoretischen Grundlagen der Ausgleichsrechnung 
zusammengestellt worden sind, soll nun die numerische Auflésung der Normal- 
gleichungen behandelt werden. Wir beschranken uns dabei auf die Diskussion 
einiger iterativer Methoden, welche unter den Begriff der Relaxationsrechnung 
fallen. Diese Methoden erlauben ein Abbrechen der Rechnung vor dem Er- 
reichen der theoretisch exakten Lésung (welche man infolge der Rundungs- 
fehler im allgemeinen ohnehin nicht erhalt). In diesem Falle kann aber der 
Rechenaufwand erheblich vermindert werden, wenn man die explizite Auf- 
stellung des Normalgleichungssystems vermeidet, was durchaus méglich ist. 

Im weiteren sollen folgende Bezeichnungen verwendet werden: Es liege das 
symmetrische Gleichungssystem 

Ay =k 


fiir die Unbekannten z vor mit der Losung 4 = A-1 k. Dann ist r = k — A zder — 


zum Versuchsvektor z gehérige Residuenvektor, und wir nennen 


f(z) = [h—2, A (h — 2) 
die Fehlerfunktion. 
Die primitivste Methode zur Auflésung von linearen Gleichungssystemen, 
welche in unseren Kreis passt, das Einzelschrittverfahren [4], bei welchem in 
jedem Relaxationsschritt nur eine Unbekannte verbessert wird, konvergiert 


ziemlich langsam, so dass sich die im Falle der Ausgleichsrechnung nicht so — 


grossen Vereinfachungen gegentiber den besseren Verfahren kaum lohnen. 
Die Verfahren bester Taktik in der Terminologie von STIEFEL [5], das 
heisst solche Verfahren, bei denen ohne Riicksicht auf die Vorgeschichte der 


Fehler in jedem Iterationsschritt méglichst klein gemacht wird, werden durch — 


Gradientenmethoden, besonders die bekannte Methode des steilsten Abstieges, 


geliefert [4]. Die Korrektur der Unbekannten wird hier proportional zum Resi- — 


duenvektor angesetzt : 
Az=)r, 


da man leicht zeigen kann, dass 7 bis auf das Vorzeichen gleich dem Gradienten 
einer gewissen quadratischen Funktion F(z) ist, welche im Lésungspunkt ihr 
Minimum hat. Der Faktor 4 wird in jedem Schritt so berechnet, dass F (als 


Funktion von A), und damit iibrigens auch die oben definierte Fehlerfunktion — 
/(z), minimal wird. Auch bei diesem Verfahren ist die Konvergenz oft sehr — 


schlecht. 


ae 
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Weniger Rechenaufwand erfordert das etwas grébere Gesamtschrittverfahren 
[4]. Bei diesem wird der Parameter A im Gegensatz zur vorhergehenden Methode 
tiber den ganzen Auflésungsprozess weg konstant gehalten. Fiir die Wahl dieses 
Wertes braucht man aber Schranken fiir die Eigenwerte von A, und diese 
stehen hier im allgemeinen nicht zur Verfiigung, da die Matrix A gar nicht auf- 
gestellt werden soll. 

Wenn man nach dem Gradientenverfahren bester Strategie fragt, das heisst 
nach dem Verfahren, das nach einer vorgegebenen Anzahl » von Schritten den 
Fehler am kleinsten macht, so kommt man, wie in [5] bewiesen wird, unter 
Zugrundelegung von f(z) = [kh — z, A (hk — z)] als Fehlermass auf die Methode 
der konjugierten Gradienten (cg-Verfahren, siehe [4], dort als n-Schritt-Verfahren 
bezeichnet ; erschépfend diskutiert in [6]). Die Methode des steilsten Abstiegs 
und das Gesamtschrittverfahren sind vereinfachte Varianten davon; der Mehr- 
aufwand beim cg-Verfahren ist jedoch im Falle der Ausgleichsrechnung relativ 
gering, da dort die beiden hauptsachlich ins Gewicht fallenden Multiplikationen 
Matrix x Vektor pro Schritt auch bei den einfacheren Verfahren erforderlich 
sind. 

Wir werden daher im folgenden die Anwendung des cg-Verfahrens in der 
Ausgleichsrechnung etwas genauer diskutieren. Vorher soll aber noch gezeigt 
werden, dass die erwahnten Strategietitberlegungen auch bei Fehlermassen, die 
sich bei der Ausgleichung natiirlicherweise aufdrangen, ebenfalls zum cg-Ver- 
fahren fihren. 

_ Im Falle der vermittelnden Ausgleichung liegt es zum Beispiel nahe, nach 
demjenigen Verfahren zu fragen, bei welchem (v,v) nicht nur monoton, sondern 
mdglichst rasch im Sinne der besten Strategie abnimmt. Die Normalgleichungen 
sind 

C’Gy=CT, 


so dass nach den obigen Bezeichnungen A = C? C, k= Cl, z= y zu setzen ist. 


Dann gilt: 
a (v, v) = (Cy —1, Cy —]) 


= (y,C7 Cy) —2(y, C7) + 0,9) 
= (z,A2z)—2(2,Ah)+ (2) 
=[h— 2,4 (h—2)]—(h, Ah) + (9) 
= f(z) + const. 
Damit ist gezeigt, dass (v,v) bis auf eine Konstante gleich der Fehlerfunktion des 


Normalgleichungssystems ist und somit das cg-Verfahren die Methode bester 


Strategie darstellt. 
Dieselbe Uberlegung hat bei der bedingten Ausgleichung mit Korrelaten 
keinen Sinn, da wir uns dann wahrend der Relaxation nicht in der Ebene F’ der 
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zulassigen Punkte bewegen, sondern in der Orthogonalebene F ”, und zwar im 
allgemeinen ausgehend von « = /, v = 0, so dass (v, v) gelegentlich eee muss 
(beim cg sogar monoton). Die Korrelaten- und Normalgleichungen lauten: 


x=BTt+1l, BBTt=—(Bil+)), 


wobei diesmal A = B B’, k= —(B1 +b), z =¢ ist. Der Lésungspunkt ware 
s=BhA+l. te ey 

Nimmt man nun das Quadrat des Fehlervektors s — x (im x-Raum, nicht im 
Raume der Unbekannten der Normalgleichungen) als Fehlermass, so ergibt 
sich aus 


(s — x, s — x) = [BT (h—2), BT (h—2)] = [h—-#, B BT (h —- 2) 
= [h — 2, A (h— 2)] = (2) 


wieder das cg-Verfahren. Zu beachten ist, dass dieses Fehlermass, im Gegensatz 
zu (v,v), wahrend der Rechnung nicht herauskommen kann, da ja die Lésung 
s nicht bekannt ist. 

Zum selben Resultat kame man auch bei der vermittelnden Ausgleichung 
mit s—x=C (h—y). 


2.2 Die Methode der konjugierten Gradienten 


Es sei wieder das Gleichungssystem 


Az 


mit der n-reihigen symmetrischen und positiv definiten Koeffizientenmatrix A 
gegeben. (Diese Bedingungen sind fiir unsere Normalgleichungen sicher erfiillt, 
sofern die Spalten von C bzw. die Zeilen von B linear unabhangig sind, was man 
iiberhaupt voraussetzen muss, falls A nicht singular werden soll. Dieser Fall 
witirde bedeuten, dass man mehr Unbekannte, als das Problem Freiheitsgrade 
hat, bzw. tiberfliissige Bedingungsgleichungen eingefiihrt hatte.) 

Wir gehen aus von einem beliebigen Anfangsvektor 2). Der 7-te allgemeine 
Iterationsschritt besteht darin, dass man in Richtung eines Gewichtswelstots b; 
fortschreitet und damit die neue Naherung 


24 =, +4; 2; 


gewinnt, wobei p; aus dem i-ten Residuenvektor und dem Gewichtsvektor des 


vorhergehenden Schrittes linear kombiniert wird. Dieser zweite Summand fallt 
im ersten Schritt weg. 


Vol. X, 1959 Iterative Lésung und Fehlerabschitzung in der Ausgleichsrechnung 27 


Die genaue Rechenvorschrift lautet (vgl. [61) ; 


7. Schrnitt: 
1 = Pyo=hk—-Aad, 
pe (Yo, %o) 
(Po, A Po) ’ 


Pp =F; Sin Piss 


i (te %) 


Bn APS 
241 = % +4; 8; . 


Die Gewichtsvektoren sind paarweise konjugiert (Name!), das heisst, es gilt: 
(pj, Ap) =O0 fir 147. 


Die Residuenvektoren bilden ein Orthogonalsystem (Beweise siehe [6]). Aus der 
letzteren Tatsache folgt aber, dass nach spatestens ” Schritten theoretisch die 
exakte Lésung erreicht wird. 

Gleichzeitig mit den Unbekannten kann die inverse Matrix A-! erhalten 
werden, wenn in jedem Schritt (¢ = 0, 1, 2, ..., m — 1) folgende Zusatzrechnung 
mitgefiihrt wird: 


tiage Pibi Start- Ant 0 12 
ji. ae of ri (Pp, A Pj) , start: ea oe ( ) 


(p p’ ist das eingangs definierte dyadische Vektorprodukt). 

Nach genau m Schritten sollte theoretisch die exakte Lésung 47*= A 
erreicht werden. 

Der Betrag des Fehlervektors h— xz sowie die Fehlerfunktion f(z) = 
[h — z, A (h — z)| nehmen bei jedem Schritt monoton ab. Dies gilt aber nicht 
fiir den Betrag des Residuenvektors 7. Mit (r,7) als Fehlermass fiir die Strategie- 
betrachtung kommt man auf eine Variante des cg-Verfahrens, welche hier 
jedoch ausscheidet, da sie bei der Ausgleichsrechnung pro Schritt eine Matrix- 
multiplikation mehr erfordert. 


—1 
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2.3 Vermittelnde Ausgleichung nach der Methode der konjugierten Gradienten 


Das eben beschriebene Rechenverfahren soll nun auf die Normalgleichungen 
der vermittelnden Ausgleichung tibertragen werden (siehe [1]; der Algorithmus 
ist dort komponentenweise geschrieben). 

Durch Substitution von 


Ave CO, RS ON z= 955 Se 


in der Rechenvorschrift (11) erhalt man direkt die Formeln (13), wobei nur 
folgende zwei Punkte zu beachten sind: 

1. Wie bisher in der vermittelnden Ausgleichung setzen wir C y — / = v, und 
speziell beim Start mit einer Versuchslésung yy: C yo —/ =v 9. Im ersten 
Schritt wird damit 


k—Az,=CTl— CT Cy, = =C y=7,;- 


In allen Schritten gilt ferner C Ay; =Av;; oder, wenn A,;e;=Ay,, dann 
Ap g; = Av; = Ax,. 

Es hangt vom speziellen Problem ab, ob man die eigentlichen Unbekannten 
y oder die ausgeglichenen Messungen x oder beides verlangt. Fiir die Abwick- 
lung des Iterationsprozesses geniigt es, v; mitzufiihren ; die x, und y; k6nnen auch — 
erst nachtraglich berechnet werden. : 

2. Die Matrix des Normalgleichungssystems, A = C’ C, tritt im allgemeinen 
Schritt zweimal auf, und zwar einerseits bei der Berechnung von /; im Nenner: 


(p;, A pi) = (é;, CaS é;) = (Ce, Ce) = (9, q;) =0 


$9 


anderseits bei der Berechnung des Residuenvektors. Aus der Gleichung — 
¥;=7;4 — A, A Pi ergabe sich 7, =7,_, — CT Av,_,. Ebensogut kann aber§ 
die urspriingliche Definitionsgleichung 7; = k — A z, verwendet werden. Diese 
liefert>7,;='—C? v, 

Es zeigt sich, dass man an den beiden erwahnten Stellen so umformen kann, . 
dass in der Ausgleichsrechnung die Matrix des Gleichungssystems dank ihrer 


speziellen Struktur gar nicht mehr explizit auftritt und sich daher ihre Vvor-— 
gangige Berechnung eriibrigt. 


ng 


ee 
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So erhalten wir die folgenden Formeln: 
7. Schritt: 


mete Pe 
Vo = ey = —C* 5 Co = (Yo: Yo) » 


A= 20 


Ayy = Age, Vy = Vo t+ Ayo, 
Ax = Aq qo » X= X%qy + Axa, 


Mis — Ag O,, 0 = 05°F Ab," 


t-ter allgemeiner Schritt (i = 1, 2, ...): 


’ 


AEs. 
Cex, 


Oe Ft By Sys 


r= —CT 0, B= ry te oh 


GZ=C-e;, 0; = (GG) » 

eee Gt 

a he 
Ay, = 4; ¢;, Viner = Vi + AY; , 
Ax, = 4; 9; » X41 = %; + AX, , 
At Ag: , t=, Ad. . 


Man koénnte, um die Stabilitat des Prozesses zu erhdhen, die v; (und damit 
indirekt die 7;) nicht rekursiv berechnen: 


U4, = 0; + Ay, , 


sondern in jedem Schritt neu. Die letzten zwei Zeilen von (13) werden dann 
ersetzt durch: 
pe CH; 44 ets 


Die y,; miissen jetzt mitgefiihrt werden. 

Diese Variante erfordert eine zusdtzliche Matrixmultiplikation pro Schritt, 
und es fragt sich, ob der Mehraufwand durch die Verbesserung der Resultate 
aufgewogen wird. 
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2.4 Bedingte Ausgleichung nach der Methode der konjugierten Gradventen 


Genau gleich wie im vorherigen Abschnitt wird die Rechenvorschrift fiir die 
bedingte Ausgleichung hergeleitet. Diesmal lauten die Substitutionen, mit 
welchen in (11) zu gehen ist: 


ABB hea -(Bl4d), =i eee 


Da fiir die Korrelaten ¢ im allgemeinen keine Ausgangsnaherung vorliegen wird, 
starten wir den Algorithmus mit f, = 0 und damit x) =/, v9 = 0. 
So wird im ersten Schritt das Ausgangsresiduum 
t= khk= (Blew): 
Wie oben wird das Skalarprodukt im Nenner von /,; umgeformt: 
(p, A p,) = (¢;, B Be) = (BT @,, B® ¢) = (9, g,) = 0; - 


Die Residuenberechnung geschieht im 7-ten allgemeinen Schritt wie folgt: 


r,=k—-Az,=—(B1+))—B Bt =—([B(B i+) 45] 


=—(Bx,+3). 


Die ¢;, welche zwar die Unbekannten des Normalgleichungssystems sind, aber — 
nur die Rolle von Hilfsgréssen spielen, werden kaum interessieren, ebenso sind 
die v; fiir die Abwicklung des Iterationsprozesses nicht notwendig. 

So erhalt man folgenden Algorithmus: 


1. Schritt: 
% == —(Bl+5), De = (Fortale 
qo = BY ey, 59 = (Yo %o) » 
Ee Cy 
ro * eri) (14) 


Axy = Io Io X= X%y + Ax, 
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i-ter allgemeiner Schritt (i = 1, 2, ...): 


oe: 

Sia Or.’ 

See T Soy Oem 

q; = B? ¢., OG; = (Yj Qi) » — 

Aa 

t 0; 

(At, = Ay, tar = + At) 
At. = hg. , i414 = %j + Ax;, 
(Av; = A; q; , U5441 ae ¥; VF Av,) 


Die Neuberechnung von v; (und damit indirekt der Residuen) in jedem 
Schritt wiirde bei der bedingten Ausgleichung darin bestehen, dass die letzten 
beiden Zeilen von (14) ersetzt wiirden durch: 


=. Ae se 
Vip = Bb, Mar =a, tl. 


Auch hier hatte man eine Matrixmultiplikation mehr, und die ¢; miissten bei der 
Rechnung mitgefiihrt werden. 


2.5 Berechnung der Gewichtskoeffizienten 


Fiir die Berechnung der in Abschnitt 1.5 mit Q und P bezeichneten Ge- 
wichtsmatrizen legen wir Gleichung (12) zugrunde, nach welcher zum Aufbau 
der Inversen von A in jedem cg-Schritt der Anteil 


Di Pi. 
Ae ee 
a (Pi, A Pi) 
zu bestimmen ist. So erhalt man unter Beachtung der in (13) und (14) ein- 
gefiihrten Bezeichnungen folgende Ausdriicke fiir die AQ und AP: 


Vermittelnde Ausgletchung: 
Q=(C7C)*=A™, 


AQ;=AA;", 
e;eF 
AQ; = os : Start: Og= 0. (15) 


a 


ZAMP X/17 
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P2C(C? OC)" Cex An ew 
AP, = AIGA # C= CAA, Cas 


T 
A= Gide Start: P= 0n (16) 


0; 


Bedingte Ausgleichung: 
P2B= BBB) B= ESB Aa oe 
AP, = A{E — B? A™ Bl; = = Bi AAguiae 


1 
AP, = #4, Start: P, = E. (17) 


Da man bei der Berechnung von A~! keine gute Konvergenz erwarten darf, 
ist die Bedeutung von so erhaltenen Naherungswerten bei einem fritheren Ab- 
brechen der Iteration sehr zweifelhaft. Man beachte allerdings, dass die Nahe- 
rungen fiir die Diagonalelemente von P bei der bedingten Ausgleichung, von 1 
ausgehend, monoton abnehmen (die o,; sind Vektorquadrate) und dass somit 
wenigstens obere Schranken geliefert werden, die besonders interessieren. Wird 
iiberdies beim selben Problem gleichzeitig auch vermittelnd ausgeglichen, so 


nahert man sich denselben Werten von unten; denn die beiden P, welche man — 


nach vollstandiger Durchfiihrung der beiden Iterationsprozesse erhalten wiirde, 
stimmen natiirlich theoretisch exakt tiberein. 

Die obere Schranke 1 fallt beim Q der vermittelnden Ausgleichung dahin, 
was verstandlich ist, da diese Werte, im Gegensatz zu den Koeffizienten von P, 
mit einer physikalischen Dimension behaftet sind. 


3. Abschatzungen fiir die Quadratsumme der Verbesserungen 
3.1 Allgemeine Bemerkungen 


Wenn zur Auflésung der Normalgleichungen ein Naherungsverfahren an- 
gewendet wird, so ist es von Interesse, wahrend der Rechnung Abschatzungen 
fiir die Quadratsumme (v,v) der Verbesserungen zu gewinnen. Denn damit hat 
man einen gewissen Anhaltspunkt fiir die Genauigkeit der erreichten Naherung. 
Das (v,v) geht tiberdies in die Ausdriicke fiir die mittleren Fehler ein, fiir welche 
besonders eine obere Schranke erwiinscht ist. 

Fiir unsere Abschatzungen brauchen wir gleichzeitig Naherungspunkte aus 
F’ wie auch aus F’", Wir miissen somit entweder in einem der beiden Raume 


ausgleichen und vom letzten Naherungspunkt aus einen méglichst nahen Punkt | 
im anderen Raume aufsuchen, oder aber die beiden dualen Arten von Ausglei- — 


chung von Anfang an parallel nebeneinanderfiihren. 
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Zur praktischen Durchfithrung des erwahnten Uberganges ware im einzelnen 
zu bemerken: Bei der vermittelnden Ausgleichtng bewegen wir uns standig im 
Raume fF’ und haben somit die Aufgabe, von einem Punkte x’ einen Punkt x” 
in F" aufzusuchen, das heisst einen Punkt, dessen v-Vektor orthogonal ist zu 
den Spaltenvektoren der Matrix C des Fehlergleichungssystems, also C v’ 
= C (x" — 1) = 0 («duale Bedingungsgleichungen »). Bei einer Aufgabe mit 1 Mes- 
sungen und y Unbekannten darf somit iiber (n — y) Komponenten von x” frei 
verfiigt werden; das Naheliegendste ist natiirlich, die entsprechenden Kompo- 
nenten von x’ zu tibernehmen. Fiir die restlichen y Komponenten liefert die 
obige Bedingung ein Gleichungssystem. Nun wird man allerdings im allgemeinen 
kaum, nur um Abschatzungen fiir (v,v) zu erhalten, nochmals ein Gleichungs- 
system vom gleichen Umfang wie dasjenige der Normalgleichungen auflésen 
wollen. Jedoch ist in manchen praktischen Fallen die Matrix C sehr speziell ge- 
baut und hat viele verschwindende Koeffizienten. Wenn zum Beispiel die Un- 
bekannten y einfach die y ersten Messvariabeln sind, dann ist der obere Teil von 
C gerade die Einheitsmatrix, und die Bestimmung der fehlenden Komponenten 
von x” wird trivial. Im 4. Teil wird ein einfaches Beispiel dieser Art ausgefiihrt. 

Bei der bedingten Ausgleichung ist der umgekehrte Ubergang ganz entspre- 
chend auszufiihren: vom letzten Punkt «” aus F” ist ein Punkt x’ aufzusuchen, 
der in F’ liegt und somit die Bedingungsgleichungen erfillt: Bx’ + b= 0. 
Auch hier gilt das oben Gesagte, das heisst, die Berechnung von x’ lohnt sich 
nur, wenn dies dank der speziellen Struktur von B ohne grosse Miihe geht. Im 
4. Teil ist ein Beispiel aus der Vermessungskunde gegeben, bei welchem die Be- 
stimmung der restlichen 6 Komponenten (f = Anzahl der Bedingungsglei- 
chungen) eben noch ohne Auflésung von Gleichungssystemen méglich ist. 

Hierbei ist zu beachten, dass man beim verfriihten Abbrechen einer be- 
dingten Ausgleichung in einem Punkt aufhért, der die Bedingungsgleichungen 
nicht erfiillt. Falls man dies nicht in Kauf nehmen will, ware auch aus diesem 
Grund der Ubergang zu einem x’ erwiinscht. 

Die im folgenden beschriebenen Methoden beruhen auf der Ubertragung 
eines Gedankenkreises, der im Zusammenhang mit der Behandlung von par- 
tiellen Differentialgleichungen der Physik entstanden ist, auf das Gebiet der 
linearen Algebra. Es handelt sich um die von SyNGE und anderen entwickelte 
«Method of the Hypercircle». Wir werden die Theorie so wenden, dass sie fiir 
unsere Zwecke brauchbar wird. (Fiir eine knappe Einfiihrung siehe [7], aus- 
fiihrlich in [2]}.) 


3.2 Grobe Abschatzung 


Es mégen die beiden folgenden dualen Probleme vorliegen: 


I: Gesucht der Punkt in F’ mit dem kleinsten Abstand vom festen Punkt / 
aus fF". 
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Il: Gesucht der Punkt in F” mit dem kleinsten Abstand vom festen Punkt 7 
ausil. 
F’ und F"” sind dabei zwei total senkrechte Ebenen im oben definierten 
Sinne. Die gemeinsame Lésung von I und II ist dann der einzige Schnittpunkt s 
von F’ und F" (Figur 3). 


Fe 


Figur 3 
Schematische Darstellung der beiden total senkrechten Ebenen F’ und F”. 


Sind nun x’ und x" zwei beliebige Punkte aus F’ und F'", so ist («’ — s, x” — s) 
= 0, das heisst, s liegt nach dem Satze von THALES auf einer Hyperkugel tiber 
x’ —x" als Durchmesser. Daraus kénnen zwei Nutzanwendungen gezogen 
werden: 

1. Nimmt man den Mittelpunkt c = («’ + x")/2 der Hyperkugel als Nahe- 
rungslosung fiir unser Problem, dann weiss man, dass der Betrag des Fehlers 
genau gleich dem Radius @ = |x’ — x"|/2 ist. 

2. Betrachtet man die Hyperkugel iiber f — /, so liefert nach dem Satz von 
PYTHAGORAS 


P= S75) XE 80 = Sy eae | 


jede untere Schranke fiir den zu minimalisierenden Abstand in Problem I eine 
obere Schranke fiir Problem II, und umgekehrt. 

Wir wollen diesen zweiten Punkt etwas naher ausfiihren und zugleich ver- 
suchen, eine brauchbare Formel zu finden. Vorerst soll x’ die Rolle von f 
spielen, x’ und x” seien zwei Naherungslésungen fiir s. x" gibt sofort eine obere 
Schranke fiir IT: 

(¥ = 'S, 4° — "sy a a ee 


Ferner ist 
(i — s,l—s) (g's = s\a ( pae 


und Subtraktion der beiden Zeilen ergibt 


@—s,l—-s)\2(0— 2) low) = (= 2", #2"). (18) 
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Damit ist (/ — s, / — s) von unten eingeschrankt (Problem I!). Es ist vielleicht 
nicht tiberfliissig zu bemerken, dass die rechte Seite dieser Ungleichung im all- 
gemeinen nicht durch den naheliegenden Ausdruck (J — x", 1— x") ersetzt 
werden darf. Die obere Schranke liefert x’: 
(i—s,l—s)S (l—x',l—x'). (19) 
Denken wir nun wieder an die Ausgleichsrechnung, so erhalten wir fiir das 
(v,v) des ausgeglichenen Punktes, indem wir (18) und (19) zusammenfassen, 


(v’, v') — (x’ — x", — x") Ss (v, v) S (v’, v’'), (20) 


mit den Bezeichnungen s —/ =v, x’ —l=v’. 

Entsprechende Ungleichungen kénnten fiir das aquivalente Problem II an- 
gegeben werden; diese interessieren jedoch hier nicht, da es sich bei beiden Aus- 
gleichungsaufgaben lediglich um verschiedene Formulierungen von Problem I 
handelt. 


3.3 Feinere Abschatzung 


Mit entsprechendem Mehraufwand lassen sich noch bessere Schranken an- 
geben: 

Es sei ein Satz von »’ linear unabhangigen Vektoren z;, parallel zu F’, 
gegeben. Die z/ spannen, von / aus abgetragen, eine v’-dimensionale Teilebene 
G’ von F’ auf (v’ < Dimension von F’). Genau dasselbe mit den zweifach 
gestrichenen Bezeichnungen gibt eine Ebene G” durch /. Es gilt somit 

(z;, 2;) =0 (21) 
fiir alle 7 und 7. G’ und G” werden im allgemeinen keinen Punkt gemeinsam 
haben. 

Wir suchen nun die beiden Punkte w’ in G’ und w” in G" mit dem kleinsten 
Abstand | u’ — u” 

Mit dem Ansatz 


watt J) % 2%, (22) 
Ww =d+ a ae ee (23) 
lautet die Minimumbedingung: 
ae ee wu all Bertin!) = 0, GH U2) c4 >’), (24) 
it (u’ —u",u' —u") = —2 (u' —u",2z;)=90 @ =I, 2, eee), 


te FARA MA, 
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Durch Einsetzen von (22) und (23) in (24) und (25) erhalt man 
2% ( zi, z,) = (L—f,-2,) (dere yo, wee als 


oa | es) = Sah ii 2.) (fe 12 cera . 


ein Gleichungssystem fiir die « und «/’. Setzen wir noch voraus, dass die z; und 
z! je ein Orthogonalsystem bilden: 
/ UA 
(z;, z;) = 0 | 
(et, 7) =0 | 


so wird die Auflésung trivial, und wir erhalten: 


fury 7%; 


(7 = fea) 
w=f+y4 7 4 (26) 
ee 2 ; se zi. (27) 


Bemerkung: Fiir einen beliebigen Punkt x’ aus F’ gilt: 


OF | 
ferner  (u” — wu’, z/) =0 G=125088 
0 


also (u" — x’, 2) = 


Das heisst aber, dass «” die Projektion von F’ auf G” ist. Ebenso ist w’ die Pro- — 


jektion von F” auf G’. 
Wie schon festgestellt, liegt s auf der Hyperkugel 


(f—*,1—2)=0 (28) 


uber f — / als Durchmesser. 


Der Punkt s liegt aber auch auf sémtlichen Hyperebenen (Hyperebene — 


= (n — 1)-dimensionale Ebene) : 


(u’ —x,2f)=0 (6 =1,2,...,9%), (29) 


(u" — %, 2) =O P= 12s ae (30) 


Deren Schnittgebilde ist eine (m — »’ — y")-dimensionale Ebene H, welche 
sowohl auf G’ als auch auf G” senkrecht steht. 


ee ee eee 


ee ee eee 
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Also liegt s auf dem Hyperkreis, den die Ebene H aus der Hyperkugel (28) 
herausschneidet. (Wir halten uns hier an die Syngesche Terminologie, wonach 
die Hyperkugel der geometrische Ort der Punkte mit konstantem Abstand von 
ae festen Zentrum, der Hyperkreis das Schnittgebilde Hyperkugel-Ebene 
ist. 

Es lasst sich leicht einsehen, dass fiir einen beliebigen Punkt x des Hyper- 
kreises (u’ — x, wu" — x) =0 gilt; das heisst, w’ — w” ist ein Durchmesser: 
Wegen der Orthogonalitaét von G’, G’ und H 


folgt aus (u’ —f, wu” —1l)=0 
und (x —wu', u" —1l)=0 


Zusammen mit (x«—f,/—x) =0 
ergibt sich ie — f mie = = . 
Mit (f—w', u” —x)=0 
schliesslich (x — uw’, u" — x)= 


Ahnlich wie friiher ziehen wir zwei Folgerungen : 

1. Nimmt man den Mittelpunkt c = (u’ + u")/2 des Hyperkreises als Nahe- 
rungslésung fiir unser Problem, so ist der Betrag des Fehlers genau gleich dem 
Radius 0 = | uw’ — u" |/2. 

2. Fiir (s — 1, s — 1) — oder, wenn wir wieder die Bezeichnungen der Aus- 
gleichsrechnung einfiihren, fiir (v, v) — gilt folgende Abschatzung: 


(u” —Lu" —l)<v,v)S (w’—Lu—l). (31) 


Die Richtigkeit der ersten Ungleichung folgt daraus, dass s — w" in H liegt und 
damit auf w” — I senkrecht steht; die zweite ist trivial. Die dualen Beziehungen 
gelten natiirlich fiir (s — f, s — /). 

Die einfachste Anwendung von (31) ergibt sich, wenn wie im Falle der 
«groben Abschatzung» in jedem der Raume F’ und F” nur je ein Naherungs- 


; punkt x’ bzw. x” als bekannt vorausgesetzt wird. Die Ebene G’ degeneriert 


TAK APSA AXP Oa 


ee TN 


dann zum Punkt x’ = uw’, G” ist die Gerade durch x” und /. Mit v’ = x’ —], 
vy” = z"= x" — 1 ergibt sich aus (27) sofort: 


a L A soln" ae 
und damit 
(v’, we es ee ‘ F 32 
(an v") = (v, v) — (v ’ v’) a ( ) 
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Wir sind insofern berechtigt, (31) und speziell (32) eine «feinere Abschat- 
zung» zu nennen, als wir hier unter Voraussetzung der durch die Vektor- 
systeme z; und z/ aufgespannten Ebenen G’ und G" die beiden optimalen 
Punkte x’ und x” —ndmlich eben w’ und u” — angeben, mit denen eine Abschat- 
zung nach (20) gemacht werden kann: w’ liefert als Projektion von / auf G’ tat- 
sichlich die kleinste obere Schranke von allen Punkten aus G’. Gehen wir 
anderseits von einem beliebigen Punkt x’ in G’ aus, so ist wu" dessen Projektion 
auf G” und somit (x’ — u’, x’ — wv’) der kleinste Betrag, den wir von (x’ —1, x’ — 1) 
abziehen kénnen, das heisst, uw” gibt die grésste untere Schranke von allen 
Punkten aus G”. Wegen der Orthogonalitat von x’ — wu" auf G" ist aber 
"a1, x = 1) — («' =u", x — we) gerade gleichn” — Ta == 7); 

Man kann ferner zeigen, dass bei Hinzunahme eines neuen 2; bzw. z; die 
obere bzw. untere Schranke tatsdchlich verbessert wird: Sei w” der alte, w/ der 
neue Punkt. Da die Glieder der Summe (27) voneinander unabhangig sind und 
zudem orthogonal vorausgesetzt wurden, ist «!. — uw” proportional zu dem zu- 
sitzlichen z/ und steht senkrecht auf uw” — /, also 


(x 


(u" —1,u" —l) > (w"—-Lu"—-l). 


Entsprechend gilt 
(u", —f,w, Se. (Hi Af pays 
Da aber 
(u’ — 1, u' —1) + (w' —f, wv’ — f) = (1—f,l — f) = const 


ist, gilt fiir die obere Schranke 
(uw, —1,u', —1l) << (vw —i,u'—]). 


Damit sind die letzten zwei Behauptungen erwiesen. 


3.4 Zusammenfassung der Abschdtzungen 


Es liege eine Ausgleichungsaufgabe vor (vermittelnd oder bedingt), welche, 
in die Sprache der Geometrie des R,, tibersetzt, folgendem Problem aquivalent 
ist: Gesucht ist der Punkt s in der Ebene F’ mit dem kleinsten Abstand von 
einem festen Punkt / ausserhalb F’. Mit F” werde die zu F’ total senkrechte 
Ebene durch / bezeichnet. Man méchte nun, ohne den exakten Lésungspunkt s 
zu kennen, Schranken fiir das Quadrat seines Abstandes von J, (v, v) haben (in 
der Ausgleichsrechnung: Quadratsumme der Verbesserungen). 

a) Grobe Abschatzung: Aus F’ und F" liege je ein Naherungspunkt x’ und 
x" vor; x’ ergabe sich als Naherungslésung bei einem iterativen Verfahren der 
vermittelnden Ausgleichung, x” bei einem solchen der bedingten Ausgleichung. 
Bezeichnet man die Abweichungen von / mit v = 5 — I, v' = x’ — 1, so gilt fiir 
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das (v, v) des Lésungspunktes die Abschatzung 
(v', 0) — (x — &" eh = x") S (0, v) S (v', v’). (20) 


b) Feinere Abschatzung: Liegt aus F’ und F” je ein Vektorsystem vor, wo- 
durch die beiden Teilebenen G’ und G” bestimmt werden, dann erhalt man die 
engsten Schranken (20), indem man fiir x’ und x” die beiden Fusspunkte w’ 
und w” der kiirzesten Verbindung von G’ und G" setzt. (20) wird in diesem 
Falle zu 


(w" Lu" VS 0) Sw 10D] (31) 


und speziell, falls wie unter a) nur je ein Punkt x’ und x” vorliegt, als Verbes- 
serung der unteren Schranke von (20) 


yf ay tVo 
(z » U re < (v, v) < (v’, v’) 


(v", v my 


(32) 


mit 
v= xI— 1, vv =x" =. 
Die Berechnung der Punkte w’ und #”, welche im allgemeinen die Auflésung 


von linearen Gleichungssystemen erfordert, reduziert sich bei orthogonalen 
Vektorsystemen auf die Auswertung der einfachen Formeln (26) und (27). 


3.5 Anwendung auf das Verfahren der konjugierten Gradienten 


In diesem Abschnitt kommen wir nochmals auf das Verfahren der konju- 
gierten Gradienten zuriick und untersuchen, welche Folgerungen sich im An- 
schluss an die bisherigen Ausfiihrungen fiir dieses spezielle Verfahren ergeben. 

Es wurde bereits darauf hingewiesen, dass die Gewichtsvektoren des cg- 
Verfahrens ein konjugiertes System bilden: (f,, A p;) = 0 fiir 7 + 7. In der 
Ausgleichsrechnung fiihrt dies auf: 


(g,,C" €¢;) 

= (4: 9;) = 0. (33) 
(eq. B B* @;) 

Diese Beziehung sowie die Orthogonalitat der Residuen 


(r;,7;)=0 fir 1 #7 (34) 


- kénnen als Rechenkontrollen verwendet werden, indem man zum Beispiel in 
- jedem Schritt (g;, g;-1) und (7;, 7;,) bildet. 
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Wegen Ax, = A, q; bedeutet aber (33), dass auch die Korrekturvektoren im 
x-Raume ein Orthogonalsystem bilden: 


(Ax;,Ax,)=0 fir 747. (35) 9 


Die Tatsache, dass die Gewichtsvektoren beim allgemeinen cg-Verfahren 
ein konjugiertes System bilden und dass daraus bei der Anwendung auf die 
Ausgleichsrechnung eine Orthogonalitét im héherdimensionalen Raume der 
Messvariabeln wird, zeitigt immerhin einige Konsequenzen : 

Zuniachst ergibt sich daraus unmittelbar der folgende 


Satz: Bei der bedingten Ausgleichung nach der Methode der konjugierten — 


Gradienten wachst (v, v) monoton. 


Dabei ist natiirlich vorausgesetzt, dass mit %) = / gestartet wird. 


Es seien im folgenden x’, x” beliebige Punkte aus F’, F”; x;, x/ hingegen © 


Naherungspunkte in der vermittelnden bzw. bedingten Ausgleichung nach 
cg. Dann gilt wegen (35): 
(x; —s,x; -—l=0, 


aber auch 
(s — x’, x/ —1) =0 
und damit 
(x; — x’, x/ —l) =0 (36) 
und speziell 
(x; —f, x7 —]l) =0. (37) @ 


(f ist, wie erinnerlich, der feste Punkt aus F’, der zum ersten Male bei der Ein- 
fiihrung des linearen Zusammenhanges zwischen Unbekannten und Messvaria- 


beln bei der vermittelnden Ausgleichung auftrat und der bei der bedingten — 


Ausgleichung Anlass zum konstanten Vektor b in den Bedingungsgleichungen 
gibt. Es wird im folgenden vorausgesetzt, dass die vermittelnde Ausgleichung 
in f, das heisst mit yy = 0 beginnt.) 

Entsprechend gilt natiirlich 


(}—x", df) =0 
und speziell 


(x; -—lx,-—f)=0. (38) 


(37) und (38) reduzieren sich auf 


(x — f, 0) =0 (39) 


und konnen wiederum in jedem Schritt als Rechenkontrolle verwendet werden. — 
Die geometrische Interpretation der beiden Gleichungen fassen wir wieder — 


zusammen als 
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Satz: Bei der vermittelnden und bedingten Ausgleichung nach der Methode 
der konjugierten Gradienten bewegt sich der Punkt x stindig auf der H yperkugel 
tiber f — 1 als Durchmesser. 


(Bei homogenen Bedingungsgleichungen ist / = 0.) 
Zwei Bemerkungen zur vermittelnden Ausgleichung wiiren noch anzufiigen : 
1. Im tiblichen Falle f = 0 fiihrt (39) wegen 


, 


(u;, v;) + (2, v;) = (vj + J, v!) = (x), v!) =0 
auf die bekannte Schlusspr obe 
(v, v) = —(l, v) , (40) 


welche aber hier schon im Verlaufe der Rechnung standig erfiillt sein muss. 
2. Nach dem Satz von PytHacoras folgt aus (38) : 


(vo; 0) = (=Ef-)—(@, -—hm—f.- 


Der letzte Term dieser Gleichung wachst aber infolge der Orthogonalitaét der 
Ax; monoton. Somit ware nochmals verifiziert, dass (v;, v;) monoton abnimmt. 
Des weiteren soll untersucht werden, was bei der Anwendung des cg-Ver- 
fahrens aus den beiden Abschatzungen fiir (v, v) wird: 
a) Grobe Abschaétzung: Wird vermittelnd ausgeglichen und von der letzten 
Naherung x; der Ubergang zu einem x” gemacht, so gilt mit v = s—/ nach 


Gleichung (20): 


(vf, v!) — (a" — xi, x" — a1) Sv, 0) S (0, 0). 

Die untere Schranke kann nicht durch (x” — 1, x" — 1) ersetzt werden, da dieser 

Wert ohne weiteres grésser als (v, v) sein kann (siehe erstes Beispiel im 4. Teil). 
Im Falle der bedingten Ausgleichung bis x/ mit Ubergang zu einem x’ 

wird aus Gleichung (20): 


(v’, o') — (x; — 2’, 4, —x) S50, 9S (v’, 0’). 
Die linke Seite der ersten Ungleichung ist aber wegen (36) gleich 


1 Ui 
(4 1,4; = 1) = (0, 8) , 
also darf man hier ersetzen. 
Zusammenfassend gilt somit folgender 


Satz: Es werde nach dem Verfahren der konjugierten Gradienten ausgeglichen 
und vom letzten Naherungspunkt (x; bzw. x!) aus der Ubergang zwm Orthogonal- 
vaume (x" bzw. x’) gemacht. Im Falle der vermittelnden Ausgleichung lefert x; 
die obere Schranke fiir (v, v), die untere ergibt sich aus der groben Abschatzun gs- 
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es. mit x’ und x". Bei der bedingten Ausgleichung gibt x; die untere Schranke, 
x' die obere; die Abschatzungsformel liefert keine weitere V erbesserung. 


b) Feinere Abschitzwng: Fiir die feinere Abschatzung mit dem Hyperkreis 
hatten wir in Abschnitt 3.3 orthogonale Vektoren z/ und z/’ vorausgesetzt. Solche 
Orthogonalsysteme liefert nun zwar das.cg-Verfahren in der Gestalt der Ax; 


| 


bzw. Ax, jedoch stellen wir leider fest, dass gerade infolge der Orthogonalitat | 


der Korrekturvektoren die u; mit den x; bzw. die u;’ mit den x/’ zusammenfallen 
und damit Gleichung (31) keine engeren Schranken als die schon bekannten 
(vj, vj) und (v/, v/) ergibt. Das Verfahren bester Strategie liefert in diesem 


Sinne nicht gratis noch zusatzliche Information: 


Satz: Es werde gleichzeitig vermittelnd und bedingt nach dem Verfahren der 
konjugierten Gradienten ausgeglichen. Die Korrekturvektoren Ax; bzw. Ax; Sig 


nen Tetlebenen G’ bzw. G" von F’ bzw. F" auf. Die Fusspunkte u' bzw. u" der — 


kiirzesten Verbindung zwischen G’ und G" fallen dann mit den letzten Naherungs- 
punkten x; und x; zusammen. 


So gibt zwar im Falle der vermittelnden Ausgleichung mit Ubergang von 
x; za einem x” (32) die beste untere Schranke, hingegen ist bei der bedingten 
Ausgleichung mit (32) nichts zu gewinnen. 


Der Beweis fiir das Zusammenfallen von u; und w/ mit x; und x’ folgt leicht — 


aus den bisherigen Ausfiihrungen. Mit ganz entsprechender Begrindung v wie 
fiir (36) gilt auch: 


(4; x, Ane) 0 Tot 9 eee, 2a ee (41) 


wobei Ax/_, = x; — xj_,. Da aber die Teilebene G" von F" (siehe Abschnitt 3.3) . 
durch die Ax! aufgespannt werden soll, ist x,’ nach (41) die Projektion © 


von x’ auf G” und damit x; = uj’. Ebenso gilt dual xj = ui. 

Diese Koinzidenzen kénnen immerhin fiir eine zusatzliche Rechenkontrolle 
ausgentitzt werden: Nach (26) gilt in jedem Schritt (unter Weglassung der 
Indizes): 

Aidt soe a ae 


(Zaeeca) oe 


Da aber 2’ = Ax’ = Aw’, wird bei Beniitzung von (13): 


jee a) ee Ge Es ad 
(2, 2’) (a’ q’, A’ q’) WAGs: Ge = ee 
oder 
C-fh@)=(r.r), (42) 


und entsprechend bei der bedingten Ausgleichung: 


—(l ee i gq") == % r") ; (43) 
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Lees STO ND es 


Falls schliesslich beim vermittelnden Ausgleichen die y mitgefiihrt werden, 
diirfte noch folgende Rechenprobe willkommen sein: 


(ro ¥) = (CTL yy = (, Cy) = (, x’). (44) 


Entsprechend gilt in der bedingten Ausgleichung fiir die Korrelaten ft: 


— (6, 1) = (Bf, t) = (7, BT t) = (f, 0") = («", 0"). (45) 


Fihren wir die Widerspriiche w = B/ + 6 ein, so folgt aus (45) durch Subtrak- 
tion von (B/, t) = (/, v") wieder eine bekannte Schlussprobe: 


—(w, t) = (v", v"), (46) 


welche infolge der Orthogonalitat der Korrekturvektoren beim cg ebenfalls 
schon in jedem Schritt erfiillt sein muss. 


Zusammenstellung der Rechenproben 


Vermittelnde Ausgleichung | Bedingte Ausgleichung 
(9) =O (@ #}) (33) 
(nr) =0 (#)) (34) 
(x; ,0;).= 0 (39) 
(Ag) = 7) (42) —0-f4/) = 1) (43) 
(75, 93) = (l, x) (44) ~ (6,4) = (xt, of) (45) 
| =@.t) “= @f, 0) (46) 


4. Beispiele 
4.1 Nivellementsausgleichung (vermittelnd) 


Alle drei Beispiele wurden vom Verfasser auf der elektronischen Rechen- 
maschine der ETH (ERMETH) nach dem im 2. Teil beschriebenen Verfahren 
der konjugierten Gradienten durchgerechnet. Dabei wurde von allen im Ab- 
schnitt 3.5 angegebenen Rechenkontrollen Gebrauch geinacht. : 

Im ersten Beispiel werden Nivellementsmessungen vermittelnd ausgegli- 


chen. Es seien sechs Punkte A bis F so im Gelande verteilt, dass infolge der 


Sichtverhaltnisse nur 12 Héhendifferenzen zwischen je zwei Punkten direkt 
gemessen werden kénnen (Fig. 4). (Wenn alle méglichen 15 Differenzen in die 
Rechnung eingehen, wird die Sache zu «regelmassig», das heisst, die Eigen- 
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werte der Matrix C? C fallen bis auf einen zusammen, und der cg-Algorithmus 
fiihrt schon in zwei Schritten zur exakten Lésung, was bei unserem Beispiel 
nicht erwiinscht ist.) 


Figur 4 


Beispiel einer Nivellementsausgleichung. 


Als Unbekannte wahlen wir die Héhen 7 bis 5 der tibrigen fiinf Punkte 
iiber A, welches der Anschlusspunkt zu einem anderen Netz sei. 
Die Beziehung x = C y zwischen den 12 Messgréssen x und den 5 Unbekannten 
y wird hier sehr einfach. Die Matrix C heisst, unter Berticksichtigung der in 
Figur 4 durch Pfeile angedeuteten Orientierung der Differenzen: 


1 0 0 0 0 
0 1 0 0 0 

0 0 1 0 0 

0 0 0 if 0 

0 0 0 0 1 

Sal 1 0 0 0 
Caled 0 1 0 0 
0 al ia 0 0 

0 =H 0 1 0 

0 0 | af 0 

0 0 =f 0 i 

0 0 0 =a) 1 


(Praktisch wiirde man nicht mit den Hdhen selbst, sondern mit Korrekturen 


rechnen, ferner waren die Fehlergleichungen mit verschiedenen Gewichten zu 
versehen.) 


Fiir unser Zahlenbeispiel nehmen wir: 
P=ii —-2,5 0 = 34 0,7 =A 6 152 Gre aaa eae — 8,2}. 


Die exakte Lésung 


ll 


p= i= {—13, 0.1 0,905 20% 
Hoe $= {— 1,300.1, 0,9..0,5. ~2,9. 14022) 0 ,5icl0;A a0 dues, BOs 


wurde nach 5 cg-Schritten erreicht. 
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Fir 4; = (v;, v;) ergaben sich folgende Werte (gerundet) : 


Schritt 7 Ne 

0 (Start) 266,190 = ((Z,.2) 

1 224,277 

2 221,874 

5 221,390 

4 221,285 

i 221,260 (exakt) | 


Diese Zahlen liefern obere Schranken fiir (v, v) = 75. Fiir den Ubergang zu 
Punkten x” aus F” hat man zu beachten, dass diese die Gleichung 


CT (x" —1) =CT vo" =0 (47) 


befriedigen miissen. Die Berechnung von x” wird dank dem speziellen Bau der 
Matrix C hier trivial: Wenn die letzten sieben Komponenten von x; (cg-Nahe- 
rungspunkt) tibernommen werden, ergeben sich die ersten fiinf Komponenten 
von x" direkt aus (47), ohne dass ein Gleichungssystem aufzulésen ist. Im 
dritten und vierten Schritt erhalt man so: 


Cig ety =") = 0,282) - (a — 0", = "= 0,087 


und damit nach (20) die folgenden Schranken: 


3. Schritt: 221,108 < (v, v) < 221,390 
4. Schritt: 221,198 < (v, v) < 221,285 
221,260 


Der selbe Punkt x” gibt im 4. Schritt mit der feineren Abschatzungsformel (32) 
die bessere untere Schranke 221,228. 

An diesem Beispiel ist ersichtlich, dass der durch den Ubergang gewonnene 
Punkt x” allein keine untere Schranke liefert, denn im vierten Schritt ware 
zum Beispiel (x” — 1, x” — 1) = 226, also viel zu gross. 


4.2 Triangulierungsnetz (bedingt) 


Das zweite Beispiel stammt aus dem bekannten Handbuch von JORDAN- 
EccERT [8] (Beispiel aus der Badischen Triangulation). Es handelt sich um 
bedingte Richtungsausgleichung in einem Triangulierungsnetz, und zwar 
wurden alle 12 Richtungen in einem Viereck gemessen (Figur 5). 
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Man iiberlegt sich leicht, dass das Problem acht Freiheitsgrade besitzt : in 
jedem Punkt legt eine Richtung die Orientierung des betreffenden Satzes fest, 
und durch vier weitere Richtungen ist die ganze Konfiguration eindeutig 
bestimmt. Die 12 Variabeln haben somit vier Bedingungsgleichungen zu er- 
fiillen. Da die Bedingungen fiir die Winkelsummen nur drei unabhangige 


1 12 


6 
Figur 5 
Beispiel einer Richtungsausgleichung. 


Gleichungen liefern, hat man noch eine Seitengleichung heranzuziehen (das — 


heisst eine Gleichung, welche trigonometrische Beziehungen beniitzt und erst 
linearisiert werden muss). Die Koeffizienten der Bedingungsgleichungen sind 
unten angegeben. Wir miissen, da die Aufgabe schon so formuliert ist, dass nur 
noch die Verbesserungen als Unbekannte auftreten, nach unserer Bezeich- 
nungsweise x = v, 1 = O setzen. 


=—3,02Z ~ 8,00) 4598) — 022, e250 275+ On Oem ome 
pelt 9 Ci eh St Ul i 0c eee 1s) 
ae 1-100 6} Ea 0 0 
= 1 0 0 0 0 es ae eee | 0 1 | 
| = 2,90 
tee 
1,02 
| 1,25 | 
Die Lésungen sind (auf drei Stellen nach dem Komma gerundet): 
t=h = {0,044 0,076 —0,229 — 0,430}, 
H=v=6S 


={0,222 0,148 -0,370 0,143 0,188 -0,331 0,229 0,201 —0,430 0,213 0,120 —0,333}. 


Die (praktisch bedeutungslosen) Abweichungen von héchstens 5 Einheiten der 
letzten Stelle gegeniiber den bei JorDAN-EGGERT [8] angegebenen Resultaten — 


rihren davon her, dass unser Rechenautomat 11stellig rechnet. 
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Tabelle der 1; = (v;, v;): 


—_————<—<—— 


Schritt 7 N; 

0 (Start) 0 (4% = 4, 85 = 0) 
x 0,109 

2 0,623 

3 0,814 

+ 0,821 


Diese Werte bilden untere Schranken fiir (v, v). Der Ubergang nach x’ geht bei 
diesem einfachen Beispiel noch ganz ohne Auflésung von Gleichungssystemen. 
Wir kénnen tiber 8 Komponenten von x’ frei verfiigen und iibernehmen daher 
die 1.4. und 9,-12. von x;'. Die restlichen Komponenten ergeben sich dann 
sukzessive durch Einsetzen der schon bekannten in die Bedingungsgleichungen, 


und zwar 
die 8. Komponente aus der 4. Gleichung, 


die 5. Komponente aus der 2. Gleichung, 
die 6. Komponente aus der 1. Gleichung, 
die 7. Komponente aus der 3. Gleichung. 


Dieser Ubergang wurde im zweiten und im dritten Schritt ausgefiihrt und 
damit je eine obere Schranke gewonnen. Die durch die 7; gegebenen unteren 
Schranken lassen sich, wie in Abschnitt 3.5 erklart, nicht verbessern, also: 


2. Schritt 0623 S(0, v¥S 2525 
3 Senritt 0,814 < (v, v) S 0,831 
0,821 


Falls man sich nach dem dritten Schritt auf Grund der relativ engen 
Schranken dazu entschliessen wiirde, die Rechnung abzubrechen, so hatte man 
tatsdcblich v-Werte, welche um héchstens 50 Einheiten der dritten Stelle nach 
dem Komma falsch sind, was in diesem Falle wohl eben noch gestattet werden 
k6nnte. 


4.3 Approximation durch Polynom (vermittelnd und bedingt) 


Hier handelt es sich um ein Schulbeispiel mit ganzzabligen Lésungen. Es 
sollen 16 vorgegebene Punkte (Ordinaten: Messvektor /) mit dquidistanten 
Abszissen durch eine Parabel 7. Grades approximiert werden. Die Ausgleichung 
wurde gleichzeitig vermittelnd und bedingt durchgefihrt : 


n = 16 Messvariable, 
y = 8 Unbekannte (Koeffizienten der Kurvengleichung) , 
B= 8 Bedingungsgleichungen (8. Differenzen = 0). 


ZAMP X/18 
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Ob man sich in einem konkreten Falle fiir die eine oder die andere Variante 
entscheidet, hangt von der Anzahl der Punkte und dem Grad des verlangten 
Polynoms ab. Ferner ist zu beachten, dass die bedingte Ausgleichung nur mit 
den Ordinaten arbeitet und keine Kurvengleichung liefert. Die Koeffizienten- 
matrix B der Bedingungsgleichungen ist einfach gebaut (Binomialkoeffizienten) 
und enthalt eine betrachtliche Anzahl Nullen. 

In die Matrix C der Fehlergleichungen der vermittelnden Ausgleichung 
hingegen gehen die Abszissen der Punkte wesentlich ein, da ja die Kurven- 
gleichung bestimmt werden soll. Die Matrix C wird im allgemeinen ausgefiillt 
sein und sich unter Umstanden sehr schlecht fiir die numerische Rechnung 
eignen. Wenn als Unbekannte die Kceffizienten der Parabelgleichung 
N= Cy + 6, & + Cy é% +--+ + 0c,_, &’! gewahlt werden, so erhalt man schon 
in unserem Beispiel mit y = 8 infolge numerischer Instabilitat tiberhaupt keine 
verniinftigen Resultate mehr. (Ich verdanke Herrn Professor RUTISHAUSER den 
Hinweis auf die Verwandtschaft der betreffenden GauBschen Matrizen C? C 
mit den sehr schlecht konditionierten Hilbert-Matrizen.) 

Eine wesentliche Verbesserung ist zu erwarten, wenn die gesuchte Funktion 
nicht nach Potenzen von &, sondern nach einem System von orthogonalen 
Polynomen entwickelt wird. Wir machten den Versuch mit den auf das Intervall 
{[— 1, + 1] normierten Legendre-Polynomen P, und wahlten die Abszissen der 
16 Punkte von — 0,7 bis 0,8 im Abstand 0,1. 

Es folgt eine Zusammenstellung der Daten, mit denen die Rechnung durch- 
gefiihrt wurde: 


P= Oe B07 P,(—0,7) 
|| P,(—0,6) P,(—0,6 P.(=0;6) 
C= 
P,(0,8)  P,(0,8) soe AP AOEB) 
lc. =O; a ei0525 Sele iia 051s iar 
|. Ly = O%6) ar GONG seni, 3206 fea 
1-05. ONZ5s ate e eG 2 ostee 
10/8" "0;4655 ae seein cso 
Pe 8°. 28°56 +702 56> arene 
f= 8* 28. 561970. 56 eae Oke 0 
{vas 1! = 8. 9 28h=56) P0neomwe een 
Oey. ee re 
1-8 28 —56 70, —56 28 —8 1 
= 


~~ ee ee 
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T= {471 402 1024 -95 1228 -46 301 -100 -331°~354 ~1018 95 =1214 78 -261° 130} 
Die exakte Lésung 

#=s={451 542 604 605 528 374 161 -80 -311 -494 -—598 -605 —514 —342 -121 110} 


sollte theoretisch nach 8 Schritten erreicht werden. 

Die tatsachlich erhaltenen Werte sind jedoch mit Rundungsfehlern behaf- 
tet, und es zeigt sich, dass zusatzliche cg-Schritte noch eine Verbesserung der 
Resultate bringen. 

Die Zwischenwerte der 4; = (v;, v,) sind, auf Ganze gerundet, in der folgen- 
den Tabelle angegeben: 


Schritt ¢ ;, (vermittelnd) n; (bedingt) 
0 (Start)| 5896874 = (i, /) Ogi ts Up 0) 
1 3038088 2005649 
2 2892854 2692796 
3 2757019 2745146 
+ 2747769 2769706 
S 2746849 2746532 
6 2746439 2745 534 
7 2745601 2745 565 
8 2745601 2745600 
9 2745600 


Dies sind zugleich die besten, auf Grund des cg-Verfahrens erhaltlichen 
Schranken fiir den exakten Wert (v, v) = 2745600. Es fallt auf, dass in Schritt 4 
und 5 das 7", im Gegensatz zur Theorie, grésser wird als der exakte (v, v)-Wert. 
Diese Stérung, die auf Rundungsfehler zuriickzufiihren ist, verschwindet 
interessanterweise spdter wieder. Um an dieser Stelle dennoch verniinftige 
Schranken zu erzwingen, wurde nachtraglich folgendermassen vorgegangen : 


Es wurden zwei Punkte x’ und «” bestimmt, mdéglichst nahe bei x4 und %4, 
welche aber, im Gegensatz zu diesen Naherungspunkten, wieder genau in F’ 
bzw. F” lagen, und mit ihnen die groben Schranken nach (20) berechnet, welche 
nun sicher korrekt sind: 


2745482 < (v, v) < 2747919. 


Als Beispiel fiir die Berechnung der Gewichtsmatrix P sollen die berech- 
neten Zwischenwerte des ersten Diagonalelementes ,, 4stellig bis zum 8. Schritt 
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angegeben werden. Was in weiteren Schritten noch dazukame, ware offenbar 
sehr fragwiirdig. Im 8. Schritt sollten beide Zahlenfolgen theoretisch denselben 
Wert erreichen. 


Schritt 4 P11 (vermittelnd) Py, (bedingt) 
0 (Start) 0 1 
1 0,1119 1 
2 0,1558 1 
3 0,1564 0,9991 
4 0,3119 0,9920 
5 0,3151 0,9900 
6 0,8093 0,9899 
| a 0,8810 0,9899 
8 0,9300 0,9898 
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Summary 


The paper begins with a short description of the theory of overdetermined 
linear systems (A usgleichsrechnung), using the notation of matrix-calculus. 

The solution of such systems by means of relaxation-methods is then discussed. 
In particular a broad description of the method of conjugate gradients is given. 
The algorithm avoids the use of the matrix of the normal equations. 

Lastly the application of the ‘method of the hypercircle’, developed by 
SYNGE and others, to this problem of linear algebra is studied, with the aim of 
finding upper and lower bounds for the sum of the square residuals — this sum 
being a minimum for the solution of the system. To get these bounds, we need 
approximate solutions obtained as intermediate values in an iterative process. 


(Eingegangen: 31. August 1958.) 
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Zur Matrizeninversion nach Gauss-Jordan 


Von HeEInz RuTisHauser, Ziirich) 


Die Aufgabe, ein lineares Gleichungssystem numerisch aufzulésen, wurde 
schon von Gauss 1826 vollstandig gelést. Die seitdem in grosser Zahl publizier- 
ten «neuen» Methoden zur Lésung linearer Gleichungssysteme und zur Inver- 
sion von Matrizen sind, soweit sie auf Elimination beruhen, alle nur mehr oder 
weniger wertvolle Varianten des GauBschen Verfahrens. Die bedeutendsten 
Varianten sind ohne Zweifel diejenigen von JORDAN?), welche die Unbekannten 
in einem einheitlichen Rechengang (das heisst ohne Riickwartseinsetzen) zu 
berechnen gestattet, sowie diejenigen von BANACHIEWICz®), welche das Auf- 
schreiben stark vereinfachen. 

Die sonst weniger verwendete Jordansche Variante (auch Gauss-Jordan- 
Verfahren genannt) hat mit dem Aufkommen der elektronischen Rechen- 
maschinen wieder an Bedeutung gewonnen. Insbesondere sind in neuester Zeit 
einige Abarten des Gauss-Jordan-Verfahrens bekannt geworden‘), die sich 
speziell zur Berechnung der inversen Matrix mit Rechenautomaten eignen, 
weil sie zu wesentlich einfacheren Rechenprogrammen fiihren als das ur- 
spriingliche GauBsche Verfahren. 

Im folgenden wird eine elementare, determinantenfreie Begriindung der in 
[4] und [5] vorgeschlagenen Variante des Gauss-Jordanschen Verfahrens ge- 
geben. 


1. Die Transformation J 


Sei A = (a;;) eine Matrix mit m Zeilen und ~ Spalten. Wir betrachten die 
lineare Substitution: 


Vi = DY 45 %} (¢ = 1, 2,..., m) () 
j=l 


mit den unabhangigen Variabeln x,, ...,%, und den abhangigen Variabeln 


a 
Wir kénnen nun in den Beziehungen (1) eine Rollenvertauschung vorneh- 


men, indem wir ein beliebiges Matrixelement a,, + 0 herausgreifen und dann 


1) Institut fiir angewandte Mathematik der ETH. 

2) Siehe etwa JoRDAN-EGGERT [3], § 35. 

3) Siehe [1] sowie auch ZuRMUHL [6] oder § 23 in [7]. 

4) Siehe Petrie [4], Ersov [2], ferner RUTISHAUSER [5]. 
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die p-te der Relationen (1) nach x, auflésen: 


Va Xj - (2) 


Durch Einsetzen von (2) in die itbrigen m — 1 Gleichungen (1) erhalten wir 


1 Api | 
= y Oss = y Ay; %; + Aig \—— Vo — y 56; 
y 5% io i ee Vo A, pq df 
oder 
Aig Aig &pj : F 
;= -++ Ui — “ta"2!) wee GSS ; 3} 
DA Ona Vp a | j re cok p) (3) 


Dies sind mit (2) zusammen wieder m lineare Beziehungen mit den unab- 
hangigen Variabeln °%,, %, ---;% 9-1) Vo: M@44 +s m_ Und den Mabhangiges 
Variabeln 1, Yo, ++. 1 Vp-p Xo Vou «++» Ym: Man kann diese Beziehungen auch 


vektoriell aufschreiben: 


yp = AW x) ; 
wobei 
x2) = (435 Xo, +++) Xq—1) Vor %q41) + 


yD a (V4) Vor +21 Vp-1 XQ Vopr s+ Ym) 


und die Elemente der Matrix A" wie folgt bestimmt sind: 


1 . . 
[ : i = = oie 
pa 
| — <p! LUG? = py fe ae 
a\!) x axe (5) 
aig ae . ° 
: firi + p, -7=q¢ ‘ys 
DPQ 
(a,,—“i2%es) firi#p, j#¢ (d). 
4nq 


Hn)» (4)q 


Der Elimination einer Unbekannten %_ aus den Gleichungen (1), genauer: ; 


der Vertauschung einer abhangigen Variabeln mit einer unabhangigen Varia- 


beln, entspricht also eine innere Transformation J,,5) der Matrix A, welche 


durch die Formeln (5) bestimmt ist. Das abgesehen von der Bedingung a,,+0 


frei wahlbare Element a,,, von dem die Transformation wesentlich abhingt, 


heisst Pivot der Transformation J, ,. 


°) Es handelt sich nicht um eine Ahnlichkeitstransformation. 


a 
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Die verschiedenen Berechnungsarten der Elemente der neuen Matrix 4 
kénnen wie folgt schematisch dargestellt werden: 


Elemente ali) mit ++), 7 +4, 
die nach (5, 6) berechnet werden. 


Element 
a) 
pq’ Zeile p, Berechnung nach (5, 8). 
nach (5, a) 
berechnet. 


— 
~ 
Ss 
<3 
3 
og 
g 
" 20 
Be 
Ys 
aq 
oO 
o 
=) 
o 
faa) 


Die Transformation J,, ist natiirlich nur im Fall a,, + 0 durchfthrbar, 
dann aber ist die Beziehung zwischen x und x") linear, es gilt namlich 


pl 


und damit 


Satz 1: Durchliiuft x = (x1, ... , Xp) ein lineares V ektorgebilde der Dimension r, 
so gilt dasselbe auch fiir den Vektor x") in (4). 
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2. Iterierte Anwendung der Transformation J 


Es ist nun naheliegend, in der Beziehung (4) weitere abhangige und unab- 
hangige Variabeln gegeneinander auszutauschen, um dann schliesslich Formeln 
zu erhalten, in denen mdglichst viele der unabhangigen Variabeln durch ab- 
hangige Variabeln ausgedriickt sind. Zu diesem Zweck muss man offenbar in 
der Matrix A™ wieder ein Pivotelement a@ + 0 auswahlen und dann die 
Transformation J,,,, auf die Matrix A” anwenden. (Das erste Pivotelement 
ist jetzt sinngemass mit a, ,, zu bezeichnen.) Die resultierende Matrix A® 
kann in gleicher Weise behandelt werden, usw. 

Damit die abhangigen und unabhangigen Variabeln systematisch ausge- 
tauscht werden, miissen die Pivots — abgesehen von der Bedingung a,, + 0 — 
so gewadhlt werden, dass sowohl die p, als auch die g, je untereinander ver- 
schieden sind. Die Pivots miissen also der «Turmbedingung» geniigen, das 
heisst, sie diirfen einander nicht bedrohen, wenn man sie als feindliche Tiirme 
auf einem (m Xn)-Schachbrett auffasst. Diese Bedingung soll verhindern, dass 
der Austausch der x und y wieder teilweise riickgangig gemacht wird. 


Wir wollen also allgemein nach Berechnung der Matrix A“) das ndchste 


Pivotelement nur unter denjenigen Zeilen und Spalten von A auswahlen, die 


nicht schon ein Pivotelement enthalten. Diese Einschrankung hat natiirlich zur — 


Folge, dass auf A‘) genau dann eine weitere Transformation J ausgeiibt wer- 
den kann, wenn auch ausserhalb der Pivotzeilen und Pivotkolonnen von A‘) 
noch mindestens ein von 0 verschiedenes Element existiert. Satz 3 gibt eine 
a-priori-Aussage iiber die Zahl der auf A hintereinander ausitibbaren Trans- 
formationen /. Zunachst gilt wegen (4): 


Satz 2: Sei A‘) die durch Anwendung der S Transformationen Joa, 
(k=1,2,...,s) aus A = A erhaltene Matrix, wobei die Indexwerte py, po, ..., Ps 


und i, Jo, --., Ys Je untereinander verschieden seien. Ist dann ¥ = (x1, ...,%p) ein 


beluebiger Vektor und y = A x = (yy, ..., Vm), 80 gilt 


A®) (§, en sisieaeg Sy) = (1, Ye) WO) Nm) ? | 
wobeli 
by FOr YA Ga, as ay OG 
Vo, tit, = gp(he dy. sy 
Vs Tit 4 sep) Dae eon 
fe 
Xqy, tire = p,(k =1,...,98). 


“nr ey oe ere 
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Satz 3: Hat die Matrix A den Rang r, so kann man genau r Transformationen 
J hintereinander auf A anwenden und erhiilt dann eine Matrix A” mit der 
Eigenschaft 


a = 0. firi + go5...p, und) f+ 4,,...,¢,, (9) 


das heisst, nach Streichen der Pivotzeilen und -spalten der Matrix A”) verbleibt 
nur noch die (m—r) X (n — r)-reihige Nullmatrix. 


Beweis: Hat A den Rang 1, so existiert ein (n — r)-dimensionales Vektor- 
gebilde 2 im Raum der x, ..., *,, so dass Ax =0 fiir alle xe &. Nach Satz 2 
ist, falls man yr Transformationen ausfiihren kann: 


A” (é1, seey é.) — (M1, teey fan ) 
wobei 
Men TUES FO), Ook Sys 


Es= | 
a See aM eae ae (10) 


Ota! Pepi. PP, , 
ke ee ea bee F) 


Ae 


ee 


und x ein beliebiger Vektor aus & ist. 

Nach Satz 1 muss aber der Vektor (é,, ...,€,) ein (m — r)-dimensionales 
Vektorgebilde durchlaufen, wenn x 2 durchlauft, das heisst, die Gréssen 
x; (7 = 4, ---,Q,) unterliegen in der Gleichung (10) keinerlei Einschrankung. 
Daraus folgt aber die behauptete Eigenschaft von A”. 

Wenn man anderseits nach 7 Transformationsschritten J zu einer Matrix 
A” mit der Eigenschaft (9) kommt, so dass kein der Turmbedingung geniigen- 
des Pivotelement a,, + 0 mehr gefunden werden kann, so kann man der- 
artige Beziehungen (10) konstruieren, die fiir beliebige x, (7 + 91,2, --- Jp) @I- 
fiillt sind. Daraus folgt dann umgekehrt wie oben, dass A ein Nullgebilde der 
Dimension n — ry und damit den Rang 7 haben muss, was zu beweisen war. 

Es ist genau zu beachten, dass die Aussage von Satz 3 nicht etwa nur bei 
einer besonders giinstigen, aber praktisch nie erreichbaren Wahl der Pivots 
realisiert werden kann. Vielmehr geniigt es vollauf, als nachsten Pivot jeweils 
irgendein von Null verschiedenes Element zu wiahlen, das beziiglich der bis- 
herigen Pivots die Turmbedingung erfiillt. 

Dagegen ist es dusserst schwierig, mit Hilfe von Satz 3 den Rang einer 
Matrix tatsdchlich numerisch zu bestimmen, weil infolge der Rundungsfehler 
Nullen in der Regel als kleine Zahlen erscheinen und daher kaum von den 
«echten» Nicht-Nullen unterschieden werden kénnen. 
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3. Inversion von Matrizen 


Sei nun A eine n-reihige quadratische Matrix. Falls sie nicht singular ist, 
kann man nach Satz 3” Transformationen J anwenden; von den n Pivots 
steht dann genau einer in jeder Zeile und jeder Spalte, so dass unter den 
Indexfolgen ,, po, ---, Pn UN 4, Jo, --+> Yn jede der Zahlen 1 bis m genau je 
einmal vorkommt. Damit hat die nach m Schritten resultierende Matrix A” 
nach Satz 2 die Eigenschaft, dass aus 


Ag=y AME... on) = (hy, 2s Na) 
folgt, wobei 
ea (11) 
S7= Ve Ul Ge 
(k= 1,2, .,#)2 


Ngee teal iat =D, 


Damit ist A”) im wesentlichen die Inverse von A und kann durch blosse 
Zeilen- und Spaltenvertauschungen in diese iiberfiihrt werden. 

Zuniachst ist die Beziehung zwischen den Vektoren (7, ... , 7,) und x gemass 
(11) eine blosse Permutation der Komponenten, namlich 


(1, Ya, oc) Nn) => [Pe 


wobeli 
1 fur 2 = 4, 1 = Tk >» 12) 
bis = 
: Q in allen anderen Fallen. 
Ebenso 
(&1, een) a Qy ’ 
wobei 
(13) 


1 firt= 4,3 7= Px, 
OY ae 
/ 0 in allen anderen Fallen. 
Damit ist Q offenbar die Transponierte von P, und es wird insgesamt 


t= P+ (yy, ..., Yq) = PY AWE, ig Epes PEA Pee 


Schliesslich beniitzen wir noch, dass bei Permutationsmatrizen P-1 = PT ist, 
und erhalten damit 


Antes PT Ary pr, (14) 


—- 


ANS PUG. A yeh 4h et ~eet 
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Die hier auftretende Matrix P hangt aufs engste mit den Pivotelementen 
der bei der Berechnung von A” ausgefiihrten Transformationen J zusammen: 
P hat Einsen an den Stellen i = £,, 7 = g,, also genau dort, wo sich die Pivots 
der m Transformationen J, befinden; alle andern Elemente von P sind Nullen. 
Numerisches Beispiel (die Pivots sind fett gesetzt): 


20 Uh a = et ie: 4 eT 
(Ort WeOwL Gy iat (ay 4 (0 L006)" 
A ge ee eee etl: 
iajaee. 19) 3t md! re mas 
OES ae ae ee 
eres (am) 2) ese 38 8 4(3) _ (2) 4°36 ~8 
Jsa( ) | Gane T ih A Jis(A ) = Sy ee elite 
peal RE ag | =f ls es 
| 423: 244 "| 
(4) __ (3)) _ spe, 45 
A elelAl Neh oe 3 14 ait | 
eon 1 3. 3 || 
Hier ist offenbar 
| 0 al 0 o | 
or ee eee | 
SS ge 1 0 
1 122 06) on 
und damit nach (14) 
0001 84 3) S 1dr 0 0001 1. eee 
ga—|| 1000 es ene re 1000 23 10st 
weil Be aes Gt EST Meet 0010 3 ieietaeeG He 
| 
A OO gi fhr a Toot 3—3|| |lo100 heed 4G eA S| 


In der Praxis wird man freilich A~ nicht durch effektive Matrizenmulti- 
plikation nach (14), sondern durch Zeilen- und Spaltenvertauschungen be- 
rechnen: 

Die p,-te Zeile von A™ ist die q,-te Zeile von A“, 
und (15) 
die g,-te Spalte von A™ ist die p,-te Spalte von A7?. 


Man wird zweckmiissig nach dieser Regel zuerst die Zeilen umordnen und dann 
die Spalten entsprechend behandeln. 


4. Summenkontrollen 


Sei S eine (m x n)-Matrix mit der Eigenschaft, dass die Summe aller Ele- 
mente in jeder Zeile +1 betragt, also 


> Su= +1 =1,2,...,m)- (16) 
j=l 
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Hieriiber gilt: 


Satz 4: Die Eigenschaft (16) ist invariant gegeniiber der T: ransformation y 
das heisst mit S haben auch die daraus abgeleiteten Matrizen S"), S) usw. diese 


Eigenschaft. © 


Beweis: (16) heisst, dass x = (1,1, ... ,1) den Bildvektor y =Sz=(1,1,..., 1) 
besitzt, das heisst, fiir x;= 1 ist y;=1. Nach Satz 2 ist also auch 


SOUS) = ae 


was zu beweisen war. Entsprechendes gilt fiir Matrizen S, deren Kolonnen- 
summen —1 sind, das heisst fiir Matrizen S mit der Eigenschaft 


ne ol GF 11 De. re PN (17) 


Satz 5: Die Eigenschaft (17) ist invariant gegeniiber der Transformation J. 


Beweis: Die Eigenschaft (17) aussert sich darin, dass fiir einen beliebigen 
Vektor x und sein Bild y die Beziehung 


n m 
dD, tet D5 = 0 
iat j= 


gilt. . 

Anderseits zieht die Giiltigkeit dieser Beziehung fiir alle x auch die Eigen- 
schaft (17) nach sich. Es bleibt also nur noch zu zeigen, dass auch fiir die Vek- 
toren x" und y" aus (4) dieselbe Beziehung gilt. Dies ist aber in der Tat der 


Fall, denn das Vektorpaar x"), y') hat bis auf zyklische Vertauschung die- 


selben Komponenten wie das Paar x, n, was zu beweisen war. 

Aus diesen beiden Satzen lassen sich verschiedene Rechenkontrollen fiir die 
praktische Durchfiihrung der Matrizeninversion ableiten. Sie beruhen darauf, 
dass man der zu invertierenden Matrix noch eine weitere Spalte bzw. Zeile 
derart zufiigt, dass eine Matrix mit der Eigenschaft (16) bzw. (17) entsteht 
und auf diese Matrix S m Transformationsschritte J anwendet. Es kann dann 
sowohl bei S™ (das ist die mit einer Kontrollsummenspalte bzw. -zeile ver- 


sehene Inverse) als auch bei allen Zwischenmatrizen S“ die Eigenschaft (16) 


bzw. (17) gepriift werden. 


5. Zur praktischen Durchfiihrung der Matrizeninversion 


Im numerischen Beispiel von Abschnitt 3 wurden der Einfachheit halber 
immer Einsen als Pivotelemente gewahlt. Da aber durch die Pivotelemente 


dividiert werden muss, wahlt man in der Praxis zweckmassig méglichst grosse — 


ae 
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Elemente als Pivots. Das kann dadurch erréicht werden, dass man jeweils 
unter den als Pivots in Frage kommenden Elementen das absolut grosste aus- 


wahlt. Dann verlauft die Inverson der Matrix aus Abschnitt 3 wie folgt: 


20 10 4 ‘ 
fpr 6 3 1 
4 3 2 1 
1 1 1 1 
0,05 = 0,5 = 0,2 — 0,05 
46.4 835 1 1 0,5 
0,2 1 1,2 0,8 
0,05 0,5 0,8 0,95 
|| 0,083333 — 0,333333 — 0,166667 0,083333 
42) — || 0.333333 0,166667 0,833333 — 0,166667 
: — 0,166667 — 0,833333 0,833333 — 0,666667 
— 0,083334 — 0,166666 0,666667 0,416666 
0,100000 — 0,300000 — 0,300000 0,199999 
4) — 0,300000 0,100000 1,100001 —0,400001 
? — 0,300002 — 1,099999 1,900002 — 1,600003 
0,200002 0,399999 — 1,600003 2,400003 
1,000000 — 3,000000 3,000003 — 1,000004 
4(4) — || —3,000000 10,000000 — 11,000010 4,000010 
= 2,999995 — 10,999990 14,000002 , —6,000010 || » 
— 0,999995 3,999990 — 5,999996 4,000003 


Da samtliche Pivots in der Diagonalen liegen, ist P die Einheitsmatrix und 
A) bereits die gesuchte Inverse. Dieser Fall muss bei positiv definiten sym- 
metrischen Matrizen immer eintreten, wenn man die Pivotelemente gemass £) 
bestimmt. 

Fiihrt man die Inversion mit einer elektronischen Rechenmaschine durch, 
so wird man zweckmiassig wie folgt vorgehen: 

a) Zunachst fiigt man eine Kontrollsummenzeile an, denn es ist unbedingt 
notwendig, Kontrollen in das Rechenprogramm einzubauen. 

B) Als nachstes Pivotelement wahlt man jeweils das absolut grésste der 
noch zur Verfiigung stehenden Matrixelemente. 

y) Um das Absuchen der Matrix nach dem gréssten Pivot zu erleichtern, 
werden die Pivotzeilen und -spalten dadurch gekennzeichnet, dass man ihre 
Elemente im Speicherwerk der Maschine mit einem besonderen Kennzeichen 
versieht (Q-Zeichen bei der ERMETH). Die Summenzeile wird schon beim Be- 
ginn markiert, da sie keine Pivots stellen darf. 
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6) Die Transformationsformeln (5) werden zweckmissig in leicht variierter 
Form angewendet: Man bilde zunachst die Gréssen 


oi fir 7G +], 
ee Inq 
j= 
: fir 37. = 9 
ang (18) 
asa fiir t == p , 
C; == 
il fi 7 


und lésche dann die #-te Zeile und g-te Spalte der Matrix aus. Dann wird, 
wenn wir die Elemente der derart geléschten Matrix mit aj, bezeichnen: 


al) - a, +c, 6, fiir alle 7 und 4. (19) 


Entsprechend werden die weiteren Transformationen J, , modifiziert. Auf 


diese Weise kann eine Vereinfachung des Programms und eine Beschleunigung 
der Berechnung erzielt werden. 

e) Nach jedem Transformationsschritt werden die Kolonnen- (Zeilen-) 
Summen kontrolliert. 

y) Damit am Schluss die zur Ausfiihrung der Zeilen- und Spaltenver- 
tauschung notwendige Information zur Verfiigung steht, speichert man nach 
jeder Transformation J,, die Zahl p in der Zelle N + q und die Zahl g in der 
Zelle M + ; die Zahlen in den Zellen M+1bisM+nundN+1bisN+n 
bestimmen dann die auszufiihrenden Vertauschungsoperationen (15) voll- 
standig. 
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Summary 


‘The so-called ‘Gauss-Jordan’ procedure for inverting matrices has been de- 
scribed in many papers and is well known. The present paper gives an unusual ap- 
proach to this method insofar as the inverse is obtained by repeated application of 
a certain non-similarity transformation to the original matrix. Of course this 
transformation is related to the Gauss elimination process, but formally has rather 
the character of an exchange of components in two vector spaces. 


(Eingegangen: 9. September 1958.) 


Friction and Heat Transfer of Compressible Flow 
into an Infinite Lattice of Flat Plates 
By SuAo L. S00, Princeton, N..]., USA) 


Introduction 


The problem of friction and heat transfer of compressible flow into a lattice 
of turbine or compressor blades and flow into the plate-fin arrangement of a 
heat exchanger, particularly that of interrupted fin construction cannot, in 
general, be approximated by existing solutions of flow over a body in an in- 
finite medium. The difference between the case of a grid of straight passages 
and a single semi-infinite flat plate is obvious: First, flow cannot proceed 
through this type of channel without pressure gradient. Second, the normal 
component (to the surface) of velocity must be zero at the center of each chan- 
nel. Third, transition can occur in two ways in a channel, to develop into 
turbulent boundary layer or into Poiseuille type of flow. 

In this analysis, an idealized system of an infinite lattice of flat plates o 
zero thickness was taken to represent the above applications. The solution 
presented in the following represents a generalization of the flat plate solutions 
developed by Brastus, Crocco, CHAPMAN, and RuseEsIN (1, 2, 3]?). Consistent 
with applications of practical interest, two cases are considered in this study: 
Case 7: constant wall temperature (approximating the condition in a heat ex- 
changer) ; case 2: adiabatic wall (approximating the condition of uncooled blade 
passages). Both cases involve uniform and steady velocity and temperature of 
the incoming stream. The method can be applied, however, to the case of 
steady arbitrary surface temperature distribution. 


1) Department of Mechanical Engineering, Princeton University. 
2) Numbers in brackets refer to References, page 309. 
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Studies made to date on the subject of flow through channels are presented 
in [4, 5]. Reference [4] considers the entrance condition, but with prescribed 
velocity profiles, [5] deals with compressible fluid, but neglecting entrance 
effect and heat transfer. Other studies on thermal entrance region are listed 
in [6-10]. 


Notations 
C Constant defined by equation (11); 
ce drag coefficient ; 
Cy specific heat at constant pressure ; 
Pipdatesas functions of velocity distribution ; 
Go, G,,... functions of temperature distribution ; 
Gt defined by equation (19); 
h enthalpy; 
k thermal conductivity ; 
K,, K,,... constants defined by equation (16); 
M free stream Mach number; 
M,, M,,... constants defined by equation (18) ; 
Nuz Nusselt number based on characteristic dimension 4%; 
p pressure ; 
Pr Prandtl number; 
Q heat transfer per unit width; 
q heat transfer per unit area; 
R gas constant; 
y recovery factor [3]; 
Re, Reynolds number based on gap s as characteristic dimension; 
Re, Reynolds number based on x as characteristic dimension; 
s gap between plates; 
P static temperature, absolute; 
Wh total temperature, absolute; 
u axial component of velocity; 
v normal component of velocity; 
% axial coordinate; 
y normal coordinate; 
y ratio of specific heats; 
n transformed normal coordinate; 
Lt viscosity ; 
é transformed axial coordinate; 
0 density ; 
co shear stress; 
y stream function; 
2 function defined by equation (A 20). 
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Subscripts 
0 Refers to inlet free stream condition; 
a refers to condition at the wall; 
2 refers to condition of the core. 
Superscripts 
Le iba Dimensionless quantities as defined; 
NE total derivatives with respect to ». 


Basic Equations and Assumptions 


The cases described in the introduction are depicted in Figure 1. 


Ay y 
| ler 
mee ome f “4 [3 3-9 

2 . of aghs Lo zi 

I Ih 

Po Po 
Case I: Constant wall temperature Case II: Insulated wall 

Figure 1 


Diagrammatic representation of the flow system. 


293 


According to the classical laminar boundary layer theory with pressure 
gradient only in the direction of flow, the continuity equation, momentum 


equation, and energy equation of this flow system are: 


Ox +r oy Os 
oleae tee] -— a tw HR): 


Nae ea cele ga er CA 


The equation of state of fluid is given by 
p= oT. 


= ZAMP X/19 
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The boundary conditions are: For both cases, 


u(x,0)=O= u(x, s), u(0,y) = M%, | 
: 5) 
Ou | oT : 
—— — — i 0) == Te , ee = 0 , 
oy ly =s/2 OF jonas 09 US; ( ,¥) 0 oy y =s/2 | 
and, for case I, constant wall temperature, 
D(x, 0) Te Ts) (6) 
for case II, adiabatic wall, 
oT 
eo Oe (7) 
y=s 
As long as a core exists, with velocity of core ™: 
ls du, 
of aces (8) 
For constant total enthalpy [11], 
d ad 
Cita (Cp Ty) = ce 
or 
Cy aT, ep 
Rely ode eh wade (9) 
Further, for the isentropic core, 
o(x, y) T(x, y) 01(¥) Ty(*) 
Co) nr 7 0) 
It is assumed also that 
C,=const and we=C ty Q (11) 
in general [3], and at the wall, 
k 
Es 2 ee ae (12) 


Ho ko 


Solution of the Problem 


The details of the solution based on Howarth transformation [12] and — 


series expansion from the leading edge is given in the Appendix. 


~— Same, or 
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-¥/ (9) “Y (q) */ (v) uoTyenbea uNzUsUIOUT 94} JO UOTN[OS 
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For the coordinates x and y expressed in terms of 


yo2ee| cal 


r) ce i As X Uo es. me ch 
Up Oo S” Sucre 


ge? 


and 


the velocity distribution is given by [1] 


Ug 


The velocity of the core is given by 


Ug 


th = Ll) + EKO) + BA) + BAC) + 


“1 = ut =14K,6+K, 2+ K,84+ Ky b+ 


ZAMP 


(15) 


(16) 


where the K’s are constants. These constants are determined to be K, = 3-44156, 
K, = —8-951536, K, = 158-3782, K, = — 1669-602, by series expansion from 
yn = 0 and asymptotic expansion from 7 = oo. The function /,(7) is given in [1] 
and [3]. fs, fz, and f/, are presented in Figure 2 respectively. In Figures 3 and 4, 


the velocity distribution for various € are presented. 


13 


11 


a|S 
10 
0 002 004 006 0:08 0-10 
—— 
Figure 3 
Velocity distribution in the core versus €. Velocity distribution normal to the plate at 


=e. 6 oo -* 
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“SL:0 = 4d “*5 (p) 
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The temperature distribution is given by 


= — T+ = (y — 1) M? [G,(n) + & Gi(q) + & Ga(m) + & Gal) + -} 07) 
and 


FA, Tt 214M E4+M,F+MEe+ (18) 


where M is the Mach number of the inlet stream and M,, M,, etc. are constants 
dependent on Ky, Kg, etc. in the above. Figure 6 gives the values of 7,/T7p. 


10 cz 


09 


08 


| 


(y-1)M7- 


WwW 


0 0-0 004 006 008 010 


Figure 6 
Temperature distribution in the core of various inlet Mach numbers M. y = 1:4. 


For the case of the insulated walls (case IT) the functions G(n)’s are basically 


independent of 7, and M particularly when G, is replaced by G* according to ~ 


it 


CO Care ye Ian 


(19) 
For the case of the insulated walls, G¥, G,, Gz, etc. are given in Figure 5, 
for k = 1-40 and Pr = 0-72. G¥ is identical to 2 7(y) of reference [3]. 
For the case of walls at constant temperature, the method presented in 


the Appendix is valid for any steady temperature distribution. The solution — 
is dependent on both T,,/T) and M. The case of uniform wall temperature was — 


calculated. Several combinations of (y — 1) M? and T,,/T, were calculated. 
The choice of parameter was made for the convenience of logarithmic inter- 
polation. (1) 7,,/T, = 0-3, (y—1) M%=0-3 (Figure 7), (2) T,/Ty)=3 and 
(vy — 1) M? = 10 (Figure 8), which show distinctly different trends. The tem- 
perature distribution in the latter case shows heat flow into the plate and the 


stream due to dissipation. Figure 9 shows the values of G/(0) at various — 


T,,[Ty and (y — 1) M®. This set of data has to do with heat transfer from the 


— a 
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(c) (d) 
Figure 8 


Solution of the energy equation for constant temperature plate 
T/T) = 3, (y—1) M? = 10, Pr = 0-72. 
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surface. Other cases will be presented and discussed in the following in terms 
of characteristics of friction and heat transfer, 

From these results, the fluid friction and heat transfer of this system of 
flat plates can be determined. 

All the numerical solutions of the ordinary differential equations were 
obtained with an IBM 704 digital computer. 


3 - nn 


=} 


Figure 9 


G},(0) for constant temperature plate. 
G’(0) = 0 at T,,/T, = 1-15 for (y— 1) M? = 0:3, G’(0) = 0 at T/T) = 1-43 for (y— 1) TE = 


Friction and Heat Transfer Characteristics 


The shear stress at the wall is given by 


Ou Uy Mo N 2 4! Pays 
Te = He gy = gee Re or Hi (0) + EN (0) + 8/500) 
where 

Re,= ne 

0 
The drag coefficient at x* is given by 
sm ast 
Fe! 

Ca at . 

“2 80 Uy 


(21) 
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For constant surface temperature, 


co= VCS [KO (Fe) + 2 KO (Fee) +5 BO (ee) (22) 
Le eign t. 


2 
The friction factor calculated is shown in Figure 10, plotted is Re, (= x ug 0 9/Mo) 
for various values of Re, together with the values for single flat plate (Re, ~oo) 
given in [3] for comparison. It is shown that the value of Re, affects the friction 
factor at the wall. At low Re,, the friction factor is greater than that of a single 


flat plate. 
002 10 4 ATI = 
<| Re, =10 2 
0-01 = L 
14 Re.=10 SS at 
| 
fi ESS et 
0-005 FH ie PSs =< = a 
ice fel < Also curves of Cy versus 
tal de | ||| 7S Res (Sion same scale 
Reference(s! <~ 
0002 | be | single plate 
NI 
: i 
0-001 ; Lt 
10 2 5 10° 3 5 10° 
Re, SS 


Figure 10 
Friction drag at various Re, as compared to that of a single plate. 


Due to expansion of the fluid at the inlet to the plates, the temperature 
variation along the hypothetical insulated wall is of interest. The wall temper- 
ature can be evaluated from 


= =1+ (y—1) M?[G#(0) + € G,(0) + & G0) +--+]. (23) 


Of greater interest is the relation between the wall temperature and the stagna- 
tion temperature of the free stream, or 
e 


fe 
(ies wPN 
7K a 1 y—1 ° 3 (24) 
i 2 M 


gud? mad) ee ig, Se) ee ee 
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For various Mach numbers of the inlet, T/T; was calculated for various & 
and M as shown in Figure 11. This shows that besides the reduced recovery 
at high speeds, the decrease in wall temperature from the leading edge means 
lower leading edge temperature when thermal conductivity of a real wall is’ 
taken into account. 


1:0 


Figure 11 


Temperature distribution along the wall of an insulated plate. 


The rate of heat transfer by convection from the constant temperature wall 
can be determined at every point along the wall according to 


ae D 
g=—he |, = hw (Tw Te) : 
he Ty ey 
= — fe 72_ Ral? gt (y —1) MP (G5(0) + £ G0) + GO) + =, 
SVX 
where 


Ty = T, (1+ 75~ mM) 


is the total temperature of the stream. The local Nusselt number at x is 


| =~ (y—1) M? 
Nu, = Ny % bs ee ore Re, (y ) To 


‘fe 2c (Z — To) (26) 
x [Go(0) + & G;(0) + & G,(0) ++]. 
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From the leading edge to x, the total heat transfer rate can be calculated from 
for constant surface temperature: 


jp 


g=— 6h yey 1) mt [e Gyo) +S 


ie 


3 
G,0) + 5 G0) ++]. 27 
The values of G/(0) are given in Figure 9. 
The ratio of heat transfer and skin friction is obtained by 
q _ Aho Ty «(Y= 1) MP 1Gs) 8 G0) ar aa) 


ihgiabe RI. (7,(0) + & 73(0) ---] 


The ratio of total rate of heat transfer and the friction drag of length x of 
plate is given by : 


= wee (y —1) M2 [£500 sa pase ee | 


a) 


Ca 


¥ [e700 + 5 180) -] 


Ss 
x 


The condition given in [3] applies only to the leading edge. © 
The Nusselt number over length L from the leading edge is: 


po 1 M2 
hy L 1 2 


ape Lie oe EE ; oe 1 ue] 
0 0 


Nu, = 


x P Re; [650) + | G0) + ---]. | 


The calculated values are plotted in Figure 12, also shown are data of single 
flat plate case given in [3] and [15] for comparison. 


Discussion 


From the above consideration, flow into a row of parallel flat plates may 
be visualized as flow over a single flat plate with prescribed pressure gradient 
as that in the core in the above. The variation of static pressure in the core in 
the direction of flow depends on the distance from the leading edge and the 
Reynolds number of the gap. As shown in Figure 6, the pressure drop seems to — 
be large, however, this isentropic core exists only for short distance from the 
leading edge. 

The presence of a core at increased velocity (Figure 3) from the leading edge 
results in higher friction drag (Figure 10) than in the case of a single flat plate. 
This increase is quite substantial for gap Reynolds number below 104. At Re, 


above 10° approximation with single flat plate condition would be quite ade- 
quate for most purposes. 


ai ee 


ee eee ea 
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The effect of acceleration of the core also Shows up significantly in temper- 
ature distribution along insulated plates due to lowering of the temperature of 
the core (Figure 6). The temperature of the wall is highest at the leading edge 
(identical to the single plate case) but decreases toward downstream. This 
effect is more significant at higher inlet stream Mach numbers. In the case 


a 
5 ae 


Figure 12 


Nusselt number of inlet length L for Pr = 0-72, y = 1:40,C = 1 
(Rez = L uy Polity, Nur = hy L|Ryy)- 


of finite plate thickness this effect would be more significant. This means lower 
leading edge temperature than the present case when conductivity of the plate 
is considered. 

For the case of steady wall temperature, when the fluid is heated by the 
wall, the effect of increased boundary layer thickness due to heating, is com- 
pensated by the presence of a core of lowered temperature, resulting in higher 
Nusselt number than the single plate case. When the fluid is cooled by the wall, 
the increased heat transfer rate due to reduced boundary layer thickness as in 
the single plate case, is rendered less favorable due to the core. However, the 
latter modification is not one-sided because of reduced boundary layer thick- 
ness when the gap is finite (Figure 12). 
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Figure 9 also shows that zero heat transfer in the constant wall temperature 
case occurs at 7,, = 7j nearly for subsonic flow only. At higher speeds, the 
direction of heat flow is affected by the signs of different orders of G,,. 

The high Mach number cases considered in this study are pure idealizations, 
since flow discontinuity is neglected here. 

The condition of whether flow separation or convergence to Poiseuille flow 
would occur cannot be obtained from the proposed solution. It can be said, 
however, the flow will converge to Poiseuille motion for Re, below critical [1] 
value. In the present case, however, the critical Reynolds number is affected 
by the rate of boundary layer growth or the temperature distribution. The cri- 
tical Reynolds number will be higher for the case of a cold wall, but lower for 
the insulated wall, and still lower for the case in which the stream is heated 
by the wall. 


Appendix 
Solutions of Equations 


Introducing the stream function y, such that (12) 


_ (oy =, 5 / ey 
Qeu= (se). , Qe = decay » (Al) 


which satisfy the continuity equation, and change the independent variables to 
y and x, the momentum and energy equations [equations (2) and (3)] become: 
Ou i ap 


4 = — — - — + we [news 


Ox Oo dx oy ess 


Oh u ap 1 Oh Oh Ou \2 
= adr —= — es |r 27) = 
Ox Cie dx, eee * oy [¥ aes a ne (se) at (3) 


Further, reducing these equations to the non-dimensional form by introducing 


Pere td nies teat 
Seog Sa 
a y 4 Qo Uy S 
yr = — -—, 3 = 
Qo Ho S Ug C ie Qo Mo S Up C Ae 
u we 
ut = = oh= spa 4 ie een 
U9 Qo’ # Po’ To 
equations (A2) and (A3) are reduced to 
Ou du 0 out 
Cpe fey meee 
Ar etc 5 ps + ut dye ut oor (A5) 
(NI Be» ue duty 1 0 OTs 5 ute +\2 
Ur ———— eg F gt ae ES ( hab) oF 
eg Behar TO * dat Py” dwt dwt Cy Tq \ Owt ) Sa, hel ia 


For the sake of convenience in the solution, substitute into equation (A5) 


yr at 
ree ; i Mi anaes oc £ 
yt yt? (A7) 


= 


ee ee 
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at and pt are simply replaced by x* and y*, The Momentum equation can be solved 
independently first of all. 


Out Out ”) Out 
Pee pe ha as ; 
ae TO a ge ie en 
with the boundary conditions 
ut(0,#*)=1, ut(*, 0) = ut(x*; 1) = 0, see | == 0i (A9) 
Ox* xt = 1/2 


The momentum equation and the boundary conditions can be satisfied by a 
solution of the form given by equation (14) and velocity of the core given by equa- 
tion (15), with & = |/**, and » is defined by [1] 


v= eS Felay eS Fala + EY sal) oe - (A10) 
The boundary conditions as applied to functions /, are (& = 0, 1 = oo): 
#=0,n=9, f,(0) =7,(0) =/,(0) =---=0, 
f(0) = f3(0) = (0) =--- = 0, (A11) 


y= ©, filo) = 2, floo)= 2. By, fg(00) = 2.Ky, fi(00) = 2 Ky---. 


For large 4 


we eA for 2 Ay pt Ay. fe 2 Kg + Ags oh. (A12) 
at 2% 1/2, 4 = 1/4-€, and 
1 1 IS et 
; =§|5¢+4:|+#|53 | A, ! plo ! ,| He (A13) 
Hence 
ad A a 2 A So (A14) 


The functions /,, f, etc. are governed by 
Abhi = 9, 


(A15) 
et+thh-Hh+2Khh=—4 Ky, (A16) 
(A17) 
(A18) 


Phi 2h Ahi = —8 Kes Bi fe 2h ik: Al7 
(+Ah—3hA+4hh = —12K,—12K, BK, +34 — 3h Ahh. 


Equation (A15) of /, is none other than the classical Blasius equation, whose solution 
by piecing up series expansion from 7 = 0 and asymptotic expansion from 7 = co 
is known, and is given in [1]. 
Functions f,, fs, etc. can be determined in similar manners, as depicted in 
Figure 3 as have been presented in the text of this paper. 
The condition that 


A18 


is satisfied, since 


Olpasa °  Oklemxe QoS 9M Inve 28 
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With substitution of equation (A7) 7 and € into the energy equation (A6), 


Oech OL OIL: ua Out uz ee 2 
Sheet Se ee QUE aN eee eel (ERT aa RU = ved ) 
pao + Y On (2 € u*) OE can (2 E ut) uj 0é Gg Le YOR ( 
ale ux Out \2 
ee (D + ile r oe 
Us aaa a ee | 
where 
@ = fy 2 Sie eS Ogee 6 ag ee (A 20) 
The boundary conditions are: 
a= 0 Ti =nOOr, TE ce ly 
- | 
: == A ee 
case i 70 ae T, (421) 
ort | 
: = 0 — i ; 
case lis ; on 0 


The energy equation and the boundary conditions can be satisfied by solution of 
the type given by equations (16) and (17). Since [equations (9) and (10) } 


ary uy, aut dut 

SOL ye 2 Oe pe eS ie 2 pet 

dé e lp. Uy dé (y ) M Uy dé , 

Ue 

M, 1 Cy T == (y— 1) MK, y 
5 2 1 : 

M, = —(K3 + 2 K,) “t= 5 (y -1) M2 (2K, + Ki), fog 

(orn IN 2 


M, = —(y — 1) M*(K, 4+. KK); 


edi 
ee Bee ) ue (2K, + 2K, K,+ K). 


The functions G, have the following boundary conditions 


1 at . 
Go(oo) = M*, G,(o0) =—K,, G,(o0) = —5 (2 Ky + KY), | 
G;(0co) = — (K, + K, K,),--- ; . 
together with for 
Tipe 
caseshy 9G. (0) ee om 
of ) (y nu, 1) Vee G,(0) G,(0) 0, 
case II: G5(0) = G;i(0) = G3(0) =--- = 0. 
For different orders of &, the function G, can be determined from: 
1 “ 1f il 
py tha = = aia) a ; (A23) 


1 
Sa PA , U 1 
Pr A thGi-hG.=—2AGi + BiB. (A24) 


ela ay eh 


Tes - Bile, 


yO es 
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Ik 
pe th — 2HG, = +27, 6/— 37,65 +h; | 
5 a : 1 » (A25) 
+CR FR A+ KA“ LUE+ 2K), | 
1 “ if 4 of. , eS = 2 
p, % +7, G,—3 fi Gs; = —2 fy G3 — 3 fs G3 — 4 f,Gy — 2h Gy + fg Gy + Ky fs | 


! > eaN ts. 4 a Ee 1 
+ 2K FHV R+3(K+ KK) A-SUKT A. | 


Equations (A15) to (A18) and (A23) to (A26) were solved numerically on an 
IBM 704 digital computer. The intermediate steps as given in [3] or [13] and [14] 
were shown to be valid but unnecessary. ‘ ate 

The physical velocity and temperature profiles were obtained by 


*} 


x ’ 4 
pf] Hl mare 


Since f is a function of x alone 
oe PHO) 


x me T+ = T+ > 
leading to equation (13). 
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Zusammenfassung 


Es wird eine Analyse vorgelegt, die erlauben soll, die Bestimmung der Reibungs- 
und Warmeiibergangs-Kenngréssen einer kompressiblen, laminaren Strémung in 
einem unendlich grossen Gitter flacher Platten vorzunehmen. 

Die Ergebnisse weichen fiir niedrige Reynoldssche Zahlen bedeutend von jenen 
einer einzelnen flachen Platte ab. Die Rechnungen fiihren zu vergréssertem Rei- 
bungswiderstand, verminderter Riickbildung und fiir den Fall, dass die Fliissigkeit 
die Wand erwarmt bzw. abkiihlt, zu vergrésserter bzw. verkleinerter Nusselt-Zahl. 


(Received: October 20, 1958.) 
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Zur Bestimmung von Ultraspuren von Silberjodid in Regenwasser 


Von Otro GUBELI, Zurich!) 


Silberjodid (AgI) bildet ein sehr schwer lésliches Salz von gelblicher Farbe. Das 
klassische Léslichkeitsprodukt bei 25°C betragt 1,5-10-16 Mol?/l2. Durch Belich- 
tung farbt es sich allmahlich grauschwarz. 

Die durch Fallung aus wasseriger Lésung bei Zimmertemperatur mit einem 
Uberschluss an Agt-Ionen erhaltene Modifikation bezeichnet man als «-Agl. Sie 
ist bei Zimmertemperatur bestandig, ist kubisch, Zinkblendegitter mit ay = 6,48 A. 

_ Bei 137°C wandelt sich die Zimmertemperatur-Modifikation in fB-Agl um. Sie 
kristallisiert hexagonal, Wurtzitgitter a = 4,59 A, c= 7,52 A. Man erhalt diese 
Modifikation nicht nur beim Erwarmen, sondern auch bei der Fallung von Agt- 
fonen mit einem Uberschuss an I~ bei Zimmertemperatur. Sie ist jedoch nicht 
stabil und zerfallt schon beim Zerreiben des Niederschlages in «-AgI. 


1) Anorganisch-chemisches Institut der ETH. 


a a 
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Oberhalb 146°C ist die Modifikation von y-Agl bestandig. Sie weist ein kubi- 
sches Gitter mit statistischer Verteilung der Ag+-Ionen auf. Die I--Ionen besétzen 
die Ecken und das Zentrum der Elementarzelle ay = 5,03 A. Die Silberionen wur- 
den vollkommen ungeordnet auf 42 Gitterliicken verteilt gefunden. Der Schmelz- 
punkt von Agl liegt bei 552°C. 

Die Komplexbildung der Silberionen ist sehr ausgepragt. In einer wasserigen 
Suspension von Silberchlorid sind noch merkliche Mengen an Ag+-Ionen vorhanden, 
die beim Bromid und schliesslich beim Jodid betrachtlich abnehmen. Die Werte 
der klassischen Léslichkeitsprodukte sind bei 25°C fiir AgCl 1,56-10-1° und fiir 
AgBr 7,7-10- Mol?/l?, Fiigt man zum Silberchloridniederschlag Ammoniak, so 
entsteht das Gleichgewicht mit den vorhandenen Ionen 


Agt + 2 NH; ad (Ag(NHs)2]*. 


Sind geniigend Silberionen vorhanden, so reagieren diese zum bedeutend stabileren 
Ammoniakkomplex. Silberbromid geht bereits schlechter in Lésung und das Jodid, 
weil zu wenig Agt-Ionen vorhanden sind, nicht mehr. Lasst man jedoch auf das 
Agl einen starkeren Komplexbildner, zum Beispiel Cyanid, einwirken, so gelingt 
die Auflésung des Niederschlages erneut. Die Komplexkonstanten dieser Gebilde 
sind, fiir Ammoniak 
AgCl + 2 NH, —> [Ag(NHs3).]* + Cl-, 
[Ag*] [NH3]}? 


K = —>_—_ = 6,3-10-8, 
[Ag(NHs3)F] 


und fiir Cyanid 
AgCl + 2 CN- —> [Ag(CN),.]- + CL, 


[Ag*] [CN=]? 
[Ag(CN)3 ] 


= 9/591. 0—22% 


Eine ganz erhebliche Komplexbildung tritt bereits bei folgendem Salzpaar in Er- 
scheinung: 
Agi + I- —> [Agl,], 


wobei nicht nur die erste Stufe, sondern auch héhere Komplexe, wie [AgI,]~?, autf- 
treten. Beim Verdiinnen mit Wasser zersetzen sich gewohnlich diese Komplexe, 
so dass AgI ausfallt. 

Die Silberionen zeichnen sich dadurch aus, dass ihre Reduktion zum Metall be- 
sonders leicht méglich ist. Die Reaktionswiderstande sind sehr gering, wenn das 
erforderliche Reduktionspotential erreicht ist. Das Normalpotential betragt 0,80 Vv 
bei 25°C. Man fiihrt dies teilweise auf die geringfiigige Hydratation des Silberions 
zuriick. Auch die elektrochemische Potentialeinstellung ist sehr gut. 

Durch Schmelzen der Silberhalogenide mit einem Uberschuss an Natriumkarbo- 
nat erhalt man das Metall. 


4 AgCl + 2 Na,CO, —> 4 NaCl + 2CO, + O, + 4Ag. 


Auf nassem Wege kann man schon in der Kalte, rascher beim Erwarmen mit Eisen- 
2-sulfat reduzieren. Ferner wird aus gewohnlichen und komplexen Silberlésungen 
in saurem Medium mit metallischem Zink das Silber abgeschieden: 


Rhos Chit Foie yp 2 Cla 4 Zak? 32 ABs 
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Die Sammlung des Regenwassers hatte auf diese Verhaltnisse Riicksicht zu 
nehmen. Es wurde an geeigneten Sammelstellen in Proben von 50 1 in Glasballons 
gefiillt und im Laboratorium mit Hilfe eines Diinnschichten-Schnellverdampfers 
aus Glas eingedampft. 

Ein Unterdruck einer Wasserstrahlpumpe geniigt, um mit Dampfheizung 50 1 
Wasser in 8 Stunden auf ein kleines Volumen einzuengen, das mit der Wasch- 
fliissigkeit 500 ml ausmacht. Diese Menge wurde in offenen Porzellangefassen und 
schliesslich im Porzellantiegel offen zur Trockene verdampft. Durch vorsichtiges 
Schmelzen mit Na,CO, und gleichzeitige Verbrennung organischer Substanzen 
wurde reduziert. Nach Extraktion mit Wasser und Filtration konnten gefalltes, 
metallisches Silber und lésliches Jodid getrennt werden. 

Ein quantitativer Silbernachweis von hoher Empfindlichkeit und guter Spezifitat 
ist in der Reaktion mit p-Dimethylaminobenzilidenrhodanin bekannt. In der 


Formel 
NH==CO 


sc. c=c-¢ _\_NICH), 
eae ee (CHa) 


ersetzt das Silberion den Wasserstoff der zyklischen Imingruppe unter Bildung 
sehr schwer léslicher, gefarbter Silbersalze. Das Reagens ist nur in saurem Medium 
brauchbar. Gleichzeitig reagieren noch Quecksilber-, Gold-, Platin-, Palladium- und 
Kupfer-1-ionen. 

Zur Ausfiihrung der Mikrobestimmung wird Filterpapier mit einer gesattigten 
Lésung des Reagens in Azeton getrankt. Das Papier lasst man trocknen und bringt 
einen Tropfen der schwach salpetersauren, gepufferten Silberlésung darauf. Mit 
sehr niedrigen Silbergehalten entstehen auf dem gelbbraun gefarbten Papier violett- 
rote Flecken. Die Erfassungsgrenze liegt bei 0,02 wg Silber in 0,2-n Salpetersaure ; 
die Grenzkonzentration betragt 1:2500000, Verwendung von 0,05 ml Lésung. 

Aus groésseren Fliissigkeitsmengen wurde der Nachweis so ausgefiihrt, dass das 
Reagens zur sauren Silberlésung gegeben wurde. Nun tiberschichtete man mit Ather 
und schiittelte durch. Das iiberschiissige Reagens ging in die iiberstehende Ather- 
phase. An der Grenzschicht Ather—Wasser bildete sich eine violette Haut. Die 
Erfassungsgrenze betrug 1 ug Silber in einer Grenzkonzentration von 1:5000000, 
wobei 5 ml Lésung verwendet wurden. 

Zum Nachweis neben stérenden Fremdionen sind folgende Méglichkeiten ge- 
geben: 

Neben Quecksilber und Edelmetallen wird die Probelésung mit 5prozentigem 
Kaliumcyanid versetzt und nun mit 2-n Salpetersdure schwach angesduert. Die 
Cyanide des Quecksilbers und der Edelmetalle dissoziieren so wenig, dass nur noch 
das Silberion mit dem Benzilidenrhodanin reagiert. Die Empfindlichkeit hat ein- 
gebiisst, so dass die Erfassungsgrenze bei 0,5 ug Silber mit der Grenzkonzentration 
1: 100000 liegt. 

Neben Quecksilber- und Kupfer-1-Salzen fiihrt die Maskierung mit Cyanid nicht 
zum Ziel. Man kann mit Salzsiure und Ammonchlorid versetzen, wobei die Chloride 
gebildet werden und das Silber dank der Stabilitat des Benzilidenrhodanins trotz- 
dem reagiert. Die Empfindlichkeit hat nochmals abgenommen. Die Erfassungs- 
grenze hegt bei 0,65 ug Silber mit einer Grenzkonzentration von 1:80000, 

Grosseren Schwierigkeiten begegnet der quantitative Spurennachweis von 
Jodid im aufgearbeiteten sodaalkalischen Filtrat. Zundchst wird das Filtrat ein- 
geengt, Jodid nach schwachem Ansiuern mit Salzsdure unter Zugabe von etwa 5 g 


Vol. X, 1959 Kurze Mitteilungen — Brief Reports — Communications bréves 313 


Eisen-3-chlorid oxydiert und in einem Destillierkolhen bei 90°C im Stickstoffstrom 
in eine Vorlage mit etwa 5 ml Tetrachlorkohlenstoff abdestilliert. Jodmengen in der 
Grossenordnung von 10-4 g farben den Tetrachlorkohlenstoff deutlich violett. Fiir 
Spuren in der Gegend von yg muss der Tetrachlorkohlenstoff mit 10 ml Wasser und 
etwa 0,5 g elektrolytisch dargestelltem, feinverteiltem Cadmium geschiittelt wer- 
den. Dabei geht das zu Jodid reduzierte Jod in die wasserige Phase iiber. Aus der 
wasserigen Lésung wird erneut nach dem Ansauern mit Salzsiure mit Eisen-3- 
chlorid oxydiert. Als Reagens verwendet man mit Fluoreszeinnatrium getrankte 


Zirkulation i Sone 
ee 


Regulier- 
hahn 


Kuhlwasser- 
ableitungen 


Kihlwasser= 
Zuleitungen 


(_ Wasserstrah/- 


Dampf-_- 

neizung SAUGPUMpE 
Kalben mit 

AblaBhahn fur J Destillat +10 | 

eingedampfte 

Wasserprobe Ansaugleitung 


fur Wasserprobe 


Dunnschicht-Schnellverdampfer aus Glas fiir Regenwasser. 


Papierstreifen. Beim Erwarmen der jodidhaltigen Probe entsteht auf dem dariiber- 
gelegten gelbgriin gefarbten Papierstreifen je nach Menge ein oranger bis roter 
Fleck. Durch Eichung mit bekannten Gehalten lasst sich die Menge halbquantitativ 
ermitteln. 

Das geschilderte Verfahren fiir den AgI-Spurennachweis ist im Auftrage der 
Eidgenéssischen Kommission zum Studium der Hagelbildung und der Hagelab- 
wehr entwickelt worden. Es fiihrte an zwei ausgewahlten Tagen des Versuchsjahres 
1954 bzw. an einem Tag des Versuchsjahres 1956, an denen im Grossversuch II zur 
Bekampfung des Hagels in der Magadinoebene die Atmosphare vermittelst AgI- 
Bodengeneratoren geimpft wurde, zum Nachweis des Impfstoffes im Regenwasser, 
wobei die Niederschlagsproben von je 50 1 jeweils an einer (Gordola) bzw. drei 
Stellen (Gordola, San Antonino und Taverne) des Versuchsgebietes gesammelt 
_ worden sind2), Bei Probenentnahmen an ungeimpften Tagen (total 2 Blindproben- 
-analysen) konnten keinerlei Spuren von Agl festgestellt werden. 


2) Tatigkeitsbericht Nr. 7, Versuchsjahr 1954, bzw. Nr. 9, Versuchsjahr 1956 des Grossver- 
suches II zur Bekampfung des Hagels in der Magadinoebene, herausgegeben von der Abteilung fiir 
Landwirtschaft des EVD; S. 38 bzw. 65. 
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Summary 


A process is given for determining ultra-traces of silver iodide in rainwater 
following on seeding experiments designed to influence precipitation. After con- 


sidering the physico-chemical properties of silver iodide a 50-1 sample of rain- — 


water has been vaporized in a rapid ‘thin-film’ vaporizer. The residue is reduced 
by melting with soda. The reduced silver can be separated from the iodide by 
extraction with water. The quantitative analysis is then made by spot test on 
the basis of drops placed on filter paper saturated in p-dimethylaminobenzilidene- 
rhodanine, compared with drops of known content or alternatively with non- 
seeded rainwater. The quantitative analysis of the iodide is made with the 
aqueous, soda-alkaline, silver-free extract. The iodide traces are enriched, oxy- 
dized, and made to react on strips of paper soaked in fluorescein-sodium. 


(Eingegangen: 8. Dezember 1958.) 


Deflation bei Bandmatrizen 


Von HEINz RUTISHAUSER, Ziirich?) 


Sei 4, ein Eigenwert, v, ein zugehoriger, normierter Eigenvektor der symme- 
trischen Matrix A. Dann ist bekanntlich : 


Bs A el, odo. (1) 


eine Matrix mit folgenden Eigenschaften: 

a) B besitzt nm unabhangige Eigenvektoren, die ebenfalls Eigenvektoren von 4 sind. 
b) Sind ,, Ag, 43, ..., An die Eigenwerte von A, so sind 0, Ay, As, ..., A, diejenigen 
von B. 

Man bezweckt mit diesem von HorELLinc eingefiihrten Verfahren, den «Beitrag» 
des absolut grossten Eigenwertes in der Matrix A zu eliminieren. Dadurch wird es. 
moglich, auch mit Methoden, die an sich nur den absolut gréssten Eigenwert einer 
Matrix liefern, dennoch samtliche Eigenwerte und Eigenvektoren zu berechnen. 

Die Vor- und Nachteile dieses Vorgehens sind hinreichend bekannt, so dass sich 
eine weitere Erorterung eriibrigt. Wenn man jedoch dieses Deflationsverfahren auf 
Bandmatrizen, das heisst Matrizen (a;,) mit der Eigenschaft 


a;,,=0 fir |i—k| >m (2) 


anwendet, so zeigt sich ein weiterer Nachteil: Die Matrix B ist dann namlich aus- 
gefullt, das heisst, die fiir die numerische Berechnung vorteilhafte Eigenschaft (2) 
geht durch die Deflation nach (1) verloren. 

; nes wird daher in dieser Arbeit ein von (1) etwas abweichendes Deflationsver- 
ahren entwickelt, welches die Bandgestalt einer symmetrischen Matrix nicht zer- 


stort, aber auch auf beliebige symmetrisch i 1 i 
BO ieee ge sy e Matrizen mit Vorteil angewendet 


1) Institut fiir angewandte Mathematik der ETH. 
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1. Orthogonale Hessenber¢ -Matrizen 


Ist wv ein zunachst beliebiger Vektor mit den Komponenten +1, %5, ) aso 
kann man sofort n — 1 weitere Vektoren angeben, die mit y zusammen ein Ortho- 
gonalsystem bilden. Wir setzen zundchst: 


Un = (41, %, ---, Xn_p ¥n) = Begebener Vektor, 
ferner: 
On-1 = (¥1, Noy +++ ¥y_y, Yn) ; 
wobei die letzte Komponente y, so bestimmt wird, dass v,_, und vy, orthogonal 
sind; aber im Falle ¥, = 0 setzen wir 
Un-1 = (0, 0, eee, (uy 1) . 
Allgemein definieren wir rekursiv 
Un = (¥%y) %q, «+s Her Vexp 0, 0,..., 0), (3) 
wobei ¥;,,, so gewahlt wird, dass yv, zu vy, orthogonal wird. Nur im Falle, dass 
¥p41 = O ist, setzen wir 
v;, = (0, 0,..., 0,1, 0,...,0) [= (Rk + 1)-ter Koordinatenvektor]. 
Damit wird v, offenbar von selbst auch orthogonal zu v,44,..-, Un_1- 
Dieses Verfahren bricht wie folgt ab: Falls die p ersten Komponenten (p = 0,1, 
..., 2#—1) von v verschwinden, ist der letzte nach (3) konstruierbare Vektor offenbar 
Bee (& scan Ohi Vaan oars, O).. 


Fiir die tibrigen  Vektoren des Orthogonalsystems wahlt man dann einfach die p 
ersten Koordinatenvektoren (das heisst die ersten  Kolonnenvektoren der Ein- 
heitsmatrix). 

Die so konstruierten Vektoren vj, ..., v, bilden ein Orthogonalsystem, welche 
nach Normierung die Kolonnenvektoren einer orthogonalen Matrix H bilden, die 
offenbar von sogenannter Hessenberg-Gestalt ist, das heisst 


| * # * * * OK 
|e oe * * || 
Qo * * oe 
00 *-. * ok 
ie] : | (4) 
0) 020 Omens 


Die letzte Kolonne dieser Matrix enthalt bis auf einen konstanten Faktor den 


_ gegebenen Vektor vy, welcher die Matrix H im wesentlichen eindeutig bestimmt. 


Es ist im iibrigen bemerkenswert, wie gering der Rechenaufwand fir die Konstruk- 


tion dieser Matrix H ist. Beispiel: » = 5, v = (1, —1, 0, 1, —1). Dann wird 


eal SS aia 


|| 0,707107. 0 0,408248 0,288675 ~—0,,5 

0,707107. 0 —0,408248 -0,288675 -0,5 
=i. 1 0 0 0 

0 0 -« -0;816496. —.0/288675 .. 0/5 

0 0 0 0,866025 -0,5 
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Man stellt leicht fest, dass die so konstruiette Matrix H auch in der Form 


|| cosy, sing, cosp, sing, SiN Py COS Ps 
| —sing, COS@, COSY, COSY, SIN P, COS Hs 
| 0 —sin gy, COS Pz COS Ps 
| 0 0 —sin Pei 
we ie ; (4a) 
| : 
| 0 0 0 


das heisst als Produkt von  — 1 elementaren « Jacobi-Drehungen» 
1 


| 2 O 


cosp; sing, Zeile 7 
—sing; COSg; Zeile 7 + 1 


0 


1 


: 


dargestellt werden kann. Dies ist zwar fiir unsere weiteren Betrachtungen uner- 
heblich, hat aber eine grosse Bedeutung fiir die numerische Rechnung. 


2. Deflation durch Transformation mit H 


Sei A ein Eigenwert der Matrix A und v ein zugehériger Eigenvektor. Wenn wir 
zu diesem Eigenvektor uv die zugehdrige orthogonale Hessenberg-Matrix H kon- 
struieren, so hat die Transformierte 


14 =HATAHA (5) @ 
offenbar die Gestalt é 


1A = (6) 


wobei wil wieder eine symmetrische Matrix ist, welche die iibrigen Eigenwerte der 
Matrix A enthalt. An dieser kann daher die Eigenwertbestimmung fortgesetzt 
werden, so dass man schliesslich simtliche Eigenwerte von A erhalt. 


Der bedeutendste Vorteil der Deflation nach (5), (6) gegeniiber dem Verfahren 
(1) besteht in 


; 
4 
_Theorem 1: Sei A eine symmetrische Matrix mit der Eigenschaft (2), v em | 
Eigenvektor derselben und H die aus v evzeugte orthogonale Hessenberg-Matrix. : 
Alsdann hat die aus A durch Deflation nach (5) gewonnene Matrix 14 ebenfalls die : 


cil ah Mie ndie Bh 


vse 


ie 3 ia! hE lh ah Naha a 
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Exgenschaft (2), mit anderen Worten: die Bandgéstalt bleibt durch diese Art dev De- 


flation erhalten. 


Beweis: Wir zeigen zunachst: 
a) Die Matrix B = (A — AE) H hat die Eigenschaft 


a . {k>p, wenni<k+1—m 
Be a ee wenni<k—m NG (7) 
falls vy ein Eigenvektor von A und 4 der zugehérige Eigenwert ist. 

In der Tat ist ja fiir k > p die k-te Kolonne v, von H eine Linearkombination 
der Vektoren v,, e; 4a) +++) @, (wobei e; den i-ten Kolonnenvektor der Einhéits- 
matrix bezeichnet), denn die ersten k Komponenten von y,, stimmen bis auf einen 
gemeinsamen Faktor mit denen von yv, iiberein. Damit ist die k-te Kolonne von 
B= (A —AE) H, das heisst der Vektor (A — AE) v,, offenbar eine Linearkombina- 
tion der Vektoren 


(& — AE ee = Orla — HE) efits. (4 = 2B) ex, 


also eine Linearkombination der (& + 1)-ten bis n-ten Kolonne von 4 — 4 E, woraus 
sich die erste Behauptung in (7) ergibt. 

Fiir k S p ist die k-te Kolonne von B hingegen gleich derjenigen von A — AE, 
womit sich auch die zweite Behauptung in (7) ergibt. 

b) Es wird nun noch das Produkt 14 = 4E+ HT B gebildet, wobei B die Eigen- 
schaft (7), HT dagegen die Eigenschaft 


ee SS re a 
nf, = 0 tar |= o Lae (8) 


besitzt. Damit folgt aber fiir 14 sofort 1a,, = 0 fiirk —7 > m und damit wegen 
der Symmetrie von 14 auch die Eigenschaft (2), was zu beweisen war. 
Numerisches Beispiel: 


SCOrRrFW 
OorFRFR AH 
PrP pare 
Peper OO 
Sat cS) 1) 


nm 
| 


Hier ist vy = (1, —1, 0, 1, —1) ein Eigenvektor zuin Eigenwert 2. Somit ergibt 
sich die Matrix H, die in Abschnitt 1 numerisch angegeben wurde, und damit 


| 4,5 1,414214  —0,866025 0 0 | 
1,414214 A. —0,816496 1,154701 0 
14 = || —0,866025 —-—0,816496 4,166667 —0,942809 0 
0 1,154701  —0,942809 373335339 0 
ae, 0 0 0 2 


Natiirlich sind die Elemente einer so berechneten Matrix 1A mit Rundungs- 


- fehlern behaftet, doch lassen sich diese wegen der Orthogonalitat von H relativ 


leicht abschatzen. 


2) Dabei ist » wie in Abschnitt 1 dadurch definiert, dass ¥,,, die erste nicht-verschwindende 
Komponente des Vektors v ist. 
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3. Pseudodeflation durch Transformation mit H 


Der grosse Vorteil jeder «Deflation durch Ahnlichkeitstransformation» besteht 
darin, dass die Verwendung eines ungenauen Eigenvektors wu lediglich zur Folge 
hat, dass die « Rander» der Matrix 14 nicht mit Nullen, sondern mit kleinen Grossen 
besetzt sind. Da die Eigenwerte von !4 dadurch nicht beeinflusst werden, hat dies 
lediglich zur Folge, dass man nicht mit der (n — 1)-reihigen Matrix ,A, sondern mit 
der vollen Matrix 14 weiterrechnen muss*). 

Die weiteren Deflationsschritte zur Bestimmung samtlicher Eigenwerte ver- 
laufen daher nicht gleich wie bei der exakten Deflation: Zwar bestimmt man zu- 
nachst ebenfalls einen Eigenwert ,A und Eigenvektor ,v der (w — 1)-reihigen Matrix 
A, und ebenso auch die zu yu gehdrige (m — 1)-reihige orthogonale Hessenberg- 
Matrix ,H. Die letztere wird aber sogleich wie folgt zu einer »-reihigen Orthogonal- 
matrix 1H erganzt: 


0 


Wee = 


Wenn man nun die Transformation 
BA =) (7H) ie (*A agi) 


durchfiihrt, so hat die wieder symmetrische Matrix 24 die Gestalt 


Bie / 
7A 
yi = ; 
Bg © 
€ é|,A ¢€ 
Eo yee ee ak 


wobei die ¢ kleine Gréssen andeuten, die durch genauere Bestimmung der Vektoren 
,v und v beliebig klein (praktisch jedoch nicht kleiner als Rundungsfehler) gemacht 
werden konnen. 

_In dieser Weise fortfahrend, erhalt man schliesslich, ungeachtet der Ungenauig- 
keit der verwendeten Eigenvektoren y, ,v, ,v usw., eine zu A dhnliche Matrix "A 
deren Aussendiagonalelemente klein sind, deren exakte Eigenwerte daher relativ 
leicht berechnet werden kénnen. 

_ Die Bandgestalt einer Matrix wird allerdings durch die Verwendung ungenauer 
Eigenvektoren zur Konstruktion der Matrizen H fiir die Deflation nach (5) im all- 
gemeinen zerstort. Die Eigenschaft (2) bleibt jedoch erhalten (bis auf Reste, die 
durch Rundungsfehler erzeugt werden), wenn man nur dafiir sorgt, dass der ore 


: * 4 
chai aoe dazu fiihrt die Verwendung eines ungenauen Eigenvektors v in (1) zu irre- 
parablen Fehlern, die sich bei weiteren Deflationsschritten in gefahrlicher Weise anhaufen kénnen. 
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wendete Eigenvektor v die ersten » — m Gleichungen des Systems (A —AE) v=0 
exakt erfiillt. 


Summary 


Deflation methods which are used to eliminate dominant eigenvalues of a matrix 

in order to compute also non-dominant ones, in general destroy special forms of a 
matrix, such as the property 

@;,=9 for |t—k| >m. (1) 


(In some cases even the symmetry of the matrix is destroyed by deflation.) In the 
present paper, the author describes a deflation method which maintains the sym- 
metry as well as the property (1). 


(Eingegangen: 16. Dezember 1958.) 
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10. Sitzung des Schweizer Komitees fiir Optik vom 
18. November 1958 in Ziirich 


Am 18. November 1958 fand im Physikgebaude der ETH in Zirich die 10. Sit- 
zung des Schweizer Komitees fiir Optik statt. Anlasslich dieser Sitzung wurden 
5 Vortrage gehalten, wovon zwei Berichte tiber die internationalen Kongresse von 
Briissel (Die Optik in der Metrologie) und Berlin (Elektronenmikroskopie) und drei 
wissenschaftliche Mitteilungen. Wir geben nachfolgend eine Zusammenfassung der 
drei Mitteilungen. 


Uber die Evakuierungsgeschwindigkeit von Bedampfungsanlagen. 
Von Tu. Kraus, Balzers?). 

Die optischen und mechanischen Eigenschaften von Aufdampfschichten sind 
sehr stark von der Desorption der an den zu bedampfenden Objekten und an der 
Rezipientenwand adsorbierten Gase (hauptsachlich Wasserdampf) abhangig. Das 
Ausmass der Desorption gibt sich durch die Evakuierungsgeschwindigkeit des 
Rezipienten zu erkennen. Anhand der Untersuchungsergebnisse konnte gezeigt 
werden2), dass unter gewissen vereinfachten Bedingungen die Pumpzeit ¢ folgender 


: _ einfachen Regel gehorcht: 


| ia oe 


G 
oR De 


(V Volumen des Rezipienten, S Sauggeschwindigkeit der Pumpe, p Druck im Rezi- 


; _pienten zum Zeitpunkt ¢, Pg erreichbarer Enddruck, C eine von den Anfangsbedin- 


gungen abhangige Konstante, die haufig vernachlassigbar klein wird, wenn ¢ vom 


1) Gerdtebauanstalt Balzers. 
2) Vel. Th. Kraus, Uber die Evakuierungsgeschwindigkeit von Hochvakuumanlagen, Vakuum- 


~ technik 8, 39-43 (1959). 
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Einsetzen der Diffusionspumpe an gerechnet wird). Die Grésse P4, welche die 
Dimension eines Druckes besitzt, charakterisiert den Zustand des Rezipienten in 
bezug auf die Adsorption. Ihr Betrag liegt bei Bedampfungsanlagen bei Zimmer- 
temperatur je nach «Verschmutzungsgrad» in der Gréssenordnung von 107° bis 
10-1 Torr und kann innerhalb des Giiltigkeitsbereiches der Regel (Py > p > Pr) 
jeweils aus dem ¢ — [1/(p — Pg) ]-Diagramm ermittelt werden. Damit kénnen Be- 
dampfungsanlagen auf ihre Einsatzfahigkeit gepriift und durch rechtzeitige Vor- 
kehrungen (zum Beispiel Reinigung) Ausfalle vermieden werden. 


Ein Objektiv fiir Profilprojektion. Von W. Lormar, Aarau%). 

Projektionsobjektive fiir Messzwecke miissen bei vorgegebenem AbbildungsmaB- 
stab ein verzeichnungsfreies Bild im telezentrischen Strahlengang liefern (Figur 1). 
Auf industriellen Profilprojektoren, die zur Stiickkontrolle auf Masshaltigkeit 
dienen, sind in der Regel mehrere auswechselbare Objektive mit verschiedenen 
normierten AbbildungsmaBstaben vorgesehen. In diesem Fall kommt zu den obigen 
Forderungen noch diejenige nach Abstimmbarkeit der Objektivbrennweiten hinzu, 
da aus konstruktiven Griinden die Distanz Objekt—Bild beim Vergrésserungs- 
wechsel erhalten bleiben soll. 


Bild 


Figur 1 


Telezentrischer Strahlengang. 


Objekt 


Bild 
Figur 2 
Strahlengang bei episkopischer Beleuchtung. 


Bei Profilprojektoren wird ausserdem oft gewiinscht, opake Objekte mit der- 
selben Optik in Epi-Projektion abbilden zu kénnen. Der Strahlengang ist hierbei 
nicht mehr telezentrisch, sondern das Objektiv wird von jedem Objektpunkt aus 


mehr oder weniger mit voller Apertur beansprucht (Figur 2) 


3) Kern & Co. AG. 
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Wegen seiner guten Korrektureigenschaften ist bisher fiir derartige Zwecke viel- 
fach der sogenannte Gauss-Typ verwendet worden (Figur 3). Die Abstimmung 
erfolgt durch Veranderung des mittleren Luftabstandes, wodurch jedoch unver- 
meidbarerweise die Bildebnung in ziemlich empfindlichem Masse beeinflusst wird. 

Die Suche nach einem korrekturmiassig gleichwertigen Typus mit geringerer 
Empfindlichkeit auf Brennweitenabstimmung fiihrte auf den in Figur 4 gezeigten 
Aufbau, der aus dem bekannten, in Figur 5 dargestellten Typ durch Einfiihren 
zweier sammelnder Kittflachen: in den Negativgliedern hervorgeht. Der vor- 
geschlagene Typ ist zwar fiir photographische Zwecke bereits bekannt geworden 4), 


(ea 


Figur 3 Figur 4 Figur 5 
Gauss-Typ. Aus dem Grundtyp von Figur 5 
entwickelter verbesserter Typ. 


war aber wenig leistungsfahig. Durch Verwendung moderner Glaser (Lanthan- 
Kron) lasst sich eine fiir Profilprojektion geeignete Korrektion erreichen, bei sehr 
geringer Abhangigkeit derselben vom Mittelabstand sowie guter Bildscharfe auch 
fiir Epi-Projektion. Fiir einen Bildwinkel von rund 16° und ein Offnungsverhaltnis 
von 1:2,7 betragt der Radius der Bildschale bei einem Ausfiihrungsbeispiel das 
2,9fache der Brennweite, was fiir den vorgesehenen Zweck ausreicht. 


Zur linearen Fehlertheorie optischer Prismen, insbesondere des 
Tripelprismas. Von R. STETTLER, Aarau). 

Die Optik unterscheidet bekanntlich zwischen Dispersions- und Spiegelprismen. 

Wahrend Dispersionsprismen zur Erzeugung von Spektren dienen, verwendet 
man Spiegelprismen, um Strahlengange abzulenken und die Orientierung optischer 
Bilder zu verandern. 


1. Allgemeines iiber Spiegelprismen 


Jedes ideale Spiegelprisma ist optisch aquivalent einem System von Plan- 
spiegeln, deren Lage mit derjenigen der reflektierenden Prismenflachen iiberein- 
stimmt, und einer Planplatte, deren Dicke gleich dem sogenannten Glasweg ist. 

Die gegenseitige Lage der Prismenflachen ist an gewisse Winkelbeziehungen 
gekniipft. Sind diese nicht genau erfiillt, so entstehen kleine Ablenkungsfehler, 
schwache Farbenzerstreuungen und eventuell Bildverdopplungen. Wahrend die 
beiden erstgenannten Erscheinungen meist unerheblich sind, sind Bildverdopp- 
lungen, wie man sie zum Beispiel beim Dachkantprisma beobachtet, im allgemeinen 
unzulassig. Die Bildverdopplung beruht darauf, dass Spiegelungen — als mathe- 
matische Transformationen aufgefasst — im allgemeinen nicht vertauschbar sind. 
Wie man leicht nachweist, sind Spiegelungen an zwei Reflexionsflachen dann und nur 
dann vertauschbar, wenn diese einen rechten Winkel bilden. Jedoch zeigen nicht 
alle fehlerhaften Prismen Bildverdopplungen. Solche entstehen nur, wenn das ein- 
fallende Biindel durch eine Prismenkante in zwei Teilbiindel aufgeteilt wird (wie 


* 4) US-Patent Nr. 1565205 (P. Rudolph). 
°) Kern & Co, AG. 
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zum Beispiel beim Dachkantprisma). Ist namlich der Prismenwinkel dieser Kante 
nicht exakt 90°, so besitzen die beiden reflektierten Teilbiindel einen Richtungs- 
unterschied, was dann zur Entstehung eines Doppelbildes fihrt. 

Fiir die technische Optik ergibt sich damit die Aufgabe, die zulassigen Tole- 
ranzen der Winkelfehler zu bestimmen. Das Verfahren ist einfach. Gegeben sei ein 
ideales Spiegelprisma. Ein einfallendes Parallelbiindel, festgelegt durch den Rich- 
tungsvektor s, (|| = 1), wird durch das Prisma verfolgt und trete mit dem Rich- 
tungsvektor sj (|sj| = 1) aus. Dann werden die Prismenwinkel mit den Winkel- 
fehlern «, (|«;| <1) behaftet, worauf das Parallelbiindel s) erneut durch das Pris- 


VASES LS 604 
Figur 1 


Zum Reflexionsgesetz. 


mensystem verfolgt wird. Fiir die Fehlerbetrachtung darf man dabei jedes nicht- 
ideale Spiegelprisma durch das entsprechende Spiegelsystem ersetzen. Wie man 
zeigen kann, sind die Ablenkungsfehler eines nichtidealen Spiegelprismas gleich den 
Ablenkungsfehlern des entsprechenden Spiegelsystems, multipliziert mit der 
Brechungszahl des Prismenglases (vorausgesetzt, dass die Inzidenzwinkel an der 
Ein- und Austrittsflache des Prismas nicht zu gross sind). Bei der Strahlendurch- 
rechnung ist noch zu beriicksichtigen, dass fiir Teilbiindel die Reihenfolge der Re- 
flexionsflachen eventuell verschieden sein kann; dadurch wird das urspriinglich 
einheitliche Biindel s, in k = 1 austretende Parallelbiindel aufgespalten. sj, sj, ..., sj; 
seien die Richtungsvektoren dieser Biindel. Dann liefern die Differenzvektoren 
A; = s;, — so die gesuchten Ablenkungsfehler, die natiirlich Funktionen der «, sind. 

Die Strahlendurchrechnung erfolgt vorteilhaft mit Vektorrechnung unter Zu- 
grundelegung des vektoriellen Reflexionsgesetzes (Figur 1): 


Bei Spiegelsystemen mit Winkelfehlern ist eine Strahlendurchrechnung um- 
standlich. Eine Vereinfachung ergibt sich dadurch, dass nur die 1. Potenzen der «, 
beriicksichtigt werden (lineare Fehlertheorie). 


2. Das Tripelprisma 


Die vorstehenden Ausfiihrungen sollen — mit einer bemerkenswerten Abande- 
rung — auf das Tripelprisma angewandt werden. 

Drei rechtwinklig aufeinanderstehende Spiegel bilden einen Tripelspiegel; das 
Tripelprisma entsteht dadurch, dass man aus einem Glaswiirfel eine gleichseitige 
Ecke abschneidet (Figur 2). 

Das ideale Tripelprisma hat die interessante Eigenschaft, dass jeder einfallende 
Strahl, nach Totalreflexion an den 3 Seitenflachen, parallel versetzt zum einfallenden 
Strahl austritt. Somit wird ein einfallendes Parallelbiindel in sich selbst zuriick- 
gespiegelt, wobei die einzelnen Strahlen Parallelversetzungen erleiden. 


poe) 
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Fiir das Tripelprisma soll nun die lineare ehlertheorie entwickelt werden. 
Gemass den obigen Ausfiihrungen kann man sich aut den Tripelspiegel beschranken. 
Beim fehlerhaften Tripelspiegel sind die Kantenwinkel mit den Winkelfehlern «. 
behaftet (Figur 3). 


Figur 2 Figur 3 
Ideales Tripelprisma. Fehlerhaftes Tripelprisma. 


In Abanderung der geschilderten Vektormethode erweist sich im Fall des 
Tripelprismas ein Verfahren mittels Begriffen des Ricci-Kalkiils als wesentlich 
wirkungsvoller. 

Da die Kanten eines fehlerhaften Tripelspiegels offenbar ein «schwach» schief- 
winkliges Koordinatensystem bilden, so hegt der Gedanke nahe, die Strahlendurch- 
rechnung auf dieses System zu beziehen. Man denkt sich deshalb in den drei 
Kantenrichtungen die Einheitsvektoren e,, e,, e, abgetragen (Figur 3). Dann erhalt 
man fiir den metrischen Fundamentaltensor g,, = e; e, in linearer Approximation 


Ge 1 @=1,2-3), Ge —%, bis=—%, 823 = —%, 
und daraus in bekannter Weise fiir den kontravarianten Fundamentaltensor 


Bra Dt 12) oe Ny 8, BO, 

Bedeuten s* und s* die kontravarianten Komponenten des Richtungsvektors des 
einfallenden bzw. reflektierten Strahls, ferner n’ und n, die kontra- bzw. kovarianten 
Komponenten der Spiegelnormalen, so gilt in Ubertragung des Reflexionsgesetzes in 
den Ricci-Kalkiil 


Gas 2uypstn, “ss, = ss nn, =). 


Die Spiegelnormale berechnet sich nach bekannten Regeln. Zum Beispiel gilt fiir die 
Spiegelnormale auf die 3. Flache, aufgespannt durch e, und eg, in linearer Approxi- 
mation 

B= Ot Ms = Oe Ns 1", 
und 
Te) Diy ee TE ty a ole 


Damit sind alle Hilfsmittel fiir die Strahlendurchrechnung bereitgestellt. Besonders 
zu beachten ist noch die Reihenfolge der Reflexionsflachen. 
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Ohne auf Einzelheiten der Strahlendurchrechnung weiter einzugehen, fassen wir 
die auf dem beschriebenen Weg erhaltenen Resultate wie folgt zusammen ®) : 

Man denke sich 2 Tripelprismen, ein ideales und ein fehlerhaftes, deren Deck- 
flachen parallel seien. 

Fallt dann ein Parallelbiindel s, auf das ideale Tvipelprisma, so wird es als ein- 
heitliches Parallelbiindel mit dem Richtungsvektor sj; = —s, zuriickgespiegelt. 

Fallt dagegen das Parallelbiindel s, auf das fehlerhafte Tripelprisma, unter 
kleiner Inzidenz gegen die Deckflache, so wird es im allgemeinen in k = 6 Parallel- 
biindel si, ss, ..., s4 aufgespalten. Werden die Vektoren s; bis s; von einem Punkt 
aus abgetragen, so bilden sie die Kanten einer sechsseitigen geraden Pyramide. Ihre 
Grundflache ist ein Sechseck besonderer Art: Zwei aufeinanderfolgende Seiten 
schliessen einen Winkel von 60° oder 120° ein (geeignete Numerierung vorausgesetzt), 
wobei gegeniiberliegende Seiten gleich lang und parallel sind. 

Die Pyramide ist raéumlich so orientiert, dass ihre Héhe mit s; zusammenfallt 
und ausserdem je 2 gegeniiberliegende Seiten des Sechsecks zu den Seitenflachen des 
Tripelprismas parallel sind. Die Ablenkungsfehler 6; = <x (sj, sj) der 6 Parallel- 
biindel gegeniiber dem ideal reflektierten Biindel s) berechnen sich wie folgt’) : 


ign ee! se eres 
, Sees 1 | 2 2 2 = 
Oi tO, 5 V0 m Vat + oS + a3 + Oy Hy + Oy Ay — My Oy, 
2 — E 
eh any fae 2 2 ! 
0 3 V6 1 Vat + HS + OF + Oy Hy — Hy Hy + Hy Oy, 
: De ee ee = 
Os = (Oy 3 VOM aT + Ob + Of — Oy Oy + Oy My P Oy Oy 


Hierin bedeutet 7 die Brechungszahl des Prismenglases. 
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Handbuch der Physik — Encyclopedia of Physics. Herausgegeben von 
S. Fttcce, Band 50: Astrophysik I: Sternoberfldchen — Doppelsterne (Springer- 
Verlag, Berlin 1958). 458 S., 167 Fig.; DM 98.-. 

Vor 30 Jahren erschien das 10bandige Handbuch der Astrophysik, das damals 
am Wendepunkt von der beobachtenden zur interpretierenden Astrophysik einen 
Markstein gesetzt hat. Seither sind die Bindungen der Astronomie an die Physik 
viel enger geworden, so dass es berechtigt und begriissenswert ist, dass im Rahmen 
des jetzt erscheinenden Handbuches der Physik in den der Astrophysik gewidmeten 
5 Banden 50 bis 54 jenes enzyklopadische Werk der letzten Generation eine neue, 
wenn auch weniger umfangreiche Auflage erfahrt. Neben dem Umfang unter- 
scheidet sich das neue Werk von demjenigen der letzten Generation durch die 
Sprache, indem Deutsch seine damalige Vorrangstellung an Englisch abgetreten 
hat, und in der Gliederung des Stoffes, indem auf die die ersten Bande des Hand- 


8) Wie nachtraglich festgestellt wurde, ist die Feblertheorie des Tripelprismas auf ahnlichem 
Weg von G. RoseNnTHAL, Optik 4, 391 (1948/49), entwickelt worden. Vgl. auch P. R. YopEr jr., 
J. opt. Soc. Amer. 48, 496 (1958). 


*) Fiir Biindel, bei denen die Reihenfolge der Reflexionsflachen umgekehrt ist, ergeben sich 
gleiche Ablenkungsfehler. 
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buches der Astrophysik fiillenden physikalischen Grundlagen im Rahmen eines 
Handbuches der Physik verzichtet werden konntes» So prasentiert sich das neue 
Werk, nach dem Prospekt zu schliessen, als eine Aneinanderreihung von Mono- 
graphien. Eine zu starke Aufgliederung des Stoffes wurde, sehr zum Vorteil des 
Werkes, vermieden, wobei allerdings in Kauf genommen werden musste, dass 
manche Autoren iiber Gebiete zu berichten haben, die weit ausserhalb des ihnen 
durch eigene Forschung vertrauten Stoffes liegen, 

Der vorliegende Band 50, Astrophysik I, Sternoberfldchen und Doppelsterne, be- 
handelt in 10 Beitragen die Spektralklassifikation, Sterne mit speziellen Spektren, 
die Spektren der weissen Zwerge und der planetarischen Nebel, die Molekiilbanden 
in Sternspektren, die visuellen, Bedeckungs- und spektroskopischen Doppelsterne 
und schliesst mit zwei theoretischen Beitragen iiber die Sternatmospharen und die 
planetarischen Nebel. Alle Artikel sind von kompetenten Autoren verfasst. straff 
geschrieben und dem neuesten Stand angepasst. Die technische Ausstattung ist in 
jeder Hinsicht vorziiglich; auch wurde durch das Anbringen der Literaturzitate am 
Fusse der Seiten mehr an die Bequemlichkeit des Lesers als an diejenige des Setzers 
gedacht. 

Die vier noch ausstehenden astrophysikalischen Bande werden den inneren 
Aufbau und die Entwicklung der Sterne, das Sonnensystem, die Sternsysteme (inkl. 
Kosmogonie) und die astronomischen Instrumente behandeln. Mégen auch diese 
eine ebenso erfreuliche Vermehrung der astrophysikalischen Literatur werden wie 
der vorliegende erste Band. M. WALDMEIER 


Theorie schallnaher Stromungen. Von K. G. GUDERLEy (Springer-Ver- 
lag, Berlin 1957). 376 S., 125 Abb.; DM 42.-. 

Gasstr6mungen verhalten sich grundsatzlich verschieden, je nachdem die 
Str6mungsgeschwindigkeit kleiner oder grosser als die Schallgeschwindigkeit ist. 
Wenn auch die Theorie der reinen reibungslosen Unterschall- sowie der reinen 
Uberschallstromung zu einem gewissen Abschluss gebracht werden konnte, so 
bietet die Behandlung eines Str6mungsfeldes, in welchem gleichzeitig Unterschall- 
wie Uberschallgebiete auftreten, grosse Schwierigkeiten, die zur Hauptsache vom 
gemischten Charakter der Differentialgleichungen herriihren. Es ist daher ver- 
standlich, dass ein Buch iiber schallnahe Strémungen einen stark mathematisch 
betonten Charakter aufweist. Das Verdienst des Verfassers besteht darin, eine 
ziemlich liickenlose Zusammenfassung der weit zerstreuten Arbeiten zu bringen, 
in welchen er neben der mathematischen Darstellung immer wieder den physi- 
kalischen Vorgang auffiihrt. Leider enthalt das Buch sozusagen keine experimen- 
tellen Resultate, und wenn mehrmals auf die gute Ubereinstimmung zwischen 
Theorie und Wirklichkeit hingewiesen wird, so bleibt der Verfasser den Beweis 
schuldig. Dem mathematisch interessierten Aerodynamiker wird dieses Buch gute 
Dienste leisten. P. DE HALLER 


The Fundamental Principles of Quantum Mechanics. Von E.C. 
KeEmMBLE (Dover Publications, Inc., London und New York 1958). 611 S.; $2.95. 

Quantum Electrodynamics. Herausgegeben von J. SCHWINGER (Dover 
Publications, Inc., London and New York 1958). 424 S.; $2.45 

Die Tatsache, dass unser Zeitalter neben andern Eigenschaften durch eine 
Uberschwemmung mit billigem, bedrucktem Papier gekennzeichnet ist, hat auch 
ihre guten Seiten. Davon zeugen die beiden vorliegenden Bande. 

Der erste Band macht die bekannte griindliche Einfiihrung von E. C. KEMBLE 


in die Prinzipien der Quantentheorie jedermann zuganglich. Das Werk selbst hat 
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sich in den 20 Jahren seines Bestehens gentigend eingebiirgert und braucht kaum ~ 
noch eine Empfehlung. ' 

Von anderer Art sind die von J. SCHWINGER herausgegebenen 34 Arbeiten 
iiber Quantenelektrodynamik. Hier werden viele der wichtigsten Vorarbeiten zur 
modernen Quantenelektrodynamik (angefangen bei P. A. M. Drracs erstmaliger 
Quantisierung des Maxwell-Feldes im Jahre 1927) zusammengestellt. Darauf 
folgen, was sehr zu begriissen ist, zwei experimentelle Arbeiten, wovon die eine 
von W. E. Lams und R. C. RETHERFORD unmittelbar den Ansporn zur weiteren 
Entwicklung bildete. Daran schliessen sich die grundlegenden Arbeiten von 
J. Scuwincer, R. P. Feynman und F, J. Dyson an. Den Abschluss bildet eine 
Arbeit von G. KALLEN mit dem Resultat, dass nicht alle Renormalisierungskon- 
stanten endlich sein konnen. 

Die Einleitung von J. ScHWINGER gibt auf wenigen Seiten eine angenehm 
lesbare Einfiihrung in die Entwicklung der Quantenelektrodynamik. 

Natiirlich wird eine solche Auswahl von Artikeln den persdnlichen Stil des 
Herausgebers widerspiegeln. Sie lasst sich daher kritisieren. Doch ist ihr Nutzen 
und ihre Wiinschbarkeit so gross, dass jegliche Kritik nebensachlich wird. 

Es ist gut, dass jeder interessierte Student nun die Mdglichkeit hat, die 
Originalarbeiten selbst, die oft in schwer zugadnglichen Zeitschriften zu suchen 
waren, nicht nur zu lesen, sondern auch zu besitzen. R. Jost 


Physics and Mathematics. Band 2: Progress in Nuclear Energy. Von 
D. J. Hucues, J. E. SanpERs und J. HorowitTu (Pergamon Press, London 1958). 
S75) Set 90S: 4 
Der zweite Band in der Biicherserie «Physics and Mathematics» ist als Fort-_ 
setzung von Band 1 der Reaktor- und Neutronenphysik gewidmet. Nach einem © 
Uberblick iiber die experimentellen Daten, die simtliche Reaktionsméglichkeiten 
von schnellen Neutronen mit den schweren Kernen Thorium, Uran und Pluto- 
nium umfassen, wird ebenfalls zusammenfassend iiber die thermischen Einfangs- — 
querschnitte bei schweren Kernen unter Beriicksichtigung der Resonanzeffekte 
referiert. In die Neutronenphysik gehGren eindeutig auch die Aufsatze iiber die 
Standardisierung von Neutronenquellen und iiber Diffusionsmessungen mit 
gepulsten Neutronenquellen. Dazu kénnte man auch die Arbeit tiber die Ausmes- 
sung und Theorie der Reaktorspektren rechnen. Betrachtungen iiber den in — 
thermischen Reaktoren wichtigen Resonanzeinfang wahrend der Verlangsamung 
gehéren nun eindeutig in die Reaktorphysik. Methoden fiir die Berechnung von ~ 
heterogenen Reaktoren und das mehr in die angewandte Mathematik weisende ~ 
Kapitel tiber eine Monte-Carlo-Methode fiir Transportprobleme beschliessen das — 
Buch. Es handelt sich hier um ein typisches Spezialistenbuch, das vor allem durch 
brillante Zusammenfassungen einen grossen Wert bekommt. Allerdings sind in — 
einer Ausgabe so divergierende Stoffgebiete behandelt, dass selbst ein Spezialist 
nur einzelne Kapitel herausgreifen wird. Es ware deshalb zu begriissen, wenn in 
einer weiteren Buchfolge eine homogenere Zusammenstellung der Aufsatze er-_ 
reicht werden kénnte. P. STOLL ~ 
. | 
Fundamentals of Gas Dynamics. Band 3: H igh Speed Aerodynamics | 
and Jet Propulsion. Herausgeber: H. W. Emmons (Princeton University Press, | 
Princeton 1958). 749 S., 313 Abb.; $20.00. ; 
Das vorliegende Werk ist der dritte Band des von TH. von K&rmdn, H. L.7 
DRYDEN und H. S. Tayvtor herausgegebenen Handbuches iiber Aerodynamik der ; 
hohen Geschwindigkeiten und Strahlantrieb. Diesem Rahmen entsprechend wird die 7 
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Gasdynamik im weitesten Sinne behandelt.“Die Autoren der verschiedenen 
Abschnitte sind alle bekannte Spezialisten ihres-Faches, und ihre Beitrage be- 
stehen nicht nur aus einer Zusammenfassung der bisherigen Arbeiten, sondern 
enthalten wichtige eigene Betrachtungen, die dem Buche ein besonderes Geprage 
erteilen. Die acht Kapitel sind: 

A. H.S. Tsien: Allgemeine Gleichungen dey Gasdynamik. 

B. Luict Crocco: Eindimensionale stationdve Gasdynamik. Austfiihrliche 
Behandlung der Str6mungsvorginge in Leitungen und Diisen, mit stationdren 
Stéssen, Reibung, Verbrennung und Warmeiibergang. Auch der Einfluss einer 
Abweichung vom idealen Gas wird untersucht. 

C. ARTHUR Kantrowi1z: Eindimensionale instationévre Gasdynamik. Dieser 
Abschnitt behandelt die Stosswellen im idealen Gas, unter spezieller Beriicksich- 
tigung der zylindrischen und kugelsymmetrischen Explosion und Implosion. 

D. Wattace D. Hayes: Allgemeine theoretische Grundlagen der gasdynami- 
schen Stésse. 

E. H. PoLracuek und R. SEEGER: Gegenseitige Beeinflussung von Stosswellen. 

F. H. G. STEvVER: Kondensationsvorginge bei hohen Geschwindigkeiten. 

G. TH.von KArmAn, H.W. Emmons, R. S. TanxkIn und G.I. Taytor: Gas- 
dynamtk der Verbrennung und Detonation. 

H. S.A. ScHaar und P. CHAMBRE: Strémungen im hohen Vakuum. 

Das ausserordentlich inhaltsreiche Buch kann dem Fachmann und dem Stu- 
denten bestens empfohlen werden. P. DE HALLER 


Jahrbuch der Wissenschaftlichen Gesellschaft fiir Luftfahrt e. V. 
(WGL) 1957 (Friedr. Vieweg & Sohn, Braunschweig 1958). 504 S., 623 Abb.; 
DM 58.-. 

Den Hauptinhalt dieses Jahrbuches bilden die Vortrage der vom 9. bis 
12. April 1957 in Essen (Ruhrgebiet) gemeinsam von der WGL und der Deut- 
schen Versuchsanstalt fiir Luftfahrt (DVL) durchgefiihrten Tagung. Die Vortrage 
behandeln die Themata: 

Entwicklungsstand der Turbinentriebwerke — Probleme und Ergebnisse der 
angewandten gasdynamischen Forschung — Der Uberdruckwindkanal in Aachen 
— Ein schnellanzeigendes elektronisches Multimanometer — Vereinfachte Bestim- 
mung der Druckverteilung an Fliigeln beliebiger Spannweite — Str6mung an 
K6rpern nicht mehr kleiner Streckung in auftriebsloser linearer Unter- und Uber- 
schallstrsmung — Auftrieb und Widerstand langer Ringfliigel in Uberschallstré- 
mung — Stromungsuntersuchungen an rotierenden Schaufelgittern — Akustische 
Steuerung der turbulenten Anfachung im Freistrahl — Beitrag zur Frage des 
Lufteinlasses bei Strahltriebwerken fiir Uberschallgeschwindigkeiten — Der Ein- 
fluss der Fluggeschwindigkeit auf die Wirtschaftlichkeit von Durch- und Aus- 
strémtriebwerk — Die Wirkung der Schleierkiihlung bei hocherhitzten Bauteilen — 
Zur Frage der Fahrwerkslastannahmen — Ortungsgenauigkeit von Funknaviga-_ 
tionsverfahren — Uber das Springen von Funkpeilungen — Biologische Strahlen- 
wirkung — Beitrage zu Nachwuchs- und Ausbildungsfragen aus den Erfahrungen 
bei der Flugzeugfiihrerauslese — Uber die Berechnung des induzierten Geschwindig- 
keitsfeldes von Tragfliigeln — Untersuchungen zur Aerodynamik der Fliigel- 
Rumpf-Anordnungen — Str6mungsuntersuchungen an Tandem-Schaufelgittern - 
Uberschallstromung um Rumpf-Fliigel-Anordnungen — Grenzschichteffekte beim 


- Rohrwindkanal — Grenzschichtablésung und Totwasserbildung — Aeronautical 


1% Vir 


Studies in the Aeroballistics Range — Uber den Entwicklungsstand der Freiflug- 
anlagen — Beitrag zur Frage nach den Flugeigenschaften der Segelflugzeuge — 
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Grenzschichtmessungen im Segelflugzeug — Die schrag angeblasene Luftschraube 
und andere Fliigel-Triebwerk-Systeme im Langsamflug — Bemerkungen zur 
Schraganblasung am Rotorblatt eines Hubschraubers — On the Response of 
Structures Having Multiple Random Inputs — Durch die Fliigel verursachte 
Spannungserhéhung in Mittelscheiben axialer Turbolaufer — Ein Beitrag zur 
Statik des Ringfliigels — Design of ‘Passenger Tie-Down’. Some Factors for 
Consideration in the Crash-Survival Design of Passenger Seats in Transport 
Aircraft — Analog-Rechengerate und Simulatoren in der Luftfahrt — Uber einige 
neuere Ergebnisse der Regelungstheorie und ihre Bedeutung fiir die Lagenrege- 
lungen von Flugzeugen — Dreiachsen-Flugregelung fiir Hochleistungsflugzeuge 
mit Integrations-Wendekreiseln als Hauptrichtgeber — A High Resolution 
Hydraulic Servo — Sichtbarmachung von Flugzeugbordgeraten — Uber gasdyna- 
mische Mischungsvorgange; Schuberhéhung durch Strahlbeimischung — Ange- 
naherte Berechnung der mittleren Zustandsfelder in ebenen oder runden turbu- 


lenten Diftusionsflammen und dissoziierenden freien Gasstrahlen — Fliegende Priif- — 


stinde fiir Strahlturbinen — Probleme bei Versuchen an Ringbrennkammern — 
Turbinen-Flugtriebwerke, Leistungen und Gewichte — Die Herstellung von Blech- 
schaufeln fiir Gasturbinen und Kompressoren — S. A. P.—Der hochwarmfeste Alu- 
minium-Sinterwerkstoff — W&armeiibergang und Grenzschichtausbildung bei 
Schwitzbaustoffen — Die Arbeiten des Flugmeteorologischen Instituts der Deut- 
schen Forschungsanstalt fiir Segelflug (DFS). Freie und zweckgebundene For- 
schung — Neue Untersuchungen iiber die Eisteilchenbildung in unterkihlten 
Wolken — Wolkenhéhenmessungen bei Tage mit laufender Registrierung. 

Das Jahrbuch stellt somit wieder eine wertvolle Sammlung dar aus allen heute 
aktuellen Gebieten des Flugwesens: Aerodynamik, Flugmechanik, Triebwerks- 


forschung, Luftfahrzeugfiihrung, Werkstoffe und Festigkeit, Hubschrauber, ; 


Funknavigation, Meteorologie und Medizin. 
Den Schluss des auch drucktechnisch sehr gut ausgefiihrten Bandes bilden 


Angaben tiber die Tatigkeit der WGL im Jahre 1957, und mehrere Nachrufe ~ 


bringen das Wirken bedeutender Personlichkeiten in Erinnerung, die in der vom 
Jahrbuch erfassten Zeitperiode zum Teil in hohem Alter, zum Teil aber auch nur 
allzu jung verstorben sind. G. DATWYLER 


se 


Dislocations and Mechanical Properties of Cristals. Von J.C. FISHER, © 


W.C. Jonnston, R. THomson und T. VREELAND. Eine internationale Konferenz 


in Lake Placid vom September 1956 (John Wiley & Sons, New York 1957). 


634 S., 347 Fig.; $15.00. 


Die Festigkeitseigenschaften der Kristalle sind um Gréssenordnungen kleiner | 
als die theoretisch fiir Idealkristalle berechneten. Durch die Einfiihrung hypo- 


thetischer Gitterstérungen von linearer Gestalt, der sogenannten «Versetzungen», 
gelang es 1934 unabhangig voneinander G.I. Taytor, E. OROowAN und M. Po- 
LANYI, diese Diskrepanz zu beseitigen. Es folgten dann 20 Jahre, wahrend denen 
immer mehr Effekte der Kristallplastizitat und -festigkeit, aber auch zahlreiche 


ao 


andere Kristallphanomene (zum Beispiel gewisse Wachstumserscheinungen, be-_ 


stimmte Polymorphien usw.) sich miihelos vermittels der Versetzungstheorie er- 
klaren liessen, ohne dass aber ein direkter Beweis fiir das Vorhandensein von 
Versetzungen erbracht worden ware. In den letzten Jahren sind nun mit ganz 
verschiedenen Methoden Versetzungen in Kristallen sichtbar gemacht worden, so 
dass an ihrer Existenz nicht mehr gezweifelt werden kann. Anhand von 42 Bei- 
tragen gibt das vorliegende Buch einen Einblick in den heutigen Stand der Ver- 
setzungstheorie, ihre Anwendungen in der Festigkeitslehre der Kristalle und 


son Sonstes limeenetinedniaaaiaieneciaie 


Vol. X, 1959 Buchbesprechungen — Book Reviews — Notices bibliographiques 329 


weitere Probleme der Kristallkunde. Besondere Aufmerksamkeit ist dabei fol- 
genden Problemen gewidmet: direkter Nachweis‘von Versetzungen, Deformation 
reiner Einkristalle, Verfestigung und Erholung, Legierungswirkung auf die Streck- 
grenze, Dampfung und Ermiidung, Strahlenschiden, Whiskerkristalle (Haar- 
kristalle, welche sehr arm an Versetzungen sind und daher eine besonders hohe 
Festigkeit besitzen). Zahlreiche, zum Teil ausserordentlich gute Abbildungen und 
die Mannigfaltigkeit der Gesichtspunkte machen dieses Symposium zu einem 
wertvollen Informationswerk iiber die Versetzungen. W. EPPRECHT 


Elementary Introduction to Molecular Spectra. Von Borcr Bak 
(North Holland Publishing Company, Amsterdam 1954). 125 S., 35 Fig.; hfl.9—. 

Die Einfiihrung Baxs in die Molekularspektroskopie zeigt sich uns als kleines, 
aber sehr ansprechendes Buch von nur wenig mehr als 100 Seiten, das aber alles 
Wesentliche enthalt, das ein Nichtspektroskopiker wissen sollte, der sich des 
Verfahrens der modernen Molekularspektroskopie bedienen will, um Aufschluss 
iiber den strukturellen Aufbau der Stoffe zu gewinnen. 

Das Biichlein ist ganz auf die Bediirfnisse des Experimentes zugeschnitten. 
In den beiden ersten Kapiteln, die die Halfte des Umfanges beanspruchen, werden 
zunachst die Molekularspektren in ihren wesentlichen Ziigen erlautert, soweit 
dies fiir ein Verstehen der grundlegenden, im Molekiil sich anspielenden Vorgange 
und ihrer Wechselwirkungen notwendig ist. Hierauf folgen in den Kapiteln 3, 4 
bzw. 5 die Einzelbearbeitungen der Spektren des Mikrowellen-, des infraroten 
bzw. des sichtbaren und ultravioletten Frequenzgebietes. Zu Beginn dieser drei 
Kapitel findet sich je ein Hinweis iiber die Moglichkeiten, die das entsprechende 
Frequenzgebiet fiir die Erkundung des molekularen Aufbaues bietet. 

Gerade diese kurzen Hinweise bilden zusammen mit den beigeordneten, in 
knapper Darstellung gegebenen theoretischen Erlauterungen das wertvolle Ge- 
prage der Schrift. In erster Linie wendet sich die Monographie an Biologen, 
Chemiker und Ing.-Chemiker, die die Méglichkeiten der Molekularspektroskopie 
fiir ihre Forschungen zunutze ziehen wollen; zumal heute die notwendigen 
Spektralgerate in zuganglicher Form und einfacher Handhabung zur Verfiigung 
stehen. R. SANGER 


Der lichtelektrische Effekt und seine Anwendungen. Von H. Simon 
und R. SUHRMANN (Springer-Verlag, Berlin 1958). 747 S., 599 Abb.; DM 97.50. ° 

Dieses Buch erscheint als neue Auflage des 1932 unter dem Titel Lichtelektrische 
Zellen und ihre Anwendungen publizierten Buches. 

Das Hauptgewicht der Forschung auf dem Gebiet der Festkérperphysik lag in 
den letzten 25 Jahren bei den Halbleitern. Die erzielten Fortschritte machten eine 
Revision verschiedener Kapitel notwendig, so vor allem jener tiber den innern 
Photoeffekt, wo einige hervorragende neue Photoleiter, wie Pbs, CdSe tigsiuswe, 
zu beriicksichtigen waren. Uber die Si-Zelle wird leider noch nichts gesagt. 

Aus diesem umfangreichen Buch werden wir nur die wichtigsten Kapitel er- 
wahnen. Nach einer kurzen phanomenologischen Beschreibung des ausseren, inne- 
ren und Sperrschicht-Photoeffektes wird in Kapitel II von R. SUHRMANN der aus- 
sere lichtelektrische Effekt im Detail betrachtet. Der Physiker findet dort eine 
genaue Beschreibung der Theorien von FowLer und Du BRIDGE liber die Tempe- 
raturabhangigkeit des Photoeffektes. Weiter werden die Einfliisse der Oberflachen- 


- schichten und der Kristallstruktur auf die Energieverteilung der austretenden 


Elektronen experimentell und theoretisch untersucht. Die Eigenschaften der tb- 
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lichen Photozellen (K, Cs, K,Sb, Cs,Sb, K-Diphenylather, ...) sind ‘tabelliert und 
graphisch dargestellt. “y 

Das Kapitel III, Innere lichtelektrische Effekte, ist von einem. Spezialisten, 
K. W. Boer, geschrieben worden. BéER beginnt wirklich am Anfang und ent- 
wickelt die ganze Bandertheorie des festen Kérpers, ausgehend vom periodischen 


Potential, was ihm erlaubt, die Photoleitung und den Sperrschichtphotoeffekt zu - 


erklaren. 

In den Kapiteln IV und V werden von H. Simon die Herstellungsmethoden fiir 
Photozellen, sowohl mit 4usserem wie mit innerem Effekt, von Grund auf ange- 
geben. Wir finden die Beschreibung sogar fast zu griindlich, denn der Autor hatte 
vom Leser sicher mehr Kenntnisse voraussetzen k6nnen. Nicht nur werden die 
Herstellungsverfahren der Zellen selber erlautert, sondern auch alle dazu gebrauch- 
ten Apparate beschrieben. Nur ein Beispiel: Der Aufbau des Pumpstandes von der 
Vorvakuumpumpe bis zur Getterpumpe, iiber den Lecksucher, das Alpert-Mano- 
meter usw., alles ist der Schilderung zu entnehmen. 

Die Kapitel VI bis XII stammen von F. Ecxart, W. LEo, R. SUHRMANN und 
H. Simon und sind speziellen Effekten und Anwendungen gewidmet: Sekundar- 
elektronen-Verstarkung, Methoden der lichtelektrischen Messungen, Anwendungen 
der Photozellen in der Photometrie, als Bildwandler, in der Fernsehtechnik, auf 
dem Tonfilm usw. 

Wichtig und besonders wertvoll sind die zahlreichen Literaturangaben (bis und 


mit 1955) am Ende jedes Kapitels, die vielen Tabellen und detaillierten Zeichnungen ~ 


von sehr vielen Rohren und Zellen. P. JUNOD 


Introduction to the Physics of Many-Body Systems. Von D. TER HAaR © 


(Interscience Publishers Inc., New York 1958). 127 S.; $ 1.95. 
Das vorliegende kleine Buch hat als Thema den in letzter Zeit besonders im 


Vordergrund des Interesses stehenden Fragenkreis der Vielkérper-Probleme. Zu — 
‘diesen Problemen gehoéren unter anderem so unterschiedliche Dinge wie die Eigen- 


schaften der Kernmaterie, das Verhalten eines Bose-Gases mit abstossenden Kraf- 
ten als Modell fiir die merkwiirdigen Eigenschaften von He? bei tiefen Temperaturen 
und die von THomaAs und Fermi erfundene statistische Methode zur Behandlung 
schwerer Atome. 

Der Autor gibt eine angenehm lesbare Ubersicht iiber das Gebiet, die eher die 
mehr intuitiven Methoden und weniger die anspruchsvollen konsequenten mathe- 
matischen Durchfiihrungen betont. Etwas anderes ware ohne eine starke Vergrés- 
serung des Umfanges auch kaum méglich. 


Imposant ist das Literaturverzeichnis. Vollstandigkeit konnte aber auch hier 


nicht angestrebt werden. So vermisst man etwa die Arbeiten der Schule von L. vAN _ 


Hove. 


Als Einfiihrung kann das kleine Werk empfohlen werden. R. Jost@ 


Ergebnisse der Hochvakuumtechnik und der Physik diinner Schich- 


ten. Herausgegeben von M. AuwirTER (Wissenschaftliche Verlagsgesellschaft — 


mbH, Stuttgart 1957). 282 S., 209 Abb:; DM 52.-. 


Das vorliegende, von Dr. M. AuwARTER herausgegebene Buch umfasst 19 Son- _ 


derbeitrage, welche aus der Feder fiihrender Wissenschafter stammen. Bei der 


Zusammenstellung dieser Beitrage ist eine Systematik weder verfolgt noch gewollt. 


Jeder Beitrag ist entweder eine Originalmitteilung oder ein kritischer Zusammen- 
punesbaneht iiber ein bestimmtes Teilgebiet. Unter letzteren geh6ren die zwei 
Beitrage: Die Aufdampfung von diinnen Schichten im Hochvakuum von W. EsPpE 
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_ und Grundlagen und Grenzen vakuumtechnischer A rbeiisverfahven in dey Metallurgie 
von O. WINKLER. Beide Beitrage geben eine knappé und prizise Darstellung des 
behandelten Stoffes. Fiinf Arbeiten theoretischer und experimenteller Natur sind 

-dem Problem der Interferenzfilter und der Glasentspiegelungsschichten gewidmet. 
Besonders hervorzuheben ist ein Beitrag von M. AUWARTER, R. HarFer und P. 
RHEINBERGER tiber die Schichtdickenmessung mit Hilfe von optischen Interferenz- 
filtern. Es wird hierbei nicht nur ein neuer Weg zur Messung von Schichtdicken von 
der Grésse von einigen Angstr6m angegeben, sondern auch ein Verfahren ent- 
wickelt, mit dessen Hilfe Informationen in Form von kleinsten Massendicke- 
Unterschieden in einer Aufdampfschicht gespeichert, konserviert und sichtbar 
gemacht werden kénnen, ahnlich wie bei einer Photoplatte. Die Brauchbarkeit dieser 
Art des Photographierens wird in einem Beitrag von T. Rur und O. WINKLER iiber 
die Untersuchung der Riickstrémung von Ol-Diffusionspumpen mit Hilfe der Loch- 
kamera sehr anschaulich unter Beweis gestellt. Eine Reihe von Arbeiten beschaftigen 
sich mit den verschiedenen Aspekten der diinnen Aufdampfschichten, darunter eine 
sehr aufschlussreiche Abhandlung von G. Hass und A. F. TuRNER iiber die opti- 
schen Aufdampfschichten fiir das Infrarot-Gebiet. Auch die Abhandlung von 
H. Maver, H. THomas und A. W. ELBeEt iiber die Messung der Lichtabsorption 

“und der Photoemission an monoatomaren Kalium-Aufdampfschichten sei noch 
erwahnt. 

Ohne Zweifel weist das Buch von AUWARTER eine Sammlung von sehr wertvollen 
Abhandlungen auf, und das Studium desselben wird allen Lesern, die sich fiir die 
Physik diinner Schichten und fiir die Hochvakuumtechnik interessieren, neue An- 
regungen bringen. E. Bas 


Electronic Digital Computers. Von Franz L. Att (Academic Press, Inc., 
New York 1958). 336 S., 8 Fig.; $ 10.-. 

‘Die Nachfrage nach Information tiber das Gebiet der elektronischen digitalen 
Rechenmaschinen wachst gleichermassen wie die Produktion und der Bedarf an 
solchen Rechenautomaten. Der Autor will mit dem vorliegenden Buch in tiber- 
blickartiger Form eine Einfiihrung in die wesentlichen Aspekte dieses Gebietes 
geben. 

Die Einleitung, das erste von 5 Kapiteln, enthalt die Beschreibung der elemen- 

- tarsten Begriffe und des Aufbaus des Buches. Im Kapitel 2 werden dem Leser in 
einem interessanten historischen Uberblick die wichtigsten, vor allem amerikani- 
schen Maschinen vorgestellt und ihre charakteristischen technischen und organisa- 
torischen Merkmale erwahnt. Von den verschiedenen Teilen einer Maschine werden 
lediglich ihre Funktionen beschrieben, so dass fiir das Verstandnis keine technischen 
Kenntnisse notwendig sind. Alle Maschineneigenschaften, die das Programmieren 
betreffen, kommen im Kapitel 3 zur Sprache. Es werden Programmierungsbeispiele 
fiir verschieden organisierte Maschinen behandelt und ganz kurz das Aufstellen 

- von Flussdiagrammen, Unterprogrammen und die Verwendung der Maschine zu 

_ Programmierungszwecken gestreift. 

Das Kapitel 4 iiber die Analyse der Probleme darf nicht mit einem Text tiber nu- 

- merische Verfahren verglichen werden. Es will nur eine Zusammenstellung und Er- 

__ klarung der landlaufig beniitzten numerischen Methoden geben fiir die Berechnung 

~ expliziter Funktionen, das Differenzieren und Integrieren, die Lésung von gewohn- 
- lichen und partiellen Differentialgleichungen, algebraische Gleichungen, Eigenwert- 


_ Theorie formal eingefiihrt. Im folgenden Kapitel werden die Booleschen Funktionen 
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; 
und anderen Gebieten, wobei die wesentlichen Eigenschaften der je einem bestimm- — 
ten Zweck am besten angepassten Maschinen abgeschatzt werden. Am Ende folgt | 
eine Betrachtung tiber die Organisation eines Rechenbetriebes. Umfangreiche Lite- 


’ raturangaben vermitteln den Zugang zur Spezialliteratur. 


FRANZ L. Att wendet sich mit seinem Buch in erster Linie an Leute, die genotigt — 
sind, sich in das Gebiet der elektronischen Rechenmaschinen einarbeiten zu missen, 
und die die Prinzipien, welche dem Entwurf und dem Funktionieren von Rechen- — 
maschinen zugrunde liegen, kennenlernen méchten, also an Physiker, Chemiker, 
Ingenieure und ahnlich Beschaftigte, die Maschinen zur Loésung ihrer Probleme 
verwenden kénnen. Er will dem Beniitzer einer Maschine das Werkzeug in die Hand 
geben, das ihm die Verstandigung mit dem Rechenmaschinen-Spezialisten erlaubt — 
und damit die Méglichkeit gibt, in kurzer Zeit auf irgendeiner Maschine Probleme ~ 
selbstandig zu lésen. Das Buch entspricht einem wirklichen Bediirfnis, und das vom _ 
Autor gesteckte Ziel wurde in jeder Hinsicht erreicht. A. SCHAI 


Logical Design of Digital Computers. Von MontTGoMERyY PHISTER, JR. : 
(John Wiley and Sons, New York 1958). 408 S., 128 Fig., 147 Tab.; $ 10.50. 1 
Der Autor teilt die Entwurfsarbeiten an einer Rechenmaschine in 3 Gebiete ein: — 
Den Entwurf des ganzen Systems, den logischen Entwurf und den Entwurf der 


"7 


- Grundschaltungen. Das vorliegende Buch behandelt ausschliesslich nur den logi- | 


schen Entwurf von digitalen Rechengeraten. In 12 Kapiteln werden eingehend die 
Methoden und Techniken beschrieben, die von verschiedenen Fachleuten angewen- 
det werden. Das Buch vermittelt dem Leser die notwendige Information und das 
Werkzeug, um den vollstandigen logischen Entwurf einer Rechenmaschine aus- 
fiihren zu k6nnen; spezielle Kenntnisse werden keine vorausgesetzt. Mit Ausnahme 
von folgenden Einschrankungen wird das behandelte Gebiet recht tiefgehend dar-— 
gestellt: Die Untersuchungen werden auf synchrone Schaltungen beschrankt, das © 
bewegliche Komma kommt nicht zur Sprache, und Division und Parallelmaschinen — 
werden nur ganz kurz gestreift. § 

Wegen der grossen Bedeutung der Booleschen Algebra wird diesem Gegenstand — 
ein ganzes Kapitel gewidmet. Die Boolesche Algebra wird darin als eine deduktive — 
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mit verschiedenen Hilfsmitteln vereinfacht (Minimisierung): Mit der Methode von | 
Quine, der Harvard-Methode und der Veitch-Diagramm-Methode. Die Anwendung — 
der Booleschen Algebra auf die Behandlung von Flipflop-Speicherelementen fiihrt — 


‘auf sogenannte «difference equations». Es handelt sich dabei im wesentlichen um die _ 


logische Untersuchung von Elementen, deren Zustand von friitheren Zeitintervallen 
abhangt (Kapitel 5). Das Huffman-Moore-Modell einer digitalen Rechenmasohindl 
wird in den Kapiteln 6 und 9 erklart und mit Erfolg auf verschiedene Probleme © 
angewendet. Dieses Modell ist nichts anderes als eine etwas vereinfachte Turing- _ 
Maschine. Schlussendlich wird in Kapitel 11 das Leitwerk einer Maschine sowie das _ 
se beim logischen Entwurf einer vollstandigen digitalen Maschine be- _ 
schrieben. 


Von besonderem Wert dieses Buches seien die Abschnitte tiber die ee 


©. 


gramm-Methode, die Aufstellung von «difference equations» fiir Speicherelemente _ 
und das Huffman-Moore-Modell hervorgehoben. Dem Aufbau des Buches wurde _ 
die Form eines 1—2semestrigen Kurses gegeben. Dementsprechend wurden zwischen | 
die einzelnen Kapitel Ubungsaufgaben eingeflochten sowie die wichtigsten, das — 
jeweilige Kapitel betreffenden Literaturangaben. | 
; Fiir das Gebiet des logischen Entwurfes von digitalen Gerdten stellt das vor-_ 
liegende Buch einen wertvollen Beitrag dar. ; A. ScHA 
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«Der Verfasser gibt hier ein neues Kriterium fiir die Vertraglichkeit eines Geschwin- 
digkeitsfeldes mit einem gegebenen Gleitlinienfeld. - Das Buch bildet eine sehr 
wertvolle Erganzung zu dem vom gleichen Verfasser stammenden Lehrbuch der 
Plastizitatstheorie. » _~ Internationale Mathematische Nachrichten, 1956 


«Cet ouvrage, d’un niveau élevé, constitue pour les ingénieurs spécialisés et pour 
les laboratoires d’essais et de recherches un exposé d’ensemble des questions de 
plasticité qui les aidera 4 mieux percevoir les possibilités d’application de cette 
propriété des matériaux. » Travaux, 1956 


«Das angenehm zu lesende Buch erscheint als Einfiihrung in das neuerdings immer 
wichtiger werdende Gebiet der Plastizitatstheorie besonders geeignet, weil es, zum 
Teil unter Zuhilfenahme neuer Modellvorstellungen, gerade die Grundvorstellungen 
so ausfiihrlich entwickelt.» ; VDI-Zeitschrifi, 1956 


«Es ist sehr zu begriissen, dass der Verfasser dieses einfiihrende Buch in die Plasti- 
zitatstheorie herausgegeben hat, gibt es doch eine dusserst anschauliche Darstellung 
der Grundlagen und der fiir die praktische Anwendung wichtigen abschatzenden 
Verfahren. Heute, wo auch der Stahlbau sich bemiihen muss, plastische Verfor- 
mungen zuzulassen, wird das Buch grossem Interesse begegnen.» 

Zeitschrift fir Flugwissenschaften, 1956 


Zu beziehen durch Ihre Buchhandlung — Obtainable from your bookseller 
Commandes & votre libraire 


BIRKHAUSER VERLAG: BASEL UND STUTTGART 


2 Ly CAS 


ARCHIV DER MATHEMATIK 
Archives of Mathematics — Archives Mathématiques 


Herausgegeben in Verbindung mit dem Mathematischen Forschungsinstitut in 
Oberwolfach von H. KngsEr - W. Stss f 


Beirat: G. Bot, E. Bomprani, J. DiEUDONNE, CH..EHRESMANN, W. FELLER, ea 
H. Gortter, H. Hapwicer, H. Horr, H. Kénic, S. MacLane, W. Macnus, 
_ T. NacEtz, Cur. Pauc, G. Pickert, K. REIDEMEISTER, P. ROQUETTE, J. A.ScHOU- | 
TEN, H. Serrert, E. SPERNER, E. STIEFEL, P. TURAN 

Schriftleitung: E. LAMPRECHT 


Jahresabonnement (6 Nummern): sFr./DM 66.—, Einzelnummer sFr. 14.— 
Annual subscription (6 issues): Sw. Fr./DM 66.-, single copy Sw. Fr. 14.— 
Abonnement annuel (6 numéros): Fr. s./DM 66.—, le numéro Fr. s. 14.— 


HELVETICA PHYSICA ACTA 


Redaktionskomitee — Comité de rédaction — Comitato di redazione 
P. HuBEr, Basel; J. Rosser, Neuchatel; A. MERCIER, Bern; A. PERRIER, Lausanne, 
P. SCHERRER, Ziirich; J. WEIGLE, Genéve 
' Redaktor: Prof. Dr. M. Frerz, Physikalische Anstalt der Universes Basel: 
Klingelbergstrasse 82 ; 


Jahresabonnement (6-8 Nummern pro Jahr) sFr./DM 64.—, Einzelhefte sFr./DM 14.— 
Annual subscription (6-8 issues a year) Sw. Fr./DM 64.-, single copy Sw. Fr. 14.— 
Abonnement annuel (6-8 numéros par an) Fr.s./DM 64.-, le numéro Fr. s. 14.— 


ELEMENTE DER MATHEMATIK 


" Revue de mathématiques élémentaires — Rivista di matematica elementare 


Zeitschrift zur Pflege der Mathematik 
und zur Férderung des mathematisch-physikalischen Unterricht 
Organ fiir den Verein Schweizerischer Mathematik- und Physiklehrer 


Redaktion: Dr. L. Locurr-Ernst, Direktor und Professor am Technikum, 
Nussbaumweg 4, Winterthur 

Dr. E. Trost, Professor am Technikum Winterthur, Basteiplatz 3, Ziirich 1 

Dr. P. Bucuner, Prof. an der Universitat, a. Rektor des Math.-Naturw. Cyptieaate 

Realpstrasse 71, Basel 


Jahresabonnement (6 Nummern): sFr./DM 14.-, Einzélnummer sFr. [DM 2.50 
Annual subscription (6 issues): Sw. Fr./DM 14.—, single copy Sw. Fr. [DM 2.50 
Abonnement annuel (6 numéros): Fr. s./DM 14.-, le numéro Fr. s. 2. 50 


Verlangen Sie kostenlose Probenummern 


‘ Birkhauser Verlag - Basel und Stuttgart 


